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La reprise d'un ouvrage en sous-œuvre imposait souvent 
des travaux lourds et coûteux. Parmi les techniques moins 
invasives mises au point par les spécilistes du secteur l'amé- 
ligration des sols par injections de résine expansive se révèle 
particulièrement efficace 

Employée pour consolider les sols des fondations d'un ou- 
vrage l'injection de résine expansive est une technique dé- 
sormais bien connue de l'ensemble des professionnels de la 
construction intervenant dans les travaux de préservation, 
de restauration ou de réhabilitation. 

Néanmoins — afin de mieux appréhender tout le potentiel 
de cette technique inventée voici près d'une vingtaine d'an- 
nées — les mécanismes de diffusion et d'interaction avec 
le sol de la résine qui y est infectée (résine polyuréthane) 
restaient encore à préciser. 

À cette fin, deux chercheurs de l'École polytechnique de 
Turin, Mario Manassero et Andrea Dominijanni, ont conduit 


Andre 


une étude théorique et expérimentale dont les résultats sont 
publiés ici; financée par la société Uretek, cette étude porte 
sur les effets de l'injection de résine expansive, dans les sols 
fins comme dans les sols grossiers. 

A la suite des récents progrès réalisés dans la modélisation: 
théorique et expérimentale des milieux poreux et particu- 
laires,les auteurs proposent notamment dans ce livre un pre- 
mier cadre théorique général pour la modélisation analyt 
et numérique de ce type d'injections; dans leurs recherches 
il se sont par ailleurs appuyés sur le procédé Uretek Deep 
Injections*, seul procédé d'amélioration des sols faisant l'ob- 
let d'un avis technique du CSTB (avis technique 3/15-796). 
C'est ainsi qu'après avoir étudié différentes configurations 
d'injection, les auteurs montrent comment réallserune opé- 
ration de consolidation à l'aide de la technologie Uretek Deep 
Injectians* dans le but d'atténuer les effets induits par une 
activité sismique. 


Deminijanni, chargé de recherche en géotechnique environnementale et enseignant, exerce à l'École polytechnique de Tu- 


rin. Spécialiste des phénomènes de transport et du comportement mécarique des sais ins, ilest l'auteur de nombreuses publications 
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Consolidation de sols et stabilisation 
d'une maison individuelle à l'aide du 
procédé URETEK Deep Injections® (avis 
technique du CSTB 3/15-796). 
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du procédé URETEK Deep _ injections" 
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Propriétés physiques 
et mécaniques 
des résines polyuréthanes 
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Les polymères d'uréthane, ou polyuréthanes, forment une grande famille de polymères 
qui se caractérisent par diverses propriétés physiques, notamment leur capacité d'expan- 
sion (gonflemenc). 


Résines polyuréthanes expansives 


La réaction chimique pour la synthèse des polyuréthanes a été découverte en 1947 par 
Bayer. Les résines polyuréthanes expansives sont produites lors d'une réaction exother- 
mique entre un polyol et un isocyanate, mélangés dans des proportions volumiques éta- 
lies en fonction de caractéristiques de production spécifiques. Au cours de la réaction, 
une grande quantité de dioxyde de carbone est générée, ce qui provoque l'expansion volu- 
métrique du mélange et la formation d'une structure spongieuse, dans laquelle les bulles 
de gaz sont piégées. La production de dioxyde de carbone implique nécessairement la pré- 
sence d'eau, laquelle réagir avec le groupe isocyanate; sans eau, il est nécessaire de recou- 
riräun agent d'expansion liquide, chimiquement inerte et avec un faible point d'ébullition, 
qui se vaporise en exploitant une partie de la chaleur de polymérisation. 


Dans un laps de temps très restreint (de quelques secondes à quelques minutes), le 
mélange durcit passant de l'érar liquide à l'éar solide. Le temps de réaction, qui dépend 
tiques de la résine et de l'utilisation, entre autres, de catalyseurs, est influencé 
par la température des composants mélangés. En conséquence, le contrôle de la tempéra- 
ture des composants permer d'accélérer ou de retarder le temps de réaction. 


des caractér 


La pression exercée pendant le gonflement er la densité finale de la résine dépendent de la 
capacité d'expansion que possède le gaz contenu dans Les bulles, avant le durcissement, La 
structure à «cellules fermées» de la résine expansive est représentée sur la figure 1.1, qui 
regroupe les images obtenues au microscope électronique par Buzzi et al. (2008). Les 
images reproduisent la structure microscopique de La résine utilisée par Uretek (Canteri, 
1998) dans le domaine géoechnique, dans des conditions de gonflement libre, auxquelles 
correspond une densité égale à 37 kg/m%. Dans ces conditions, le volume expansé est égal 
à trente fois le volume initial du mélange (Dei Svaldi et al, 2005; Buzzi et al, 2008; 
Buzz er al., 2010). 

La densité du mélange à l'érat liquide est égale à 1.070 kg/m? er est très proche de celle de 
l'eau (= 1000 kg/m). La résine expansée est, au contraire, caractérisée par une densité 
sensiblement plus faible, qui dépend de la pression de gonflement. 


De nombreux essais en laboratoire ont été effectués sur la résine Uretek afin de détermi- 
ner ses principales propriétés physiques er mécaniques (Favaretri er al., 2004). 
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Figure 1.1 Image obtenue au microscope électronique d'une résine polyuréthane Uretck expansée sans 
confinement (densité égale à 37 kg/mo. () Agrandissement x 100 () agrandissement x 200, 
Buzz e al, 2008). 


En particulier, la relation entre la pression de gonflement et la densité (ou poids spéci- 
que) a été étudiée en conditions œdomérriques, à l’aide d’un dispositif spécial qui permet 
d'injecter une résine à l'intérieur d’un cylindre métallique rigide, muni d'un piston. 
Immédiatement après l'injection, la résine commence son expansion, qui, en raison de la 
forte rigidiré du conteneur, ne peut se faire que dans le sens vertical. Lors de son expan- 
sion, le mélange déplace un piston vers le haut, qui, après une course de quelques centi- 
mètres, est bloqué par un élément transversal équipé d'un manomètre, La pression de 
gonflement à été déterminée comme celle devant être fournie par l'élément transversal 
pour empêcher le piston de se déplacer vers le haut. La dépendance de la pression de gon- 
flement vis-à-vis de la densité de la résine, obtenue expérimentalement, est représentée sur 
la figure 1.2. 
À partir des résultats expérimentaux, on à déduit la relation empirique suivante: 
P= expf 0,23: (y, 0,36)]-1 Un 

où: 
+ P= pression de gonflement, déterminée dans des conditions œdométriques (exprimée 

en MPa); 


= Y,r= poids par unité de volume (ou poids spécifique) de la résine expansée (exprimé 
en KN/mÿ). 


Résines polyuréthanes expansives 


"2 


Relation empirique 
10 P=exp 10,23: (y,-0,36)] - 1 L 


Pression de gonflement mexmale, P (MPa) 


Poids volumique de la résine, y 


Figure 1.2 Pression de gonflement de la résine Uretek en fonction du poids par unité de volume 
(données expérimentales de Favarett et al, 2004) 


Dans les essais effectués, les valeurs de pression de gonflement mesurées étaient comprises 
entre 0,40 et 10,00 MPa, dans le domaine étudié du poids par unité de volume (y,= 2.00 
à 10,00 KN/m). En outre, Buzz er al. (2008) ont rencontré des conditions de gonfle- 
ment libre de la résine, où pression de gonflement nulle, pour un poids spécifique 


36% , où g est l'accélération de la pesanteur (= 9,81 m/s). La rcla- 


tion empirique de l'équation [1.1] a été obtenue en minimisant l'erreur d’approximation 
des chiffres de Favaretti er al. (2004) et en imposant 4, = 0,36 pour P = 0 kPa 

D'après ces résultats il apparaît que la pression de gonflement diminue au fur et à mesure 
de la diminution de La densité et de l'augmentation de l'expansion. Cette propriété revêt 
une importance capitale pour les applications dans lesquelles les résines polyuréthanes 
sont injectées à l'intérieur du sol et, par conséquent, soumises à un effec de confinement. 
D'après l'équation [1.1], il est possible de déterminer le rapport entre Le volume expansé de 
la résine (volume de la résine finale V,) et le volume initial du mélange injecté à l'état 
liquide à l'intérieur du sol (volume de la résine initiale V, ). En effet, la masse de résine M, 
ne varie pratiquement pas lors du passage de l'état liquide initial à l'état final de solide 
expansé et, par conséquent, les volumes Ver V: peuvent être exprimés de la façon suivante: 


[2] 


U2b] 
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masse de la résine, inchangée lors du passage de l'étar liquide à l'état solide; 


+ v,=1070 8. g = 105 N° - poids spécifique de la résine à l'éar liquide; 
nm m 


ique de la résine expansée. 


des équations [1.1], [1.2a] et [1.2b], le rapport de gonflement, défini comme 
VlVin est égal à: 


1 
(7 036 in [l+ PP] 


[1.3] 


La relation entre Ver V obtenue avec cette procédure est illustrée sur la figure 1.3. 


Comme cela a déjà éré observé, si la relation de polymérisation a lieu dans un environne- 
ment exempe de confinement (air libre), la résine se solidife avec un poids volumique égal 
20,36 KN/m° avec un taux d'expansion volumétrique égal à environ 30 (Dei Svaldi et al, 
2005). 


La relation reproduite sur la figure 1.3 revêt une importance capitale pour l'évaluation 
des injections de résine dans les sols (Foti & Manassero, 2009). 
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Figure 1.3 Relation expérimentale entre Le rapport de gonflement et la pression de gonflement 
de la résine Uretek 
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2. Propriétés mécaniques et hydrauliques 
des résines polyuréthanes 


Les résines polyuréchanes, la suite des processus de gonflement er de durcissement, déve- 

loppent des caractéristiques mécaniques (à savoir: rigidité et résistance) très élevées, en 

fonction du poids spécifique Anal, ÿ,: 

Sur les résines Uretek, on a déterminé de manière expérimentale (Favaretti et al., 2004; 

Dei Svaldi er al, 2005; Foti & Manassero, 2009) : 

+ le module de rigidité (ou module d'Young), £, correspondant à une déformation de 
1% (figure 1.4); 

+ la résistance à la compression monoaxiale, 0, (figure 1.5). 
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Figure 1.4 Module de rigidité en fonction du poids par unité de volume final dela résine Uretek. 
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Figure 1.5 Résistance à la compression monoaxale en fonction du poids par unité de volume final 
de la résine Uretek 
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Ces deux paramètres ont été évalués en fonction du poids spécifique y,; dans un inter- 
valle de valeurs comprises entre 0,5 et 3,5 KN/m3. 

En ce qui concerne le module de rigidité, les éprouvettes ont permis d'identifier des 
valeurs comprises entre 10 er 80 MPa, comparables aux modules caractéristiques des rer. 
rains densifiés. 

Les résistances à la compression mesurées sont significativement élevées (0,25 à 6,50 MPa). 
À partir de celles valeurs, il est possible de calculer la résistance au cisaillement de la résine, 
supposée purement cohésive, au moyen de la relation 


t,= U4] 


2 


Il en résulte que la résistance au cisaillement varie entre 0,125 et 3,25 MPa, en fonction 
de la densité finale ou du degré d'expansion de La résine. 

Des essais d'imbibition dans l'eau ont également été effectués sur les résines Uretek. Ils 
montrent des variations de poids spécifiques comprises entre 2 et 3 % à une température 
de 23 °C. (Gabassi, 2011). 

Buzzi #r al. (2010) ont, en outre, effectué des essais de conductivité hydraulique à la fois 
sur les résines Urerek pures et sur des échantillons constitués d'un mélange sol-résine. 
Pour la densité minimale (0,36 kN/m°), ou gonflement maximal, les valeurs de condue- 
ivité hydraulique mesurées étaient comprises entre 1X10-? er 1X10-% m/s sur la résine 
pure et indiquaïent l'existence d’une porosité réduite interconnectée. Sur les échantillons 
de sol injectés, la conductivité hydraulique s'est révélée égale à environ 101 m/s. 

À partir de ces résultats, il est possible de conclure que la structure de la résine est presque 
entièrement constituée de pores fermés et imperméables à l'eau. La conductivité hydrau- 
lique extrêmement faible, mesurée sur des échantillons de sol injectés, est due à la pré- 
sence de micro-défauts aléatoires. Dans les conditions naturelles, le poids spécifique de la 
résine ne varie pas selon les conditions d'immersion dans l'eau. 

Ces résines peuvent donc être considérées comme presque imperméables et pratiquement 
étanches à l'eau. 


[3 _ Injection de résines Uretek dans les sols 


Les résines polyuréthanes peuvent être injectées dans les sols à diverses profondeurs par 
des forages de petits diamètres, ce qui évite de fragiliser les ouvrages de fondations et/ou 
Les dallages traités qui les surmontent. Une description précise de la technologie Uretek 
dénommée Deep Injections a été effectuée par Canteri (1998). Avec cette technologie, il 
est possible de réaliser des percements d'un diamètre inférieur à 30 mm (en général, de 12 
à 26 mm), avec un entraxe minimal compris entre 45 et 150 em. Des tubes d'acier sont 
insérés dans ces forages pour l'injection de la résine polyuréthane. Les interventions sont 
effectuées par une équipe de techniciens spécialisés, avec un camion-atelier équipé. Des 
camions de différentes dimensions sont disponibles selon le type d'intervention et selon 
les conditions d'accès au chantier Dans les cas exceptionnels, il est possible d'acheminer 
le matériel ec la résine par d'autres moyens de transport: trains, avions, bateaux ou télé- 
phériques. À l'intérieur du camion-atelier, en plus de divers outillages, se trouvent les 


LN 


Injection de résines Uretek dans les sols 


réservoirs des deux composants chimiques nécessaires à la production de la résine. Ces 
deux composants sont mélangés avant leur entrée dans le tube d'injection et ont un temps 
de réaction variable allant de quelques secondes à quelques minutes. L'injection est effec- 
tuée par un opérateur à l'aide d'un pistolet d'injection prévu à cet effet. 


Critères de conception 


[1 Introduction 


(Au cours des dernières années, les cechniques d'injection à l'aide de résines polyuréthanes 
se sont perfectionnées (Foti & Manassero, 2009); et aujourd'hui, elles sont employées 
dans de nombreuses applications (figure 2.1), notamment: 


+ le comblement et la stabilisation de cavités souterraines: 

+ le déplacement d'eau et la réduction de la conductivité hydraulique: 

+ l'allègement des charges pour des couches de rerrain situées sous Le volume de sols crai- 
rés, afin d'en limiter l'affaissement; 

+ la consolidation des sols, pour augmenter la rigidité er la résistance au cisaillement; 

+ le relevage des fondations et des dallages, pour la reprise des afaissements. 

Le présent chapitre décrit brièvement les critères de conception des interventions de 

consolidation effectuées à l'aide d'injections de résine polyuréthane expansive, en fonc- 

tion de leur finalité, 
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Figure 2.1 Application d'ijections de résine polyuréthane (adaptation d'après wunwuretekusa com) 
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Pour exécuter un projet de traitement adapté à l'aide d'injections de résine polyuréthane 
expansive, il est nécessaire de réaliser une étude détaillée du site qui tienne compte de sa 
structure géologique, des aspects géotechniques, de l'hydrogéologie et des conditions 
limites, 

Après avoir examiné les différentes solutions possibles, et après avoir enu compte de l'es- 
timation des coûts de réalisation et du calendrier des travaux, le choix du meilleur traite- 
ment par injections est arrêté. 

Une fois la zone à consolider identifiée, il convient d'effectuer des études spécifiques por- 
tant sur la définition des caractéristiques mécaniques et hydrauliques du sol er d'établir les 
objectifs du traitement. Une analyse approfondie des résultats issus des reconnaissances et 
essais in situ, des essais en laboratoire et des calculs réalisés à l'aide de modèles techniques, 
il est possible de mettre en place le projet de rraîtemen du sol à l'aide de résines expansives. 


L. Remplissage des cavités 


La technologie mise au point par Uretek, dénommée Cavity Filing”, a pour objectif d'ob- 
tenir Le remplissage complet et la stabilisation de cavités souterraines et de tout autre gale- 
rie et vide rechniques (vide sanitaire par exemple) à l'aide du soufflage de billes d'argile 
expansée, à granulométrie contrôlée, complété par une injection de résine fortement 
expansive (polyuréthane). L'action combinée des deux matériaux employés, donc les 
caractéristiques techniques sont présentées en figure 2.2, apporte une solution complète 
au problème. 


Agile expansée Résine Uretek Geoplus® 
Mas volumique du maté € 20 ENI | Psion de gonflement 10000 kPs 

k __ (500 kg/m | maximale (100 kg/em?) 

Angle de fouementimeme | = 40 Temps d'expansion Burèmement court 
Absorption d'u (EN 13059) | <75% | Rpandenvohuméiquedel 2; 20 

Module de déformation Md | >25MPa | Module éstque comparable | 102 180 Ma 

avec Da > 80% (250 kgem | à cclui d'un sol de fondation _ | (100 à 1 800 ge?) 
Résisance À à ruprure > 1300kPa Poids volumique delarésine | 113 Nm 

EN 13055) G 13kgfem) | epansée {20 à 300 kg/em) 


Figure 2.2 Caractéristiques techniques de l'argile expansée et de a résine Uretek Geoplus utilisées 
dans la technologie Cavity Filing. 


La force d'expansion élevée de la résine a pour objectif de remplir le volume des vides inter- 
granulaires, de compacrer les grains constituant le squelerre et de permeure une mise en 
compression des parois de La cavité, conformément aux objectifs préalablement déterminés. 
Le remplissage de la caviré avec de l'argile évire l'apparition de foncis (eHondrements du 
sol en surface) liés à des effondrements soudains du toit de la cavité, tandis que la mise en 
compression produite par la résine expansive empêche l'apparition d'un affaissement en 
surface dû à des déformations des parois de la cavité 


Réduction de la conductivité hydraulique 


L'intervention consiste en une première phase de vérification des dimensions et de la géo- 
métrie de la cavité souterraine; à la suite de quoi, un maillage de percement est réalisé 
pour le soufflage de l'argile expansée et l'injection de la résine. Le souflage de l'argile 
expansée au moyen d'une buse directionnelle permet d'obtenir un remplissage optimal de 
la cavité, normalement jusqu'à 98 % de son volume initial. Par la suite, l'injection de la 
résine expansive entraîne le remplissage total de la cavité générant une mise en compres- 
sion des parois de ladite cavité par clavage, et ce conformément aux objectifs projetés. Les 
principales phases du remplissage des cavités souterraines sont schématisées sur la 
figure 2.3. 


Pros 1 
Remplssage avec de billes 
d'argile expansés 


Phaso 2 
Inioton de résine xpansive 
Urstek Gaoplus” 


Figure 2.3 Principales phases de la technologie Uretek Cavity Filing. 


L'intervention présente de nombreux avantages, notamment la restauration de la conti- 
nuité entre la cavité comblée er Le sol, ainsi qu'une distribution uniforme des contraintes 
à l'intérieur du sol proprement dit. 

De plus, l'expansion de La résine détermine la compression du sol environnant, ce qui pro- 
voque une amélioration des caractéristiques de rigidité et de résistance du sol. Étant 
donné que l'argile expansée, tout comme la résine polyuréthane, est caractérisée par un 
poids spécifique inférieur à celui de l'eau, la particularité de cette méthode consiste à ne 
pas surcharger le sol en dessous de la cavité proprement dite, 


[3 Réduction de la conductivité hydraulique 


Les résines Uretek ont une suucture presque entièrement composée de cellules fermées et 
imperméables à l'eau, une caractéristique qui tend à faire de ces résines un produit presque 
éanche et imperméable à l'eau. Les injections de résines expansives de polyuréthane 
peuvent être utilisées pour remplacer ou éliminer l'eau du sol et, par conséquent, réduire 
la conductivité hydraulique du milieu traité. 


Critères de conception 


D'après les résultats expérimentaux obrenus sur des argilesfissurées par Buzzi er al. (2008), 
la résine Urerek Geoplus® est en mesure de pénétrer dans des fissures d’une largeur égale 
41/10 mm, et l'interface argile-résine atteint une épaisseur comprise entre 1 et 3 mm. La 
pénétration de la résine dans les macropores de l'argile entraîne une diminution de la per- 
méabilité, allant essentiellement de pair avec l'état de la fissuration (nombre, longueur et 
ouverture des fissures) 

D'après l'analyse des résultats des essais de perméabilité (figure 2.4) réalisés sur des échan- 
tillons de terrain sous-jacent aux injections de résine polyuréthane, Buzzi er al (2010) ont 
observé que l'injection de résine dans un sol est susceptible de réduire la perméabilité 
selon un facteur d'environ 50. D'autres tests en cours de réalisation (cf. Uretek France) 
tendent à confirmer ces résultats. La perméabilité intrinsèque K (m?) est corrélée à la 


conductivité hydraulique # (m/s) de la relation 4 = K =, où ju, (Pas) est la viscosité 
LA 


dynamique de l'eau et y, (Pa/m) le poids spécifique de l'eau. 
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Figure 2.4 Évolution de la perméabilité intrinsèque estimée par des essais de perméabilité à l'air Les cercles 
vides représentent des échantillons de sol non traité, tandis que les cercles noirs (pleins) représentent 

des sols traités; le carré vide représente la perméabilité intrinsèque d'un sol non fssuré, Profondeur de 
lissuration: 700 mm; profondeur d'injection: 250 mm (modifié par Buzzi et al, 2010) 


Les injections de résine expansive constituent une solution particulièrement intéressante, 
permettant de réduire la conductivité hydraulique et de combler Les fissures dans Les sols 
sensibles aux phénomènes de rerrait-gonflement. 


Dans l'étude menée par Pasquetto er al. (2008), il a été démontré que l'augmentation de 
la densité des sols fissurés, induire par l'expansion de La résine dans les fissures, protégeait 
contre le risque de variations importantes et successives de volume. Par ailleurs, la substi- 
tution de la résine à l'eau interstiielleréduie les ffets produits par la variation de la teneur 
initiale en eau. 

La densification du sol et la réduction de la teneur naturelle en eau induites par un traite- 
ment d'injections réalisées à l'aide de résine Urerek Gcoplus” entraînent une diminution 
considérable du potentiel de retrait du sol en cas de variations des conditions climatiques, 
telles que les périodes de sécheresse prolongée. 


Alégement 


Dans le cas des interventions de consolidation des ouvrages de fondation, les injections 
réalisées à l'interface sol-fondation permettent de remplir les vides présents directement 
sous es fondations et d'imperméabiliser l'interface proprement dite. 

Les résines polyuréthanes expansives peuvent, en outre, être utilisées dans d’autres appli- 
cations (par exemple, dans le cas d'ouvrages de soutènement) pour remplir les vides à l'in- 
erface entre la structure et le sol, en réduisant la possibilité de formation de parcours 
préférentiels d'infilerations le long de ces discontinuités 


[4 Allègement 


Les injections de résine Uretek peuvent se révéler particulièrement cfficaces pour l'allège- 
ment des couches de sols très compressibles qui se trouvent sous des couches de sols gros- 
siers sensibilisées par la présence de la nappe. 


La figure 2.5 représente un exemple de cas dans lequel une couche de sables limoneux est 
situé au-dessus d'une couche d'argile normalement consolidée, très compressible et sou- 
mise à des déformations visqueuses, er dans laquelle Le niveau de la nappe a une incidence 
sur La couche de sables plus superficielle, En raison des charges transmises par la couche 
de sables limoneux, les argiles normalement consolidées qui se crouvent en dessous sont 
soumises à des affaissements significatifs. 
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Figure 2, Stratigraphie de référence. 


L'injection à l'aide de résine expansive permettra d'obtenir une réduction des charges 
transmises aux argiles et, en conséquence, d’en limiter l'affaissement progressif. 


En efft, de par leur poids spécifique significativement inférieur à celui de l'eau, les injec- 
tions de résine polyuréthane permettent d'alléger la couche sableuse et de réduire les charges 
ransmises aux argiles de fondation. Les injections se propagent dans le sol par imprégna- 
ion, expansion de cavité et fracturation subvertcale pénétrant dans les pores du sol entou- 
rant Les fractures, en refoulant l'eau qui était initialement contenue dans lesdits pores. 


Grières de conception 


En conséquence, les injections d'allègement permerrent d'obrenir: 

+_le remplacement de l'eau contenue dans les pores de la couche sablonneuse, qui se carac- 
cérise par un poids par unité de volume y, = 10 KN/m tandis que la résine polyuré- 
thane peut atteindre des valeurs de poids par unité de volume y,;= 1,2 — 1,5 KN/m°; 

+ le déplacement du sol (solides er eau), la résine remplaçant ces derniers en entraïnant 
un allègement significatif du poids total appliqué sur les argiles sous-jacentes, 


[5._ Amélioration des propriétés mécaniques 


L'expansion de la résine à la suite de son injection dans le sous-sol peut être assimilée à 
l'expansion d'un volume de forme irrégulière, occupant selon une première approxima- 
don une cavité sphérique ou cylindrique, en fonction de l'emplacement et du mode 
d'exécution des injections. L'effet sur le sol environnant se caractérise par une augmenta- 
tion des deux paramètres d'érat de référence de tout milieu particulaire (c'est-à-dire, la 
densité er l'état de contraintes de confinement) avec une amélioration consécutive des 
caractéristiques de rigidité er de résistance d’un volume significarif de sous-sol situé sous 
les fondations de la structure endommagée, Ce sype d'approche conceptuelle peut être 
totalement assimilable à celles adoptées pour les interventions faisant appel à des injec- 
tions de compactage, ou injections solides (compaction grouting), consistant à injecter sous 
pression un mortier très visqueux dans le sous-sol (Fori 8: Manassero, 2009). 

Pour les interventions destinées à améliorer les caractéristiques mécaniques des sols de 
fondation, la technologie employée est la technologie Uretek Deep Injections* (figure 2.6), 
une méthode de consolidation non invasive, adaptée tant aux sols granulaires que cohé- 
sifs, qui s'applique également à tout type de fondations (semelles isolées, semelles filantes, 
dalles, remblais, voiries.…), quel que soit son matériau de construction. 

L'intervention prévoit tout d'abord de percer des trous d’un diamètre inférieur à 30 mm 
et avec un entraxe de 50 à 150 cm, à travers des fondations, de manière à atteindre avec 
précision le volume de sol à traiter; des tubes sont ensuite insérés à l'intérieur de ces 
orages afin d'injecter la résine Uretek Geoplus* dans le sol. 


Les injections sont effectuées alors que la résine est encore en phase liquide mais déjà en 
phase d'expansion. Le volume de la résine augmente de cinq à quinze fois en quelques 
secondes, er la force de gonflement développée peut atreindre une valeur de 10000 kPa, 
qui s'autorégule en fonction de la résistance rencontrée. 

L'expansion de la résine se poursuit jusqu'à ce que le terrain traité soit densifié au point de 
refuser une compression supplémentaire et de contraindre La résine à poursuivre son expan- 
sion dans Le sens de moindre résistance, ce qui provoque un soulèvement de l'édifice. Des 
récepteurs laser dotés d'une sensibilité au demi-millimètre près, sont installés sur l'ouvrage et 
permettent d'exercer une surveillance continue de l'injection: dès que l'objectif de relevage 
est atteint, l'injection est interrompue er Le traitement est arrété en toute sécurité 


Il convient de souligner le Fait que les injections sont concentrées dans un bulbe de pres- 
sion significatif, que lon peut évaluer dans un tout premier temps selon la théorie de 
Boussinesq (comme représenté sur la figure 2.7), c'est-à-dire dans le volume de terrain 
qui est sollicité majoritairement par les contraintes induites par la charge sus-jacente, En 
outre, la rapidité de la réaction d'expansion et de solidification permet de parfaitement 
contrôler la diffusion de la résine lors de l'injection. 


Amélioration des propriétés mécaniques 


Figure 2. Schéma d'intervention avec injection destinée à améliorer les caractéristiques mécaniques 
du sol sous les fondations. 


ÉTAT DESTENSIONS SOUS LA FONDATION 
ET SCHÉMA D'INTERVENTION 


Figure 2,7 État de contraintes sous la fondation et schéma d'intervention, 


Critères de conception 


L'efficacité de l'intervention est vérifiable non seulement par le début de soulèvement 
observé lors de chaque injection, mais également par des essais géorechniques in si, 
comme, par exemple, des essais comparatifs au pénétromètre ou au pressiomètre. 

Le procédé Uretek Deep Injections” peut constituer une alternative particulièrement inté- 
ressante aux technologies traditionnelles, telles que les plots béton, les micropieux ou le 
procédé jer-grouting.…, pour les interventions de consolidation des sols d'assise des bâti- 
ments endommagés par des affaissements différentiel lés à de mauvaises caractéristiques 
mécaniques du terrain d'injection. Les injections de résine se révèlent souvent avanta- 
geuses aussi bien en termes de coût qu'en termes d'invasivité limitée et de durée des tra- 
vaux réduite. 


Un exemple intéressant d'intervention effectuée avec la technologie Uretek Deep Injec- 
tions” est la consolidation du sol des fondations réalisée sur le bâtiment historique (repré- 
senté sur figure 2.8) de la Punta della Dogana à Venise (Gabassi er al, 2011). 


Figure 2. Punta della Dogana, Venise: édifice affecté par un aflaissement de ses fondations 
(Gabassi et al, 2011). 


En mai 2003, le bâtiment a été rouché par un affaissement majeur et soudain de ses fon- 
dations au niveau de la partie restreinte de l'édifice donnant sur le Grand Canal, en rai 
son d'une série de travaux d'entretien du lit du canal, et de l'apparition d'un phénomène 
de siphon lors de l'assèchement de la nappe phréatique des fondations du mur du canal. 
Des études approfondies ont été effectuées, validées par un nombre important d'essais 
gévtechniques in sit et en laboratoire, mais aussi par une surveillance de précision. Après 
avoir évalué les différentes technologies possibles et idoines pour stabiliser Le bâtiment, 
à été décidé de recourir à des injections de résine expansive au moyen de la technologie 
Uretek Deep Injections*. 


Amélioration des propriétés mécaniques 


Il a été jugé approprié de mettre en œuvre des interventions destinées uniquement à annu- 
ler la tendance aux déformations en constante évolution (figure 2.9), et, dans la mesure du 
possible, de compenser partiellement ces effes sur les structures impliquées avec un soulè- 
vemenr limité des parties fortement afaissées. 
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Figure 2 Évolution au fl du temps de l'affissement de a façade sur le Grand Canal de Venise 
(Gsbassi et al, 2011). 


Il a donc été choisi d'intervenir dans le volume situé immédiatement sous les fondations. 
Les injections ont été effectuées dans des volumes de sol affectés par le phénomène de 
siphon, à des intervalles réguliers sur un certain nombre de niveaux à une profondeur 
comprise entre 4,0 m minimum et 7,0 m maximum (figure 2.10). 


Il convient, par ailleurs, de souligner Le fait que les injections ont été réalisées en plusieurs 
étapes, afin de permertre La dissipation complète des surpressions interstitielles induites 
par les injections (Gabassi es 44, 2011) 

L'intervention a été divisée en trois phases différentes, réalisées dans un laps de temps de 
six mois, et a permis de restaurer la stabilité et la fonctionnalité de l'édifice, mais aussi 
d'arrêter l'évolution des affaissements, en rétablissant les conditions originelles du sol. 
Pendant toute la durée des travaux mais aussi entre les diverses phases, le comportement 
de la structure a été mis sous surveillance, en plus du système de surveillance déjà en place, 
un relevé altimétrique de cibles solidement fixées aux murs, à l'aide de niveaux optiques 
de précision. 
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Section tramsvurese Section longitudinale 


Figure 2.10 Coupe type représentant la répartition des injections — phase 1 de l'intervention, à sept niveaux 
(Gabassi et al, 2011) 


Les effets positifs induits par les injections de résine dans le sol ont été mis en évidence, 
autrement qu'à travers les résultats du système de surveillance (figure 2.11), par la réali- 
sation d'essais de pénétration dynamique comparatifs pendant er entre les phases. En plus 
de ces essais, des sondages ont été effectués afin d'identifier le matériau injecté et d'en 
effectuer des prélèvements pour en déterminer les propriétés mécaniques. 


sors 


Figure 2.11 Évolution de l'titude d'un point du bâtiment suvelé au fl du temps, correspondant au paint 
ayant subi des aflaissements majeurs - façade du côté du Grand Canal (Gabessi et a, 2011). 


Relevage 


ls. Relevage 


La technologie Uretck Floor Lift" (figure 2.12) est de plus en plus utilisée pour relever et 
mettre à niveau des dallages (rigides er flexibles) des bâtiments industriels, de magasins et 
d'entrepôts, d'immeubles résidentiels, de routes er de pistes d'aéroport concernés par des 
phénomènes de subsidence, 
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Figure 2.12 Technologie Uretek Floor LI (mmwureteki). 


La technologie Floor Lif” prévoit une première phase de relevage dans laquelle des instru- 
ments laser seront utilisés dans les zones affectées par l'afaissement, afin de mesurer avec 
une précision millimétrique les affaissements survenus. De plus, au niveau des joints, on 
se sert d'un comparateur d'une précision de l'ordre du centième de millimètre. 


Ensuite, le dallage est percé avant la mise en place de tubes métalliques de petit diamèrre, 
disposés selon un maillage régulier. 


La résine est ensuite injectée jusqu’à provoquer le soulèvement grâce à la capaci 
sion de la résine logée entre la dalle et la plate-forme d'assise (figure 2.13). 


expan- 


La résine st njociéo à raw des tubes de 6 à 12 mm au niveau de l'merfsce so dalle 
tance un mon tring ler canstams urant l'invention. 


Figure 2.13 Intervention à l'aide de la technologie Uretek Floor LiR® (uwuuretckit). 


Grtères de conception 


Le relevage est contrôlé à l'aide d'un système de surveillance laser qui permet d'avoir une 
précision de l'ordre du demi-millimètre, Le résultat en est un dallage avec une planéité 
restaurée et un plan d'assise homogène ayant de meilleures caractéristiques de stabilité. Si 
l'on applique de surcroît le procédé Uretek Deep Injections', il est possible, en outre, 
d'augmenter la portance de l'asise en profondeur. 

Un exemple significatif de soulèvement du dallage d'un bâtiment industriel est présenté 
ci-après (figure 2.14) 


Figure 2.14 Planimätie du bâtiment industriel soumis à l'intervention avec la technologie Uretek Floor Lift®. 


L'intervention à été effectuée avec Les technologies Uretek Floor Lifi* et Uretek Deep 
Injections" et a concerné l'injection, dans la couche d'appui du dallage, d'une résine poly- 
uréthane à pression de gonflement élevée, qui a permis de rétablir Le contact entre la plate- 
orme et Le dallage, mais aussi d'obtenir un relevage de ce dernier. 

Avant l'intervention, différents essais au pénétromèure ont été effectués dans le voisinage 
immédiat, à l'intérieur d’un domaine d'essai d’une surface de 60 m°. 

Les injections ont été effectuées en sous-face du dallage et dans le fond argileux, selon un 
maillage avec un entraxe d'environ 1,0 m, après avoir percé la dalle en béton armé consti- 
tuant Le dallage proprement dit. Des injections ultérieures ont été effectuées Là où il était 


Relevage 


nécessaire de récupérer un affaissement majeur er dans les points où, au cours des pre- 
mières injections, le relevage de la structure ne s'était pas révélé homogène en raison de la 
présence de discontinuités structurelles ou d'une hétérogénéité du sol. 

Des percements d'un diamètre de 18 mm ont été effectués à travers la semelle formant le 
dallage, et leur longueur au-delà de l'épaisseur de la dalle n'excédait pas 1,1 m. Le plan 
d'appui du dallage s'est révélé se trouver À environ 25 em de profondeur par rapport au 
niveau de la surface du sol. 

Les niveaux d'injection se situaient à une hauteur de 0,1 m et 1,1 m par rapport au plan 
de pose du dallage (figure 2.15). Par ailleurs, les injections ont été tout d'abord effectuées 
dans les niveaux les plus superficiels, puis dans les niveaux plus profonds. 
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Figure 2.15 Sdhématisation des injections. 


La consommation totale enregistrée pendant la phase d'injection a été de l'ordre de 
14800 kg de matériau expansé, pour un total d'environ 2 600 m° de dallage relevé. 

Au cours des travaux, Le dallage a fait l'objet d'une surveillance au moyen d'instruments 
laser de précision, lesquels ont permis de relever constamment les déplacements verticaux. 
L'interruption de l'injection de résine a été généralement effectuée en fonction de la dérec- 
tion du relevage préalablement fixé pour Le dallage et, dans certains cas, en fonction du 
seuil de pression de la pompe d'injection de la résine. 

La figure 2.16 (page suivante) représente Les courbes de niveau du dallage avant er après 
l'injection, ce qui permet de valider l'intervention. 


Grtères de conception 


| 
4 


& 


A o 


Figure 2.16En 
après intervention. 


courbes de niveau du dallage avant intervention; en bas: courbes de niveau du dallage 


RCE . 
Autres typologies 


d'injections 
et comparaisons 


[1 _ Introduction 


Le présent chapitre prend en compte les possibilités d'utilisation des rechniques d'injec- 

ion de résines polyuréthanes expansives en faisant, en outre, quelques comparaisons avec 

Les techniques d'injection traditionnelles les plus courantes actuellement. 

En général, les traitements des terrains par des techniques d'injection peuvent être classés 

de manière synthétique en 

+ injections par imprégnation (ou permeation grouting) pour la réduction de la conduc- 
rivité hydraulique er le renforcement des terres; 


+_ injections de compactage (ou compaction grouting) pour la consolidation des sols; 

+_ injections de compensation (ou compensation grouting) pour le contrôle et la récupéra- 
tion des affaissements de bâtiments er d'autres structures, 

Les finalités er les objectifs des trairements, en plus des caractéristiques hydrogéologiques 

et géotechniques du sous-sol, constituent les éléments fondamentaux pour faire Les choix 

Les plus opportuns et les plus judicieux au regard du type de système et des composés d'in- 

jection à utiliser. 

Les techniques traditionnelles d'injection font généralement appel à un mélange d'eau ct 

de ciment {avec des additifs Auidifiants), à des composés de silicare ou à des résines (non 

expansives) 


n présence des mêmes caractéristiques rhéologiques des mélanges, les techniques rradi- 
tionnelles d'injection permettent le contrôle opérationnel du traitement uniquement à 
uravers Le réglage de la pression et des débits à la sortie du système de pompage. 

Par contre, comme cela a déjà été montré dans Le chapitre 1, les résines polyuréthanes 
se caractérisent par de remarquables propriétés de gonflement qui se développent durant 
le processus de solidification. Par conséquent, à la suite de la phase d'injection et de pro- 
pagation dans le sol en conditions fluides, l'expansion de la résine continue jusqu'à 
atteindre l'équilibre des contraintes avec Le sol environnant au moment où la pression de 
gonflement de la résine, dérivant de la réaction de synthèse, coïncide avec la pression 
moyenne de confinement du matériau environnant. 

Cette caractéristique particulière des résines polyuréthanes offre des possibilités opéra- 
tionnelles intéressantes de réglage et de contrôle du traitement, mais nécessite, pour sa 
mise au point et son calibrage, des connaissances et des expériences approfondies en ce 
qui concerne les caractéristiques rhéologiques du mélange d'injection et celles géotech- 
niques du terrain environnant. 
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Dans les comparaisons avec les rraîrements traditionnels, la technique du jer-grouting n'est 
pas prise en compte, parce que, prévoyant le mélange avec Le sol dans lequel est effectuée 
l'injection, elle se caractérise par des différences significatives vis-à-vis des applications 
basées sur l'utilisation des résines expansives. 


L. injections par imprégnation 


Les injections par imprégnation sont employées, en règle générale, pour réduire la conduc- 
divité hydraulique naturelle du sol, ou pour créer des éléments de renfort à l'intérieur 
même du sol qui interagit avec Les structures à stabiliser. 

Dans la technique traditionnelle du permeation grouting, le choix des composants du 
mélange est très important, car ce dernier doit pénétrer à l'intérieur des pores sans altérer 
la structure et l'état de contraintes du sol. Ce choix dépend aussi bien de la conductivité 
hydraulique du sol que de la technique urilisée pour l'injection, et a une grande inuence 
sur Le résultat final du traitement. Si la viscosité du mortier est trop élevée ou que la gra- 
nulométrie des éléments en suspension est trop grossière, la pénétration complète et uni- 
forme est impossible, ec il ne sera donc pas possible d'arreindre les objectifs du projet. 

La capaciré de gonflement des résines polyuréchanes expansives empêche leur utilisation 
pour réaliser de simples injections d'imprégnation, puisque, inévitablement, dans la phase 
de solidification et de gonflement, elles modifient l'état de contraintes et de déformations 
du sol. Quelques analogies avec les principes de fonctionnement des injections par impré- 
gnation peuvent néanmoins être identifiées dans le cadre de la technologie Caviry Filling® 
développée par Uretck pour le remplissage de cavités souterraines. 

En effer, cetre rechnologie n'a pas non plus pour but d'améliorer les caractéristiques méca- 
niques er de modifier l'éar de contraintes et de déformations du sol, sauf pour ce qui est 
de Le ramener dans les conditions antérieures à la formation de La cavité. L'expansion de la 
résine polyuréthane, dans ce sype particulier d'applications, est surtout exploitée pour 
obrenir le remplissage complet de la cavité et non pas pour provoquer le compactage du 
sol environnant la cavité elle-même. 


Enfin, il peut aussi être observé que le remplissage de cavités souterraines ou de fissures 
qui se sont formées à l'interface sol-structure, produit une réduction de la conductivité 
hydraulique à grande échelle, en empêchant le passage de l'eau à travers de possibles voies 
préférentielles de filtration. 


[3._ injections solides (injections de compactage) 


Les injections solides constituent un outil efficace pour améliorer les caractéristiques 
mécaniques du sol, grâce à la densification des grains provoquée par la pression transmise 
par le matériau injecté dans le sol. 

L'injection solide, ou sompaction grouting, wraditionnelle est la technique la plus utilisée 
pour densifier (compacter) Le sol, dans l'optique d'en augmenter la rigidité et la résistance 
au cisaillement, Cette technique consiste à injecter un mortier très visqueux et à même de 
former un bulbe, dont l'expansion provoquera la déformation et la densification du sol 
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environnant. L'expansion du bulbe est contrôlée par la pression d'injection du mortier et 
par les caractéristiques mécaniques du sol dans lequel s'effectue l'injection. 

La méthodologie la plus courante pour la validation de la technique d'injection solide 
(compaction grouting) est la comparaison des essais réalisés in situ avant et après le traîte- 
: par exemple, des essais de pénétration statique (CPT, Come Penetration Ti) ou de 
pénétration (SPT, Standard Penetration Tes). Le concepteur définit en principe pour le 
paramètre de comparaison une valeur minimale, considérée comme le niveau minimum 
d'amélioration requise. 


Même les interventions avec des résines polyuréthanes expansives, en particulier avec les 
technologies Urerek Deep Injections’, se basent principalement sur des principes de fonc- 
tionnement analogues. 

Un effet significacif de compacrage peut être obrenu au moyen des traitements Uretek 
Deep Injecions* dans les sols à grains grossiers (sable er gravier), dans lesquels, après une 
pénétration initiale de la résine, se forme un véritable bulbe à l'intérieur même du sol. 
L'expansion de la résine dans la phase de solidification provoque l'augmentation du 
volume du bulbe, ce qui se traduit par un compactage du sol environnant. 

Dans le cadre de a modélisation du projet, on peut généralement supposer que l'expan- 
sion de La résine est, dans un premier temps, comparable à l'expansion d'une cavité cylin- 
drique ou sphérique, et que le processus s'arrête quand un état d'équilibre entre la pression 
de gonflement de la résine et la tension moyenne de confinement mobilisée dans le sol 
alentour est atteint. L'évolution du rapport entre le volume de résine final et Le volume de 
résine initial peur être évaluée en fonction de la pression de gonflement de la résine er de 
la contrainte moyenne de confinement du sol, à la valeur P* conformément à la représen- 
tation de La figure 3.1 
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Figure 3.1 Évoltions qualitatives du rapport de gonflement (V,/V,) en fonction de a pression 
de gonflement de la résine et de l contrainte moyenne de confinement du sol. 
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En fonction de la pression d'équilibre, P”, atteinte à la fin du processus d'expansion, il est 
possible d'évaluer les variations des contraintes de confinement du sol à partir desquelles 
il est possible d'estimer l'augmentation de la résistance au cisaillement et de la rigidité du 
sol autour de la cavité. 


Il convient de souligner également que, contrairement aux injections avec d'autres types 
de résines, caractérisées par une plus faible capacité de gonflement, la technologie Uretek 
permet d’anindre des pressions d'équilibre, P*, ayant des valeurs considérablement éle- 
vées, comme cela est montré de manière détaillée dans les chapitres suivants 

À titre d'illustration, on peut se reporter à la figure 3.2, dans laquelle la courbe rhéolo- 
gique de la résine Uretek est comparée à celle d'une résine dotée d'une moindre capacité 
de gonflement: l'intersection avec la courbe représentant la contrainte moyenne de conf 
nement du sol apparaît, pour la résine Uretek, à une valeur de pression et un rapport de 
gonflement nettement supérieurs à ceux associés à la courbe d'une résine à capaci 
gonflement plus limitée. 
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Figure 32 Évoltions qualitatives du rapport de gonflement V,/V,en fonction de La pression de gonflement 
de la résine Uretek (courbe descendante continue) et des autres résines (courbe pointilée descendante), 
t de la contrainte moyenne de confinement du sol (courbe continue ascendante) 


Puisque l'amélioration des caractéristiques mécaniques du sol environnant le bulbe impré- 
gné par le mélange d'injection provient de la densification induite par l'expansion du 
bulbe lui-même et de l'augmentation de la contrainte de confinement, il est donc évident 
que la résine Uretek peut fournir des résultats considérablement supérieurs à ceux d'une 
résine ayant des propriétés de gonflement significativement inférieures. 
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4. Injections pour le contrôle et 
la compensation des affaissements 


Les injections de compensation ont été utilisées pour compenser les affaissements de 
structures provoqués essentiellement par l'excavation de galeries souterraines. 

La figure 3,3 montre un exemple de structures construites sur des pieux de fondation 
affectés par des affaissements dus à des excavations souterraines. Dans ce cas et dans des 
as similaires, il est possible de recourir à des interventions par injections de compensa- 
don (compensation grouting) pour limiter l'affaissement des structures (Mair er al, 1994; 
Harris et al, 1994; Harris et al, 1996; Mair 8e Taylor, 1997; Harris, 2001 ; Mair, 2008; 
Mair et al, 2009; Pigorini et al, 2009; Bezuñjen er al, 2011). 
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Figure 3.3 Exemple d'injection de compensation (compensation grouting) pour compenser es afaissements 
induits par l'excavation d'un tunnel du mêtr d'Amsterdam (modhfi£ de Bezuijen et al, 2011). 


La figure 3.4 (page suivante) montre les principaux aspects er les composants qui carac- 
térisent la technique des injections de compensation. 
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Figure 5.4 Principaux aspects et composants de la technique d'injection de compensation 
{modifié de Mar et al, 2009). 


Les principales étapes de traitement requises pour ce type d'intervention peuvent être 

résumées comme suit: 

— Des essais sont effectués sur site pour définir avec une extrême précision le comporte- 
ment des mélanges à injecter et, en particulier, pour vérifier Les mécanismes de propa- 
gation er de répartition dans le sol, l'évolution de la viscosité de la phase liquide, les 
temps de solidification/durcissement et, dans Le cas de résines expansées, le comporte- 
ment rhéologique de la phase d'expansion. 

— Préalablement à lexcavation du tunnel et après l'installation des systèmes de surveil- 
lance en profondeur comme en surface, des forages sont exécutés, et des tubes d'injec- 
tion avec vannes ou tubes à manchette sont positionnés entre a future excavation et les 
fondations de l'ouvrage surmontant Le futur tunnel. 

— Ensuite, est injectée, à basse pression, une quantité contrôlée de mortier en vue de den- 
sifier Le sol, comblant ainsi es vides dans le sol et provoquant même une augmentation 
minime de l'état de contrainte. 

Lors de l'excavation de galeries souterraines, on injecre une certaine quantité de mortier 
au moyen de tubes spéciaux (appelés tubes à mancherre, ou TAM), installé entre la gale- 
rie et les fondations des structures, pour compenser Les affaissements et la variation de 
l'état de contraintes induire par l'excavation. L'installation des TAM peut provoquer des 
affaissements, que l'on peut toutefois compenser de diverses façons. Avant de commencer 
l'excavation souterraine, on injecte un mortier de conditionnement pour resserrer Le sol, 
remplir d'éventuels vides et contrebalancer les affaissements induits par l'installation des 
TAM. Ensuite, on procède simultanément à l'injection du mélange et à l'excavation dans 
l'optique de limiter à la fois les afaissements et les déformations des structures, Avant de 
procéder à l'installation des TAM, il est indispensable d'effectuer une série d'essais sur site 
afin de rester les mélanges à injecter et, en particulier, d'en vérifier la propagation effec- 
tive. La formation de fractures horizontales, dans lesquelles se propage Le mélange, pro- 
voque le soulèvement d’une zone étendue et permer de compenser les afaissements 
induits par les excavations souterraines. 
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Mème pour ce type de traitement, les paramètres permettant de contrôler la propagation 
du mélange injecté ne sont autres que les caractéristiques physiques du fluide, notamment 
la viscosité et la pression d'injection. 

Les résines polyuréchanes représentent une bonne alternative à la technique traditionnelle 
des injections de compensation, en particulier dans les sols à grains fins (sables limoneux, 
limons et argiles) saturés, au sein desquels la résine ne forme pas de bulbe, mais où elle se 
propage à l'intérieur de fractures qui se produisent au cours de son expansion. En effet, si 
le milieu est saturé, en fonction de la conductivité hydraulique et de la cension capillaire, 
il n'est pas dit que l'on réussira à refouler l'eau présente dans les pores, et on peut procéder 
à la fracturation du sol. Une telle fracturation est d'autant plus probable que la portion de 
so fin, caractérisé par une faible conductivité hydraulique, est plus significative. 

Après la réalisation d'un essai sur le terrain par Uretek Nederland, en collaboration avec 
l'Université de technologie de Delfi (TU Delf), on a pu observer comment, en fonction 
de l'état de contraintes, de la conductivité hydraulique et de la tension capillaire du sol, 
dans les sols cohésifs la résine Uretck Geoplus ne parvenait pas à déplacer l'eau présente 
dans les pores et à former un bulbe comme dans le cas décrit au paragraphe «Injections 
solides (injections de compactage) ». Au contraire, la propagation du mélange a lieu seu- 
lement à l'intérieur des fractures (figure 3.5) préexistantes ou qui se sont formées au cours 
du processus, et la compensation des affaissements des structures de surface est obtenue 
lors de la formation de fractures horizontales dans lesquelles l'expansion de la résine peut 
se produire. 


Figure 3,5 Excavation d'inspection à la suite de l'essai effectué sur le terrain par Uretek Nederland, 
en collaboration avec l'Université de technologie de Del. Les fractures verticales et horizontales 
sont évidentes, pour a plupart de forme discoidale. 


Le schéma de la figure 3.6 (page suivante) présente les principales phases du processus 
d'injection, d'expansion et de propagation de La résine. L'injection est effectuée dans un 
sol normalement consolidé (degré de surconsolidation, OCR = 1, coefficient de pression 
des terres au repos, &, < 1), caractérisé initialement par des contraintes géosratiques hori- 
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zontles, Gjy moindres que les verticales, 0,9. Au début du traitement, il se forme une 
fracture verticale, puis se produit l'expansion de la résine à l'intérieur de cette fracture, ce 
qui en augmente l'ouverture. Cette expansion provoque une augmentation des contraintes 
horizontales avec, en conséquence, une inversion des contraintes principales (0j, = O. 
O3 = 6) dans Le sol environnant. 
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Figure 3.6 Schématisation du processus de fracuration dans un sol cohésif normalement consolidé. 


En raison de la rotation des contraintes principales, l'injection ultérieure de résine à pro» 
mité du premier point permet de provoquer la formation d'une seconde fracture dans le 
plan. L'ouverare de la fracture horizontale produit le soulèvement de la fondation et 
compense, en conséquence, les affaissements. Il convient d'insister sur la nécessité de réa- 
liser deux phases d'injection, pour s'assurer avec davantage de certitude de la formation de 


fractures horizontales et donc de l'issue satisfaisante du traitement. 


Modèles théoriques 


[1 Introduction 


Les injections de résines polyuréthanes expansives s'utlisent majoritairement dans la plu- 
part des interventions visant à améliorer les propriétés mécaniques des sols et à remédier 
aux affaissements, complets et/ou différentiels, des fondations superficielles. 

L'approche théorique adoptée dans la conception de ces interventions dépend de la gra- 
nulomérrie du sol dans lequel est réalisé le traitement de consolidation. 

Dans le cas des sols à grains grossiers (sable moyen/grossier et graviers), la conductivité 
hydraulique est sulfisamment élevée pour permettre à la résine de pénétrer dans les pores 
et d'entraîner la formation d'un bulbe constitué d'un mélange sol-résine. L'expansion de 
ce bulbe se poursuit jusqu'à arteindre des conditions d'équilibre avec les contraintes de 
confinement générées dans le sol environnant. 

Dans le cas des sols à grains fins (sable fin, limon ec argile), la résine ne réussi pas à péné- 
rer dans les pores et, par conséquent, son expansion provoque La formation de fractures, 
dont l direction dépend principalement de l'homogénéité et de l'isouropie du sol ainsi que 
de l'état de contraintes initial. La résine se propage dans les fractures, ce qui entraîne des 
variations de l'état de compaction et des déplacements significatifs dans le sol environnant. 
Pour l'analyse des effers produits par les injections de résines polyuréthane, les conditions 
de drainage sont d'une importance capitale, et sont elles aussi fonction de la granulomé- 
tie du sol traité. 

Les sols à grains grossier 
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Dans les sols à grains fins, en revanche, l'expansion de la résine se produit en condition 
non drainée, avec formation de surpressions intersttielles, done la dissipation s'effectue 
dans un laps de temps appréciable. Il s'ensuit, dans ce cas, des variations de l'étar de 
contraintes qui se produisent au fl du temps, et peuvent nécessiter l'exécution d'injec- 
tions différées dans le temps afin d’artcindre les objectifs de consolidation prédéfinis. 
Comme l'évaluation des surpressions incerstitielles générées dans le sol par l'expansion de 
la résine s'avère très difficile, on à généralement recours à une analyse en termes de 
contraintes totales, dans laquelle Le sol est assimilé à un milieu à phase unique, caractérisé 
par un critère de rupture de Tresca (matériau purement cohésif). De certe façon, il est 
possible d'estimer Les effets à court terme produits par les injections, à la fois en termes de 
variation de l'état de contraintes et en termes de déplacements induits dans le sol. 
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Différentes typologies de consolidation découlent des différents modes de comportement 
des résines polyuréthanes dans le sol, L'expansion en condition drainée du bulbe de sol 
traité détermine la densification du terrain environnant, par analogie avec ce qui se passe 
dans Le cas des injections solides. De plus, le bulbe constitue un élément de renfort, carac- 
térisé par des valeurs de rigidité et de résistance au cisaillement supérieures à celles du ter- 
rain environnant. Dans ces conditions, l'exécution d'injections verticales provoque 
surtout des déplacements horizontaux dans le sol, ce qui entraîne une densification, er des 
déplacements verticaux négligeables. S'il s'avère nécessaire de produire des déplacements 
verticaux importants pour compenser Les affaissements de structures de fondations super- 
ficelles, il est indispensable de prévoir la formation d'une discontinuité horizontale, 
moyennant le prétraitement d'une partie du sol au-dessous des ouvrages de fondation : Les 
injections successives, cffectuées au-dessus du sol traité, auront tendance à se propager 
horizontalement et produiront, lors de l'expansion de la résine, le soulèvement souhaité 
Dans les sols à grains fins, la formation de fractures horizontales permet d'expliquer le 
soulèvement du niveau de la surface naturelle (niveau du sol) lors de l'expansion de la 
résine. Les fractures horizontales se produisent lorsque les contraintes horizontales sont 
supérieures aux verticales: lorsque le coefficient de pression des terres au repos est supé- 
rieur à l'unité, cette condition peut être satisfaite à la suite de l'expansion des fractures ver- 
ticales initiales. La densification du sol est d'autant plus considérable que le réseau de 
fractures créé dans le sol est dense, de sorte que des volumes entiers de sol sont compri- 
més entre une fracture et la suivante. 


Dans le cadre du présent chapitre, les méthodes d'analyse chéorique du comportement 
des résines polyuréthanes dans Les sols seront présentées en faisant davantage référence aux 
aspects fondamentaux qu'aux applications. I s'en suivra deux approches théoriques diffé 
rentes, respectivement pour es sols à grains grossiers et pour les sols à grains fins, d'après 
Les résultats obtenus dans une série d'essais d'injection à petite échelle. Pour les sols à 
grains grossiers, il sera fit référence à la théorie d'expansion des cavités, tandis que pour 
les sols à grains fins, il conviendra d'introduire quelques concepis de la mécanique à la 
rupture dans les milieux élastiques. 


[2 Comportement de la résine dans le sol 


Les résines expansives polyuréthanes ont un comportement différent en fonction de la 
granulométrie du sol dans lequel elles sont injectées. 

Afin d'évaluer, ne serait-ce que visuellement, cette différence de comportement, on a 
mené une série d'essais d'injection ponctuelle sur un échantillon de sol à grains fins (argile 
compactée) et sur un autre à grains grossiers (sable sec, déposé par gravité), contenus res- 
pectivement dans deux moules en acier, ayant chacun un diamètre de 27,5 em et une hau- 
teur de 50,0 em. Au cours de la préparation des échantillons dans le moule, le sol a été 
introduit jusqu'à une hauteur inférieure à celle du moule, en laissant un volume libre dans 
la partie supérieure («vide») correspondant à un espace d'une hauteur égale à environ 
3,0 cm. Ce dispositif a été adopté afin d'éviter d'atteindre des pressions trop élevées au 
sein du moule, lorsqu'il est fermé en son sommet, compte tenu de la capacité de gonfle- 
ment de la résine et du degré de confinement de l'échantillon. 
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Dans chaque moule, on a injecté 4 kg de résine Urerek Geoplus” à une profondeur d'en- 
viron 15,0 em à travers le couvercle supérieur du moule, L'injecrion a été effectuée en plu- 
sieurs fois, chacune d'une durée variable de deux à quatre secondes, avec une pause de 
deux secondes entre elles. 


L'injection de la résine dans de l'argile compactée a entraîné l'expulsion d'une certaine 
quantité d'eau dont on a retrouvé là présence dans le fond du moule, démontrant ainsi 
que l'expansion de La résine provoque la formation de surpressions intersttielles dans Le 
sol comprimé. 

À la fin de la phase d'injection, on a ouvert le couvercle supérieur du moule er on à 
constaté la présence d'un bouchon de résine dans l'espace laisé libre pendant la phase de 
préparation de l'échantillon. À l'ouverture du couvercle inférieur, nous avons constaté la 
présence d'une fracture dans le so, à l'intérieur de laquelle la résine s'est propagée, comme 
représenté sur la figure 4.1. 


Figure 4.1 Ouverture du couverde inférieur du moule faisant apparaître l'argile compactée: il est possible 
de voir les traces des deux principales fractures qui se sont propagées jusque dans le fond du moule. 


Par la suite, en reiranc le sol du moule, nous avons constaté l'apparition de deux fractures 
principales (figure 4.2, page suivante), à partir du point d'injection, dans lesquelles la 
résine s'est propagée sur toute la hauteur disponible, pour atteindre en partie les parois 
latérales, au point d'épouser partiellemenc la forme du moule, ainsi qu'on peut le voir sur 
La figure 4.3 (page suivante). 

Ave l'essai d'injection sur des échantillons de sable, le résultat obtenu est sensiblement 
différent. 
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Figure 4.2 Résine solidifiée à l'intérieur des fractures qui se sont propagées dans le sol arileux. 
Vue depuis le ford du moule. 


Figure 4.3 Extraction de la résine hors du moule. 


L'ouverture du couvercle inférieur n'a révélé la formation d'aucune fracture, non plus que la 
présence de résine sur le fond. Dans ce cas aussi, on a retiré le sol, ce qui a permis de décou- 
vrir, à partir d'une certaine profondeur, un agglomérat bulbiforme de résine (figure 4.4). 


Comportement de la résine dans le sol 


Figure 4.4 Agglomérat bulbiforme de résine à l'intérieur du sable. Vue depuis le fond du moule. 


Ensuite, en retirane le sable, on a extrait l'agglomérat bulbiforme (figures 4.5 et 4.6, page 
suivante), à l'extrémité duquel un bouchon constitué uniquement de résine s'érair Formé 
dans l'espace vide supérieur. La longueur du bulbe était de l'ordre de 30,0 em (pour une 
longueur de moule de 50,0 em). 


Figure 45 Taille réelle de l'agglomérat bulbiforme de résine obtenu par injection dans un échantillon 
de sable. Dans ce cas aussi, On à constaté la formation d'un bouchon de résine en contact avec 
le couvercle supérieur, où l'on avait laissé un espace vide d'environ 3 à 4 em. 
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Figure 4.6 Agglomérat bulbiforme de résine extrait du moule, 


Les agglomérats bulbiformes ont ensuite été sectionnés à l'aide d'une scie circulaire, pour 
en observer la composition interne. La figure 4.7 représente une photographie d'un 
bulbe dans lequel deux sections transversales ont été découpées, l'une à proximité du 
point d'injection et l'autre plus en profondeur (ou dans le tronçon final de l'agglomérat). 


Figure 4.7 Photographie du bulbe reconstitué après réalisation de coupes transversales à l'aide d'une scie 
éreulaire, Deux coupes transversales ont êté réalisées: une à proximité du point d'injection et une plus 

en profondeur (Remarque : sur la photographie, le bulbe est présenté à l'envers, cestà-dir la tête en bas, 
servant de socle) 
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Les images des coupes sont présentées sur les figures 4.8 ex 4.9, d'où il ressort clairement 
que la résine a pénétré entre les grains de sable en formant un agglomérat 


Figure 4.8 Coupe transversale de l'agrégat bulbiforme, effectuée à proximité du point d'injection. La résine 
a pénétré entre les grains de sable à pari de la zone d'injection, cette dernière étant reconnaissable 

au fait qu'elle est remplie de résine pure. Aux alentours du point d'injection, la résine, par l'intermédiaire 
du processus de pénétration ou d'imprégnation entre les grains, a formé un agglomérat ave les grains 
de sable. 


Figure 4.9 Coupe transversale de l'agrégat bulbiforme, effectuée au-dessous du point d'injection. La résine 
à complètement pénétré entre les grains de sable 


D'après ces éléments, on peut donc considérer que, dans les sols à texture grossière, il existe 
un processus très rapide de pénétration de la résine entre les grains de sol qui entraîne la 
formation d'un bulbe, lequel, en se solidifiant, s'expanse en compactane le sol environnant, 
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La différence de comportement de la résine, liée à la différence de granulométrie des deux 
sols, nous conduit à envisager deux approches d'étude différentes. 

Dans le cas du sable, la formation du bulbe apparaît suite à la pénétration de La résine en 
phase liquide dans Les «pores» du sol, opération suivie par une phase d'expansion lors de 
la solidification de La résine. La modélisation de l'expansion du bulbe peut être ramenée à 
l'étude du problème de l'expansion d'une cavité (sphérique ou cylindrique) à l'intérieur 
d'un milieu élasto-plastique (annexe À) 

En revanche, pour évaluer les conditions de déclenchement des fractures dans l'argile 
compactée, il est possible de faire appel à la mécanique de la rupture en milieu élasto-plas- 
tique purement cohésif, en se référant à des conditions non drainées (annexe B). Tou- 
tefois, étant donné l'impossibilité théorique de prévoir le nombre, la direction et l'ampleur 
des fractures, la capacité de prévision des fers induits par la réalisation d'injections dans 
les sols à grains fins, pouvant être obtenue avec des modèles théoriques, est forcément 
limitée. 


3. Modèles théoriques pour les sols à grains 
grossiers 


L'étude de la propagation des résines expansives Uretek dans les sols à grains grossiers 
(sable moyen/grossier et gravier) peut être fondée sur la théorie d'expansion des cavités 
dans les milieux élasto-plastiques, en condition drainée, En raison de la conductivité 
hydraulique élevée de ces sols, lors de la phase d'injection, il existe un effet initial signifi- 
catif de pénétration, qui dépend de la viscosité de la résine ainsi que du nombre et de la 
taille des vides du sol. 

On soulignera également que, dans le cas des sols à grains grossiers, toutes les corrélations 
rapportées sont obtenues en référence à des analyses en condition drain. 
quent, les contraintes indiquées doivent être comprises comme étant effectives. 


3.1. Effet de pénétration initial 


Pour les sols à grains grossiers, en raison de la granulomérrie qui leur confère des valeurs 
de perméabilité élevées, le phénomène initial de pénétration de la résine, qui se produit 
suite à l'injection, est important. 

En fonction de la viscosité dynamique de la résine, 1, er du coefficient de conductivité 
hydraulique du sol imprégné par la résine, €, il est possible d'estimer Le rayon du bulbe 
initial 4 (qui est supérieur au rayon du trou d'injection 4), selon que l'on considère une 
imprégnation dans Le cas d’une cavité cylindrique ou sphérique, comme représenté sur la 


figure 4.10. 
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Figure 4.10 Imprégnation dans le cas d'une cavité sphérique (à gauche) et dans le cas d'une cavité 
chindrique (à droite). 


D'après la loi de Darcy et le bilan massique de la phase liquide, on peut établir les équa- 
tions ci-après pour Le rayon du bulbe initial, 4 [m], qui sont utilisées respectivement dans 
le cas d’une injecrion avec une source sphérique [4.1] et avec une source cylindrique [4.2] 


BrakAe , 
FE Aer 


RAENVEN EL ST 


A6)" 6) 


+4, = coefficient de conductivité hydraulique du so! à la résine [m/s]; 
différence de charge engendrée par la pression d'injection p, [ml]; 
porosité du sol [-]; 


Ba Kopor 
m4, 


+4 Hal 
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= y, AH = pression d'injection de la résine [Pa]; 


+ YA = poids volumique de la résine à l'érat liquide [N/m À; 
SD 


+ _K= perméabilité intrinsèque du sol [m’ 
+, = viscosité dynamique de la résine (Vinson & Mitchell, 1972); 
+_4ç= rayon du trou d'injection [m]; 

= temps de réaction de la résine [s]. 


Le coefficient de conductivité hydraulique du sol imprégné par la résine et la perméabilité 
intrinsèque du sol sont liés par la relation suivante: 


Y 
=xt 4.3] 
k rs [4.31 
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Érant donné que la perméabilité intrinsèque, &, est uniquement fonction des caractéris- 
tiques du squelette solide du so, il est possible d'établir une corrélation entre le coelficient 
de conductivité hydraulique du sol à la résine, k, et le coefficient de conductivité hydrau- 
lique du sol à l'eau, 4 


BY 144] 
_ [AA 
. coefcient de conductivité hydraulique du sol à l'eau [m/s]: 
+ 4, = coefficient de conductivité hydraulique du sol à la résine [m/s]; 


30 mm, et des temps de réaction de la résine 
60 s, s; = 120, on a évalué le rayon du bulbe 
initial, @g) en fonction de la conductivité hydraulique à l'eau, 4. Les valeurs attendues 
pour les différents paramètres sont les suivantes: 

+ = 0,465 Pass; 

001 Pas; 

+ = 10,5 KN/m°; 

2 Va 9,8 KN/m?; 

+ n=03; 

01,3 kPa, 202,7 kPa. 

Les résultats obtenus dans le cas d'une cavité sphérique sont présentés sur Les figures 4.11 
à 4.16, tandis que ceux obtenus dans le cas d'une cavité cylindrique sont présentés dans 


les figures 4.17 à 4.22 (pages 43 à 46) 


CAVITÉ SPHÉRIQUE 
6 mm, p,= 101,3 kPa 


PotEG DOG ECS GNU EG OIGEE 1e 
ke Im/s] 


Figure 4.11 Diférence entre le rayon du bulbe initial et le rayon du trou di 
de la conductivité hydraulique du sol dans le cas d'une cavité sphérique (a = 6 mm) 
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Figure 4.12 Diflérence entre le rayon du bulbe initial et le rayon du trou 
de la conductivité hydraulique du sol dans le cas d'une cavité sphérique (a, 
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Figure 4.13 Diflérence entre le rayon du bulbe initial et le rayon du trou 
de la conductiité hydraulique du sol dans le cas d'une cavité sphérique (a; 
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CAVITÉ SPHÉRIQUE 
q=6 mm, p,= 202,7 kPa 
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Figure 4.14 Diflérence entre le rayon du bulbe initial et le rayon du trou 
de la conductivité hydraulique du sol dans le cas d'une cavité sphérique (a; 
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Figure 4.15 Diflérence entre le rayon du bulbe initial et le rayon du trou 
de la conductiité hydraulique du sol dans le cas d'une cavité sphérique (a, 


Modèles théoriques pour les sols à grains grossiers 


CAVITÉ SPHÉRIQUE 
4 = 30 mm, p, = 202,7 kPa 


1-30: 
a 
a 
Dal t=105 


IGGEGÉIGEGS GE LOG IE IRON 


KIm/s 


Figure 4.16 Diférence entre le rayon du bulbe initial et le rayon du trou d'injection en fonction 
de la conductivité hydraulique du sol dans le cas d'une cavité sphérique (a, — 30 mm). 
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Figure 4.17 Diférence entre le rayon du bulbe initial et le rayon du trou d'injection en fonction 
de la conducthité hydraulique du sol dans le cas d'une cavité cylindrique (ay = 5 mm). 
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CAVITÉ CYLINDRIQUE 
4=15mm, p, = 101,3 kPa 
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Figure 4.18 Diférence entre le rayon du bulbe initial et le rayon du trou d'injection en fonction 
de la conductivité hydraulique du sol dans Le cas d'une cavité cylindrique (a, = 15 mm). 
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Figure 4.19 Diférence entre le rayon du bulbe initial et le rayon du trou d'injection en fonction 
de la conductiité hydraulique du sol dans Le cas d'une cavité cylindrique (ay = 30 mm). 
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Figure 4 20 Diférence entre le rayon du bulbe initial et le rayon du trou d'injection en fonction 
dé la conductivité hydraulique du sl dans Le cas d'une cité cylindrique (a = 8 mm). 
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Figure 4221 Diférence entre le rayon du bulbe initial et le rayon du trou d'injection en fonction 
de la conductivité hydraulique du sal dans le cas d'une caié cylindrique (a = 15 mm). 
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Figure 4.22 Diférence entre le rayon du bulbe initial et le rayon du trou d'injection en fonction 
de la conductivité hydraulique du sol dans le cas d'une cavité cylindrique (a, — 30 mm). 


3.2. Détermination de la pression d'équilibre 


Dans la phase d'injection dans des sols à grains grossiers, la résine polyuréthane parvient 
à pénétrer dans les «pores» du sol à l'état liquide, en formant un bulbe de sol traité dont 
La forme dépend du mode d'injection. Dans le cas d'injections ponctuelles, le bulbe prend 
une forme sphérique, tandis que dans le cas de traitements «en colonnes, le bulbe se 
caractérise par une forme cylindrique. 

Dans l'étape ultérieure de solidification, la rés 
se dilatent jusqu'à ce qu'ils atteignent un état 


et, par conséquent, le bulbe de sol rraité 
quilibre avec le sol environnant. 


Au fur et à mesure de l'augmentation de la déformation radiale du bulbe, la pression de 
gonflement de la résine diminue, tandis que la contrainte radiale de confinement du sol 
environnant augmente. Le processus d'expansion s'arrête dès lors qu'est atteint un état 
d'équilibre où la pression de gonflement de la résine est égale à la contrainte radiale de 
confinement du sol. 


La déformation radiale du bulbe peur être évaluée au moyen d'un rapport sans dimension: 
a 
LA 


[4] 


olume initial de la résine à l'état liquide [m4]: 
+= volume final de la résine expansée [m?] 


La relation entre le rapport V-/V, et la pression de gonflement de la résine est représen- 
tée par la loi rhéologique indiquée dans le chapitre 1: 
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A 1 


Ve 036+— 
0,23 


où P est la pression de gonflement de la résine (exprimée en MPa). 


5 
In(l+ P(MPa)) 


Afin de déterminer la contrainte de confinement générée dans le sol, il est possible de faire 
appel à la théorie d'expansion des cavités dans les milieux élasto-plastiques. Cette théorie 
est exposée en détail dans l'annexe À, où sont présentés les développements théoriques 
mis au point par Yu & Houlsby (1991) pour tenir compte des grandes déformations ren- 
contrées dans la zone plastique. 

Le comportement d'un sol à grains grossiers en conditions limites (maximales) peut être 
assimilé à celui d’un matériau régi par un critère de rupture de Mohr-Coulomb d'un sol 
frottant et cohésif. En désignant par /la fonction régissant la déformation plastique, par 
0, la contrainte radiale et par @ la contrainte rangentielle, le critère de rupture peut être 
exprimé de la manière suivante: 


J=6,-N a -Y =0 1471 
n=1tsne (al 
[4.9] 
dans laquelle: 


+ = angle de résistance au cisaillement du sol; 


+ «= cohésion effective du sol. 


Le comportement d'un sol à grains grossiers dans un état critique peut être simulé en sup- 
posant un critère de rupture de Mohr-Coulomb pour un matériau purement frottant, où 
en posant Ÿ = 0 dans la formule [4.7]. 
Les paramètres géorechniques du sol utilisés son les suivants 
+ module de Young, £; et coefficient de Poisson, v, pour la caractérisation du compor- 
tement élastique du sol; 
+ cohésion, et angle de résistance au cisaillement, @; 
+ angle de dilatance, W. 
Pour obtenir la relation entr le rapport V,/V,et la contrainte de confinement générée 
dans Le sol, il est nécessaire, en premier lieu, de déterminer le rayon de La cavité, 4, sous 
l'action de la pression interne générique P. Le long de la limite de la cavité, la contrainte 
de confinement est égale à P, pour assurer l'équilibre. Le rayon de la cavité peur être déter- 
miné au moyen des équations suivantes: 
4 P£k 
P-P}ay 
2-m.G 
4-R P>R 


R<P<P 14.10] 
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où: 
+ 49 = rayon initial de la cavité ou du bulbe de sol traité; 


+ Pÿ= pression initiale, supposée égale à la contrainte géostatique horizontale dans le cas 
d'une cv 


é cylindrique, et à la contrainte géostatique moyenne dans le cas d'une 
cavité sphérique; 
+ P,= pression d'avant la plastification, évaluée par l'expression suivante: 


DA +8 11] 


+ m = coclicient géométrique, étant supposé égal à 1 dans le cas d'une eavité cylin- 
drique, égal à 2 dans le cas d'une cavité sphérique; 
Fee 

DETR 

+R, = rapport entre le rayon de la cavité, 4 à une pression générique > P,, et le rayon 
initial de la cavité, 


+ G= module de cisaillement, évalué comme G = 


Le paramètre R, peut être déterminé en se référant à la solution de Carter et al. (1986), 
basée sur l'hypothèse de petites déformations, soit à la solution de Yu & Houlsby (1991), 
basée sur l'hypothèse de grandes déformations dans la zone plastique et de petites défor- 
mations dans la zone élastique plus extérieure. 

Dans l'hypothèse de petites déformations, R, peut être déterminée de la manière sui- 
vante: 
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voir [4.9] 
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Si l'on considère l'hypothèse de grandes déformations dans la zone plastique, l'évaluation 
de &, est plus complexe et implique d'évaluer une intégrale ne pouvant être résolue dans 
un domaine fermé (Vu 8 Houlsby, 1991). Une solution approximative a été fournie par 
Yu & Houlsby (1991) sous la forme suivante: 


À nt [ 
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À partir du rayon final de la cavité, il est possible de dérerminer le volume final du bulbe, 
Wir 


(nd Hi cavité cylindrique 
| [421] 


re 
$ 


où H est la hauteur de la cavité cylindrique, égale à la hauteur de Ia portion de colonne 
dans laquelle est effectuée l'injection. 


sphérique 


Si l'on connaît le volume initial du bulbe, donné par: 
ma}: H cavité cylindrique 
= 4.22 
Len. & cavité sphérique al 
il est possible de déterminer la variation de volume subie par Le bulbe, en fonction de la 
pression P: 
AV = Vy Vis 14.23] 
En supposant que La variation de volume du bulbe soit uniquement imputable à l'expan- 
sion du volume de la résine, on peut déduire le volume de la résine finale, V;: au moyen 
de la relation suivante: 


a +AV, 14.24] 
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En connaissant le volume de résine injecté, il est donc possible de déterminer la relation 
entre le rapport VV, et la contrainte de confinement générée dans le sol, égale à la pres- 
sion dans la cavité ?. 

À titre d'exemple, prenons Le cas d’un sable sec, avec les paramètres géotechniques suivants : 
+ poids par unité de volume du sol, y, = 19 KN/m$; 

30 MPai 


+ module élastique de Young, 
+ coefficient de Poisson, v = 0,25; 

+ cohésion effective, c'= 0 kPa; 

+ angle de frotrement maximal, gp = 32° 
+ angle de frottement à l'état critique, Qey = 32°: 

+ angle de dilarance, 1} = 0°: 

+ degré de surconsolidation, OCR = 1 

On souhaite conceptualiser une injection de résine polyuréthane expansive de type Ure- 
tk Gcoplus” à une profondeur de 3,5 m; et on suppose la présence d'une surcharge, 
4= 10 KPa, au niveau de la surface naturelle du sol. En outre, on prend pour hypothèse 
que l'expansion de la résine se produit selon une cavité cylindrique (c’est pourquoi m = 1) 
ayant une hauteur 77 = 3 m; on suppose un rayon initial du bulbe, & = 0,1 m. Le 
volume initial de la résine, Y,, est égal à 0,019 m° (soit l'équivalent d'une masse initiale 
d'environ 20 kg, pour une porosité du sol » = 0,3). 


La figure 4.23 représente la loi rhéologique er la relation entre Le rappore VA, et la 
contrainte de confinement générée dans le sol, caleulée d'après la théorie d'expansion des 
cavités, dans l'hypothèse de petites déformations et de grandes déformations dans la zone 
plastique. 
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Figure 4.23 Rapport de gonflement de a résine en fonction de la pression de gonflement de La résine 
Qigne bleue) et de la contrainte de confinement du sol, dans l'hypothèse de petites déformations 
igne verte) et de grandes déformations ligne rouge) 


Modèles théoriques pour les sols à grains grossiers 


Comme on peur le voir sur la figure 4.23 en correspondance avec les valeurs du rapport 
VV légèrement supérieures à 1, la pression de gonflement est très élevée, tandis que la 
contrainte de confinement est proche de la contrainte géostatique horizontale, Il est, par 
conséquent, possible de déterminer un érat de non-équilibre qui détermine la déforma- 
tion uhérieure du bulbe, avec pour corollaire une augmentation du rapport V,/V,. Le 
point d'intersection entre les deux courbes correspond au point d'équilibre entre la pres- 
sion de gonflement de la résine et la contrainte de confinement exercée par le sol autour 
du bulbe. La valeur de pression correspondante est, par conséquent, la valeur d'équilibre 
final. 

Le point d'intersection entre les deux courbes est obtenu au niveau d'une pression d'équi- 
libre, P = 450 kPa, ce qui correspond à un rapport VV, & 5,3, dans le cas de petites 
déformations. 


Dans le cas de grandes déformations dans la zone plastique, on obtient P & 375 kPa, ce 
qui correspond à un rapport V,4V,,= 6. 

Il convient de noter que, si l'on néglige l'hypothèse des grandes déformations dans la zone 
plastique, on à tendance À surestimer la pression d'équilibre et à sous-estimer le taux de 
gonflement, 


3.3. Variation de l'état de contraintes dans le sol 


Après avoir déterminé a valeur de la pression d'équilibre au moyen de la procédure décrite 
dans le paragraphe précédent, il est possible, avec la théorie de la cavité en expansion, 
d'évaluer la modification de l'éar de contraintes induit dans le sol entourant le bulbe. 


L'état de contraintes dans le domaine élastique extérieur au rayon plastique (c'est-à-dire 
pour r > b, où best le rayon plastique) peut être exprimé au moyen des relations suivantes 


: a jen n 23 
m+N 
Libre (vint : 
S=R Le v L 426) 


Dans la région de plastifcation (c'est-à-dire » < 4), la contrainte radiale er la contrainte 
tangentielle peuvent être calculées avec les équations suivantes: 


C+m)N[Y+(W-1}27] te 
CT 127) 
(+m}N [Y+(N 0-2] 

CERN) en [428] 


Dans l'hypothèse de perires déformations, le rayon plastique, d, est obtenu avec l'équation 
suivante: 


Re 14.29] 
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Dans le cas où l'on prend pour hypothèse de grandes déformations dans la zone plastique, 
best donné par: 


ru l430] 
n 
où a est Le rayon du bulbe après que l'expansion a eu lieu. 
À partir de La détermination de la pression d'équilibre, ?, il est possible, par conséquent, 
d'évaluer la modification de l'état de contraintes dans Le sol. Dans le cas de la cavité cylin- 
drique, l'augmentation de la contrainte verticale, AG, peut être finalement évaluée 
d'après La relation élastique: 


A6, = VAG, + AG,) [4.31] 


3.4. Variation de l'état de déformation du sol 


L'expansion du bulbe provoque, en plus de l'augmentation de La contrainte moyenne, une 
déformation volumétrique positive (c'est-à-dire de compression) à l'intérieur de la zone 
plastique. L'évaluation de cetre déformation peut être réalisée de manière analytique, en 
se référant uniquement à l'hypothèse de petites déformations. Dans ce cas, la déformation 
volumétrique &, est égale à: 


tm 432 


où: 


+6, = déformation radiale, donnée par: 
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+ & = déformation tangentielle, donnée par: 
se fs, NME _ Mn) |, 1] 134] 
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En se servant des formules [4.32] à [4.34], il est possible d'évaluer la diminution de 
volume qui se produir à l'intérieur de la zone plastique pour une valeur de pression de la 
cavité expansée égale à la pression d'équilibre. À partir de ce résultat, on peut déterminer 
l'augmentation de la densité relaive, D), déterminée par l'expansion du bulbe. 


La densité relative est fonction de l'indice des vides selon la relation suivante: 


D = 14.35] 
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et où épinet Fm représentent respectivement les indices minimum er maximum des vides 
du sol. La variation de l'indice des vides peut être corrélée à la déformation volumétrique, 
£, au moyen de la relation suivante 


[4.36] 
14e, 


Étant donné que la variation de la densité relative, AD, peut être corrélée à la variation 
de l'indice des vides, Ae, au moyen de La relation suivante: 


1437) 
on obtient: 


[4.38] 


où es l'indice des vides initial 


[4 Modèles théoriques pour sols à grains fins 


Le processus d'expansion des résines polyuréthanes injectées ponctuellement dans des sols 
à grains fins (sables fins, limons et argiles) se prête à être classé théoriquement suivant l'ap- 
proche de la théorie de la fracturation dans un milieu élasto-plastique, dans des condi- 
ions non drainées. Le processus de fracture apparaît lorsque la résine, en ne parvenant pas 
à pénétrer dans Les pores du sol, développe une pression de gonflement très élevée au cours 
de la phase d'injection. 

Les résines Ureték se propagent à l'intérieur des fractures préexistantes dans le sol et dans 
celles générées par l'expansion des résines proprement dites lors de l'injection. 

L'étude de la propagation des fractures dans le sol est d'une complexité considérable, 
assortie de difficulkés de prévision, parce quelle est conditionnée par la présence d'un 
manque d'homogénéité dans le sol. On peut seulement émettre l'hypothèse que la résine 
injecrée pénètre dans les microfissures qui se sont formées au niveau des zones où se 
concentrent les déformations de cisaillement (bandes de cisaillement), ce qui provoque 
leur propagation dans le sol. 


4.1. Critères d'apparition du phénomène de fracturation 


Dans un sol cohésif, l'injection de résines expansives Uretek provoque la formation de 
fractures dans le matériau, Le premier résultar évaluable en théorie est l'ampleur de la 
pression de gonflement que produit l'apparition du processus de fracturarion. Cette pres- 
sion est appelée «pression de fracturationv, 

Du point de vue macroscopique, le déclenchement d'une fracture peur être attribué à des 
mécanismes de rupture en traction ou en cisaillement (figure 4.24). 
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Figure 4.24 Fradturation du sol provoquée par l'injection d'un fluide sous pression dans une cavité 
@) rupture en traction dans des conditions non dainées et () rupture imputable à a résistance 
au cisaillement dans des conditions non drainées (modifié par Mitchell & S0ga, 2005). 


La rupture en traction se produit quand la contrainte efficace inférieure s'annule ou est 
égale à la résistance à la traction du matériau, &; quand elle est différente de zéro. 

Dans le cas d'une cavité sphérique ou cylindrique, cette rupture se produit dans des 
conditions élastiques, parce que là contrainte tangentielle efficace, 0’8, diminue avec 
l'augmentation de la déformation radiale de la cavité, jusqu'à ce que soit atteinte la résis- 
tance à la traction du sol (figure 4.24a). 

La fracturation est, au contraire, générée par la rupture en cisaillement quand, avant 
même que la contrainte rangentielle arteigne La valeur de la résistance à la traction, la rup- 
ture en cisaillement (figure 4.24b) est déterminée. 

D'après les indications de Mitchell & Soga (2005), à partir des Lois d'élastcité linéaire, on 
arrive à l'expression analytique suivante concernant la pression de fracturation selon le 
mécanisme de rupture en traction: 


Péts (14m): 0, mu +m 0! 14.39] 
où 
+ Pépusion = Pression à l'intérieur de la cavité qui conduit à la fracturation en traction 
sol; 
+ m = coefficient géométrique étant supposé égal à 1 dans le cas d'une cavité cylin- 
drique, égal à 2 dans le cas d'une cavité sphérique; 
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+63; = contrainte principale totale inférieure à l'état initial; 


+ 4 = pression intersitielle initiale en conditions hydrostatiques; 


+ o= résistance à la traction du sol. 


Cette équation donne une relation linéaire entre la pression de confinement et la pression 
de fracturation, avec un cocfficient de proportionnalité égal à deux dans le cas d’une 
cavité cylindrique er égal à trois dans le cas d'une cavité sphérique 

Dans Le cas d'un mécanisme de fracturation en cisaillement, Le procédé peut être déclen- 
ché par la formation de microfissures locales sur le contour de la cavité en expansion. La 
résine peut s'infltrer dans Les fractures er engendrer des efforts de traction aux extrémités 
de la fracture proprement dite (figure 4.24b). Dans l'hypothèse où le critère d'analyse du 
début de la fracturation correspond à la déformation plastique initiale de la paroi d'une 
cavité cylindrique en expansion, on peut obtenir la relation suivante rapportée à la pres- 
sion de fracturation pour le mécanisme de rupture en cisaillement: 

2m 


EM 14.40] 
m+l 


+ Pjésaitms = Pression à l'intérieur de la cavité qui conduit à la fracturation en cisaille- 


ment du sol, correspondant à la pression d'avant là plastification; 


+6, = cohésion en condition non drainée du sol. 


L'équation [4,40 à la différence de l'équation [4.39], implique une relation entre la pres- 
sion de confinement et la pression de fracturation linéaire avec un cocfficient de propor- 
dionnalité unitaire, dans le cas d'une cavité cylindrique, et égal à 4/3 dans le cas d'une 
cavité sphérique. 

Dans la région de plastification, il est possible d'avoir, localement, une concentration de 
déformations de cisaillement et que, par conséquent, des microfissures apparaissent dans 
le matériau. Dans ces conditions, il est raisonnable de supposer que les injections de résine 
Uretek provoquent l'ouverture des microfractures et favorisent leur propagation dans le 
sol traité. 

Les expressions analytiques [4.39] et [4.40] permettent d'établir, en première approxima- 
ion, si l fracturation aux alentours de la cavité est une facturation en traction ou en 
cisaillement. 


Du point de vue théorique, il est possible d'établir, en première approximation, les condi- 
tions de contrainte dans lesquelles le système de fracturation se révèle sub-verrical ou sub- 
horizontal. La fracturation se produit en direction sub-verticale si les contraintes 
horizontales, dj, sont inférieures aux contraintes verticales effectives, 019, dans des 
conditions géosratiques, c'est-à-dire si le coefficient de pression des terres au repos Ky< 1 
En situation inverse, c'est-dire si Kj, > 1, le système de facturation est sub-horizontal. 
L'extension, la fréquence er l'ouverture des fractures sont difficiles à prédire, d'un point de 
vue théorique, parce qu'elles dépendent de la présence de mierodéfauts er de fissures dans 
le sol. 
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4.2. Critères de propagation des fractures 


Si l'injection est réalisée dans le sol de surface, les conditions de diffusion des fractures 
(par rupture en traction ou en cisaillement) sont facilement vérifiées. Dans ce cas, si au 
point d'injection la contrainte horizontale est inférieure à la contrainte verticale (K, < 1), 
il se forme une fracture sub-verticale; inversement, si la contrainte horizontale est supé- 
rieure à la contrainte verticale (K, > 1), il se forme une fracture sub-horizontale. 
Les fractures, en première approximation, peuvent prendre la forme d'une pièce de mon- 
naîe («penny-shaped fracture, en anglais). Cette rypologie de fractures a été étudiée en 
théorie par Sneddon 8 Lowengrub (1969), qui ont obtenu des résultats utiles pour l'éva- 
luation de leur longueur er de leur ouverture. 
La longueur de la fracture dépend, fondamentalement, de critères énergétiques (Griffith, 
1920) ou mécaniques (Irwin, 1948). Ces critères sont difficilement applicables en pré- 
sence de matériaux non homogènes comme les sols qui présentent des propriétés phy- 
siques et mécaniques variables. En se référant au critère énergétique de Griffith, la 
demi-longueur de la fracture, £, peut être associée à la pression interne, P, au moyen 
d'une relation qui est fonction de l'énergie de surface du matériau, G, (Sneddon & 
Lowengrub, 1969): 

nE-G, 
2-7 


l441] 


où Æ est le module de rigidité du matériau. 


En variante, en adoptant le critère mécanique d'Irvin, la demi-longueur L est égale à: 
LT 


14.42] 


ap? 


où &, représente la ténacité à la fracture (en anglais, « fracture toughness»), dont la valeur 
est fonction de la résistance mécanique du matériau. 
Les relations [441] et [4.42] sont équivalentes, parce qu'il peut être démontré que 
e-1t=v Ke 
Sr he 
Dans le cas des sols, il existe peu de données de la littérature concernant les valeurs de 
l'énergie de surface, G, et de la ténacité à la fracture, &,, dans les sols. Ces paramètres 
restent donc très incertains et dificiles à déterminer. Fang (1994) rapporte les résultats 
d'une série d'essais effectués sur des échantillons d'argile compactée afin de déterminer la 
ténacité à la fracture, Dans ces essais, pour des sols différents de par leur teneur en eau et 
leurs degrés de compactage, la ténacité à la fracture a toujours été supérieure à 


14.43] 


0.08 MPa nr. Dans tous les cas, ces données à elles seules ne sont pas suffisantes pour 
tirer des conclusions sur Les valeurs pouvant êrre atteintes dans tous les sols. Il s'ensuit que 
la détermination de la longueur des fractures est us incertaine er évaluable uniquement 
en se fondant sur les résulrats des essais sur le terrain dans lesquels sont effectuées des exca- 
vations appropriées afin de reconstruire la géométrie de la propagation de la résine. 
Après avoir fixé la longueur de la fracture, il est possible d'identifier la pression d'équilibre 
au moyen des relations fournies par Sneddon & Lowengrub (1969) pour le calcul de l'ou- 
verture de la fracture, en fonction de la pression interne P dans un milieu purement élas- 
tique. En particulier, le rapport de gonflement est égal à: 
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1 sP<o, 


el mpica-vrs 

VE po) | sP> 04 
NAT: 

En condition non drainée, on suppose que v = 0,5 et E = Æ, où Fest le module de rigi- 

dité en condition non drainée. 


[4.44] 


À partir de l'intersection entre la relation [4.44] er la loi rhéologique [4.6], il est possible 
de déterminer la pression d'équilibre pour une demi-longueur spécifique de la fracture. 
Par exemple, supposons que l'on fasse une injection à une profondeur de 3,5 m dans une 
saturée, caractérisée par les paramètres physiques et mécaniques suivants : 

18,5 KN/m°; 


La contrainte verticale est égale à 64,75 kPa. S'agissant d'argile normalement consolidée, 
Le coefficient de pression des terres au repos peut être caleulé au moyen de la relation sui- 
vante: 


1445] 


ce qui donne &, = 0,577. 

Étant donné que la valeur de K, est inférieure à l'unité, on peut s'attendre à la formation 
d'une fracture sub-verticale. En supposant une demi-longueur de la fracture égale à 
0,25 m, on obtient le résultat de La figure 4.25. 
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Figure 4.25 Rapport de gonflement de a résine, en fonction de la pression de gonflement (igne bleue) 
et de la contrainte de confinement exercée par le sol environnant (igne rouge). 
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La pression d'équilibre est d'environ 570 kPa. En supposant par exemple: 
K je = 0,08 MPa: Un 


on obtient à partir de l'équation [4.42] une pression critique, associée à la demi-longueur 
de la fracture L = 0,25 m, soit: 


14.46] 


Dans ce cas, ou pour Ki = 0,08 MPa:  , les résultats ne sont pas concluants. Il 
convient, par conséquent, de prendre pour hypothèse une plus grande longueur de la 
fracture, pour rechercher la correspondance entre la pression d'équilibre er la pression eri- 
tique. Si l'on suppose une demi-longueur de La fracture £ = 0,67, on arrive au résultat de 


La figure 4.26. 
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Figure 4.26 Rapport de gonllement de a résine, en fonction de la pression de gonflement (ligne bleue) 
et de la contrainte de confinement exercée par le sol environnant (ligne rouge). 


La pression d'équilibre est, dans ce cas, d'environ 124 kPa. Pour Ke = 0,08 MPa: 
on obtient une pression critique égale à: 


Kk, fr 
pee Mio, =124kP 4 
PET +04 APS 14.47] 


On peut donc en conclure que pour l'argile considérée, en supposant une ténacité à la 


fracture Ke = 0,08 MPa -Vm , la demi-longueur L = 0,67 correspond à la valeur exacte, 
permettant la congruence de la pression critique dérivée du critère d'Irwin avec la pression 
d'équilibre. 
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4.3. Soulèvements induits par des injections de résine 


Pour produire le soulèvement du niveau de la surface naturelle du sol et des éventuelles 
superstructures qui se trouvent au-dessus du niveau des injections, il convient de détermi- 
ner la formation de fractures sub-horizontales. Même dans les sols caractérisés par des 
valeurs de coctficient de pression des terres au repos, Kg inférieures à l'unité, la formation 
initiale de fractures sub-verticales suivie de leur ouverture produit, localement, une inver- 
sion des contraintes principales, qui permet a formation ultérieure des fractures sub-hori- 


zontales (figure 4.27). 
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Figure 4.27 Formation de fractures sub-horizontales dans les sols ayant un coefficient de pression des terres 
au repos K, < 1 (phase 1): l'ouverture d'une fradture verticale provoque l'augmentation des contraintes 
horizontales, entrinant une rotation des directions principales (phase 2); l'exécution ultérieure d'une 
nouvelle injection détermine La formation d'une fracture sub-horizontale, dont l'ouverture produit 

le soulèvement du niveau de la surface naturelle du sol 


À proximité du point d'injection, l'ouverture d’une fracture verticale produit une aug- 
mentation de la contrainte horizontale du sol, égale à la pression d'équilibre nette. 
Si l'on prend pour référence le dernier exemple du paragraphe précédent, les contraintes 
géostatiques valent: 
+ 6,0 = 65 kPa 
+ 09 = 37 kPa 
La pression d'équilibre netre est égale à: 

P 09 = 124 — 37 = 87 kPa. 
I s'ensuit qu'à proximité du point d'injection, la contrainte horizontale est égale à la pres- 
sion d'équilibre, ou: 


0,= P= 124 kPa > Gg. 
Les prévisions quant à la longueur de la fracture horizontale peuvent être déterminées 
avec les mêmes critères que ceux décrits précédemment, ou en recherchant l'état d'équi- 
libre entre la pression de gonflement de la résine et la contrainte de confinement (dans ce 
cas, verticale) générée dans le sol. 

L'ampleur des soulèvements, si l'on omet la compression de la couche de sol au-dessus du 
point d'injection, peut être estimée, dans un premier temps, au moyen de la formule éta- 
blie par Sneddon & Lowengrub (1969) pour le calcul de l'ouverture de la fracture 


4-Q-v)r dE - 
Fa 


ax) = E 


14.48] 
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où: 
+= distance à partir du point d'injection. 

Le déplacement maximal est obtenu au niveau du point d'injection, ou pour x = 0. 
Pour l'application de l'équation [4.48], les incertitudes liées à l'identification de la demi- 
longueur L de La fracture sont confirmées, En outre, en fournissant une estimation pour 
une telle grandeur, le recours aux solutions d'analyse ne permet pas de renir compte de la 
géométrie réelle du problème. De ce point de vue, les solutions numériques permettent 
d'obtenir des résulrats plus précis et plus exacts. 


[5 Conclusion 


Les injections de résines expansives polyuréthanes sont différentes des autres rypes d'in- 
jections, dans la mesure où elles sont à même de développer spontanément une pression 
de gonflemen élevée, qui atteint un état d'équilibre avec la contrainte de confinement 
exercée par le sol environnant. 

Dans les sols à grains grossiers, les injections entraînent la formation d'un bulbe de sol 
traité, dont le gonflement dérermine la compaction et la densification du sol environnant. 
En outre, le bulbe proprement dit constitue un élément de renfort, susceptible de faire 
augmenter la résistance au cisaillement et la rigidité du sol traité. 

La mise en place de conditions non drainées provoque, en revanche, la formation de frac- 
rures dans les sols à grains fins. L'ouverture des formations horizontales peut être exploi- 
tée pour compenser les affaissements complets et/ou différentiels des fondations 
superficielles. En outre, le sol est d'autant plus densifié que le réseau des injections effec- 
ruées est dense, Puisque la dissipation des surpressions interscitielles, au fil du temps, pro- 
duit une modification de l'état de contraintes, il convient de tenir compte, dans ces cas, 
de la nécessité éventuelle de procéder à des injections répétées pour atteindre les objectifs 
fixés. Une évaluation de la variation des élévations induites par la dissipation des pressions 
intersticielles ainsi que de la durée du processus de consolidation ne peut se fonder sur des 
modèles analytiques simplifiés, tels que ceux décrits dans le présent chapitre. Le cha- 
pitre 5 ci-après présente un exemple d'analyse numérique aux éléments finis concer- 
nant le processus de fracturation d'un sol argileux, dans lequel sont déterminées les 
surpressions interstitielles qui se forment au cours du processus d'expansion de la résine, 
ainsi que les déplacements induits par leur dissipation. Ce n'est que par une analyse 
numérique aux éléments finis qu'il est possible d'évaluer avec précision le processus de 
consolidation er les déplacements qu'il produit. 


Introduction 


Ce chapitre décrit certaines analyses cffecruées au moyen du logiciel Plaxis 2D afin de 
simuler les effets induits dans le sol par des injections localisées de résines expansives, 
Les analyses effectuées peuvent être considérées comme des références, par exemple, pour 
la modélisation aux éléments finis des interventions effectuées au moyen d'injections de 
résines polyuréthanes afin d'améliorer les caractéristiques mécaniques du sol (compacrion 
grouting, ou injection solide) ou de reprendre en sous-œuvre des tassements excessif de 
fondations superficielles (compensarion grouting, ou injection de compensation) 
Le logiciel Plaxis 2D est un logiciel de caleul aux éléments finis à deux dimensions, per 
mettant de simuler le comportement du sol en fonction des variations de l'état de 
contraintes et de déformations en différents points du sol et dans les éventuels éléments 
de structure avec lesquels il interagir IL utilise habituellement pour effectuer des analyses 
de déformation et de stabilité dans le cadre de multiples applications géoechniques et 
permet de simuler des situations réelles se rapportant à des déformations planes (plane 
strain) ou à des déformations axisymétriques (axisymmetrie). 

Le comportement mécanique du sol peut être simulé au moyen des modèles de compot- 

tement suivants: 

— Modèle linéaire élastique — Simule le comportement d'un milieu élastique linéaire iso- 
wope. Ce modèle implique de définir deux paramètres de déformabilité dans Le 
domaine élastique: le module d'Young (F) et le coefficient de Poisson (v) 

— Modèle de Mohr-Coulomb — Simule le comportement d'un milieu élastique linéaire 
isotrope parfaitement plastique. En plus des paramètres de déformabilité dans le 
domaine élastique, F et v', il est nécessaire de définir des paramètres de résistance, à 
savoir: cohésion effective (e), angle de frottement (q'} et angle de dilatance (4) 


— Modèle Jointed-Rock — Simule Le comportement d'un milieu élasto-plastique aniso- 
trope le long des plans de discontinuité. 

— Modèle Hardening-Soil — Modèle hyperbolique de type élaste-plastique permettant 
d'examiner les variations de la rigidité du sol en fonction de l'état de contraintes ct du 
degré de consolidation, 

— Modèle Hardening-Soil Small - Modèle analogue au modèle Hardening-Soil permet 
tant de différencier le comportement du sol en petites déformations de son comporte- 
ment en grandes déformations. 

— Modèle Cam-Clay modifé — Modèle simplifié adopté pour simuler le comportement 
des sols cohésifs normalement consolidés. 
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— Modèle Soft Soil (Cam-Clay) - Modèle permettant de simuler le comportement de 
sols ou de matériaux très compressibles tels que Les argiles normalement consolidées ou 
la tourbe. 

— Modèle Soft Soil Creep — Modèle analogue au modèle Soft Soil permettant d'exami- 
ner le comportement visqueux du matériau. 

Le logiciel Plaxis permet d'effectuer différents types de calculs par la méthode des élé- 

ments Anis, en faisant la distinction entre les analyses de type plastique (plastic), les ana- 

Iyses de consolidation (consoldarion), les analyses de stabilité avec la méthode de réduction 

des paramètres de résistance (ph, « reduction) et les analyses dynamiques (dynamic), en 

conditions drainées et non drainées. 

Le logiciel Plasis a été utilisé pour simuler les effet liés à l'expansion de la résine polyuré- 

thane dans le sol qui impose une déformation volumétrique au niveau du cluster repré- 

sentatif de la zone occupée par le volume initial du Auide injecté. 

La version du logiciel adoptée pour les analyses réalisées dans le présent document est 

Plaxis 2D v 8.6. 


L. Interventions par injection solide 


Ainsi qu'il a été mentionné dans les chapitres précédents, les interventions par injection 
solide sont généralement effectuées dans le bur d'augmenter l'état de compaction du sol 
entourant le ou les points d'injection et, par conséquent, d'améliorer les propriétés méca- 
niques du matériau. 

Dans le présent paragraphe, le logiciel Plaxis 2D a servi à simuler certaines applications 
possibles de ce type d'interventions, aux fins d'évaluer la capacité du modèle numérique 
À reproduire les effets induits dans le sol environnant par l'expansion des résines expan- 
sives de polyuréthane injectées localement, 


2.1. Expansion de la cavité sphérique 


La première analyse réalisée à l'aide du logiciel Plaxis 2D dans Le cadre d'interventions par 
injcetion solide prévoyait d'effectuer une simulation des effes liés à une injection ponctuelle 
de résine polyuréthane dans un milieu granulaire homogène isotrope en l'absence de nappe 
(C'est-à-dire, sable meuble où légèrement compact). Le comportement mécanique de ce 
marériau a éé mis en œuvre à l'aide du modèle de comportement de Mohr-Coulomb, dont 


les paramètres caractéristiques supposés dans l'analyse sont présencés dans le tableau 5.1. 


Tableau 5.1 Paramètres géotechniques atwibués dans l'analyse. 


Paramètres de déformabilité Paramètres de résistance 
Module Cohésion Angle de 
de Young effective, c dilatance, 
IMPal kPa] ml 
Sable meuble ou à Pen . à = 
légèrement compact 


Interventions par injection solide 


On suppose que le bulbe d'injection, généré du fair de la pénétration de la résine à l'état 
liquide dans les vides incersitiels entre les grains de sable, est de forme sphérique et a un 
rayon de 20 em à compter du point d'injection. La loi de comportement pour un matériau 
de ce type est la même que celle d'un sable meuble ou légèrement compact (tableau 5.1). 
Le modèle géométrique de référence, représenté sur la figure 5.1, est un modèle axisymé- 
urique et a été créé avec une surface celle qu'il soit possible d'ignorer les «effets de bord 
sur les résultats de l'analyse (largeur = 10 m, hauteur - 10 m). 


Figure 5.1 Modèle géométrique de référence pour l'analyse. 


Le point d'injection de la résine a été localisé à une profondeur de 3 m à partir du niveau 
de la surface naturelle du sol. 
Afin d'évaluer la pression de gonflement caractéristique de la résine, on a d'abord com- 
mencé par définir la courbe de réaction du sol au processus d'expansion de la résine pro- 
prement dite, Cette courbe a été obrenue en appliquant au cluster du matériau représentant 
le bulbe d'injection diverses valeurs de déformation volumérrique, pour chacune des- 
quelles le logiciel de caleul a permis de calculer la valeur de là pression de confinement 
exercée par la formation sableuse au niveau de la surface de gonflement. À partir de la 
déformation volumétrique définie par l'utilisateur, il est possible de caleuler la variation 
de volume (AV) due à l'expansion de la résine au moyen de la relation suivante: 

AV = A, V, 51] 


déformation volumétrique définie par l'utilisateur; 


Modèles numériques 


+ Vi = volume initial du cluster représentant Le bulbe d'injection 


En 


3 0.033510 m° [52] 


avec: 
— = rayon du bulbe d'injection, dans le cas présent supposé égal à 0,2 m. 
Le volume initial de la résine injectée (V,) peur être calculé comme suit: 


0.010053 m° =10 litres 153] 


+ Vi = volume initial du cluster représentant le bulbe d'injection 

+ n = porosité totale du sable meuble ou légèrement compact, dans le cas présent sup- 
posé égal à 0,3. 

Enfin, pour chaque valeur de déformation volumétrique définie par l'urilisateur, il est 

possible de déduire le volume final de la résine à la suite du processus d'expansion (V.): 


Ve = Vi +4V 154] 


avec: 
+ V,= volume initial de résine injectée; 


+ AV= variation de volume due à l'expansion de la résine, 


ÈË 
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Pression de confinement, P MPa] 
Réaction du sol au processus d'expansion 


Figure 5.2 Courbe de réaction du sol au processus d'expansion de la résine polyuréthane, 
caleulée par le logiciel Plxis 2D. 


La comparaison du résultat obrenu avec la solution analytique proposée par Carter et al. 
(1986) (figure 5.3) montre une correspondance optimale, déterminée par la congruence 
entre les conditions aux limites adoptées dans les deux solutions. Dans ee cas, les prévi- 
sions effectuées au moyen de la solution analytique coïneident avec celles obtenues par 
l'approche de type numérique aux éléments finis. Cette dernière, toutefois, offre la possi- 
bilité d'envisager des géométries plus complexes et des propriétés du sol non homogènes. 
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Pression de continement, P [MPa] 


+=" Résction du sol au processus d'expansion 
— Solution analytique proposèe par Carter et al. (1986) 
Figure 5.3 Comparaison entre la courbe de répanse du sol au processus d'expansion de la résine 


polyuréthane, calculée par le logiciel Plais 2D, et celle définie à partir de la Solution anabtique 
proposée par Carter et al. (1986). 


Une fois connue la réaction du sol au processus d'expansion de la résine polyuréthane, il 
a été possible de définir la pression de gonflement caractéristique de la résine proprement 
dite pour le traitement étudié (7), en identifiant le point d'intersection entre la courbe 
décrire ci-dessus er la rhéologie propre au Auide injecté (figure 5.4). 


| 


0 O1 02 03 64 05 06 07 08 0$ 1 1 12 13 1 1 
Pression de confinament et de gontloment, P [MPa] 


_—+- Réaction du sol au processus d'expansion 
_— Courbe théologique de la résine polyuréthane 


Figure 5,4 Définition de a pression de gonflement caractéristique de la résine et du rapport de gonflement 
pour le traitement objet de l'étude. 


Modèles numériques 


À partir de la valeur de P. obtenue, nous avons procédé à l'évaluation de la déformation 
volumétrique correspondante à appliquer dans la modélisation numérique afin de cale 
ler les variations de l'état de contraintes et de déformations induites par le processus d'in- 
jection et d'expansion dans le sol environnant 


En premier lieu, la loi rhéologique de la résine a permis d'évaluer le rapport de gonfle- 
ment caractéristique de l'inrervention (Ac) associé à la pression caractéristique Pc 
figure 5.4), au moyen duquel on a calculé le volume final de la résine à la suite du pro- 
cessus d'expansion (V,) + 

2 10,5 


0364 n(+ 2) 
0,23 


83 15.5] 


où: 
+ PC = pression de gonflement caractéristique de la résine pour le traitement étudié 


[MPa], égale à 0,73 MPa. 
V = Re V, = 0038503 m° = 38 5 litres (5.61 


+ V,= volume initial de résine injecrée. 
Par la suite, on a calculé la variation de volume de la résine survenue pendant le processus 
d'expansion (AW), à partir de laquelle il a été possible de dériver la valeur de la déforma- 
don volumérrique à introduire dans l'analyse numérique (A): 


AV = VV, = 0.028450 m° 1571 


+ V= volume final de la résine à la suite du processus d'expansion; 
+ V,= volume initial de résine injectée, 


av 
7, 849% 15.8] 


avec: 
+ Vi = volume initial du cluster représentant Le bulbe d'injection. 

Les figures 5.5 à 5.10 présentent Les résultats obtenus à partir de l'analyse numérique 
mise en œuvre à l'aide du logiciel Plaxis 2D, lesquels mettent en évidence les déplace- 
ments totaux subis par Le sol entourant le bulbe d'injection de la résine du fait du proces- 
sus d'expansion, les variations de l'état de contraintes et les limites de la formation 
sablonneuse satisfaisant aux conditions de plastification 


En ce qui concerne les déplacements totaux, on constate qu'ils se produisent principale- 
ment dans la direction radiale par rapport au point d'injection et affectent la formation 
sablonneuse jusqu'à une distance maximale d'environ 1 m par rapport au point d'injection 
proprement dit. Le déplacement total maximal mesuré se situe au niveau du périmètre exté- 
rieur du bulbe d'injection, et est égal à 59,8 mm. Il est important d'observer que la géomé- 
trie du modèle influence considérablement les résultats obtenus. En effet, les déplacements 
induits montrent que le traitement à une plus grande influence sur la partie du sol située 
au-dessus du point d'injection, en raison de la présence de contraintes de confinement infé- 
rieures à celles qui agissent sur le matériau situé à des profondeurs plus importantes. 
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1 KES 


Figure 5.6 Déplacements totaux du sol (arrow, fes»). 


Modèles numériques 


Figure 5.8 Contraintes verticales (shadings, « dégradés). 


Interventions par injection solide 
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Figure 5.10 Points de plastfication du sol. 


Modèles numériques 


En ce qui concerne l'éat de contraintes, les résultats obtenus permettent de constater que 
le sol entourant le point d'injection subit une variation des contraintes horizontales 
jusqu'à une distance maximale égale à environ 1,5 m, tandis que l'effet de la variation des 
contraintes verticales est épuisé à environ 1 m du point d'injection proprement dit. En 
outre, de manière similaire à ce qui a été observé pour les déplacements, le sol au-dessus 
du bulbe est davantage affecté par la variation des contraintes verticales que celui situé au- 
dessous de lui. De plus, il est clair que la pression de gonflement exercée par la résine sur 
les parois du bulbe d'injection n'est pas parfaitement uniforme, ainsi qu'il est supposé 
dans la solution analytique proposée par Yu & Houlsby (1991), mais qu'elle varie en fonc- 
ion de l'érat de contraintes de confinement du sol environnant: les valeurs maximales 
s'observent, en effet, dans la partie inférieure du bulbe et diminuent progressivement au 
fr et à mesure que l'on se déplace vers le haut. 

Enfin, sur la figure 5.10, on peut voir que les conditions du traitement analysé sont de 
nature à provoquer la réalisation des conditions de plastfication dans une zone de sol 
assez étendue, qui va jusqu'au niveau de la surface naturelle du sol, ce qui pourrait se révé- 
ler critique en prévision de la réponse sous forme de déformation du matériau à d'éven- 
ruelles variations de charge induites en surface. 


2.2. Expansion de la cavité cylindrique 


Le second exemple d'intervention par injection solide, simulée à l'aide du logiciel Plaxis 
2D, a pour objectif d'effectuer la simulation d'un traitement «en colonne» d'un milieu 
granulaire homogène et isotrope en l'absence de nappe (c'est-à-dire sable meuble ou légè- 
rement compact) 

Le modèle géométrique de référence et les caractéristiques mécaniques attribuées aux 
matériaux sont similaires à ceux de l'analyse décrite dans le paragraphe précédent, excep- 
tion faite de la géométrie du bulbe d'injection, pour lequel on a émis l'hypothèse qi 
était cylindrique, d'une hauteur égale à 2 m er avec un rayon égal 0,2 m (figure 5.11). 
Dans ce eas-là aussi, le point d'injection est situé à une profondeur moyenne de 3 m par 
rapport au niveau de la surface naturelle du sol, à l'origine 


Par analogie avec ce qui a été effectué dans l'analyse présentée dans le paragraphe précé- 
dent, on a procédé à l'évaluation de la pression de gonfement caractéristique de La résine 
pour l'intervention étudiée. 

En premier lieu, on a défini la courbe de réaction du sol vis-à-vis du processus d'expan- 
sion de La résine polyuréthane, calculée à l'aide du logiciel Plaxis 2D (figure 5.12) 

Il ressort du résultat obrenu que, dans l'analyse, la pression limite de la formation sablon- 
neuse, de l'ordre de 0,45 MPa, a été atteinte, ce qui correspond au seuil au-delà duquel la 
résine poursuit son processus d'expansion sans qu'apparaissent dans le matériau environ- 
nant de nouvelles augmentations de la pression de confinement. 


Interventions par injection solide 


Figure 5.11 Modèle géométrique de référence pour l'analyse. 
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Figure 5.12 Courbe de réaction du sol au processus d'expansion de la résine polyuréthane, 
calculée par le logiciel Plaxis 2D. 


Modèles numériques 


En comparant le résultat obtenu avec la solution proposée par Carter er al. (1986) 
(figure 5.13), on peut observer une bonne congruence jusqu'à des pressions de l'ordre de 
0,2 MPa. Pour des valeurs de pression plus élevées, la solution numérique tend à se diffé- 
rencier de la solution analytique, notamment en ce qui concerne la géométrie de la zone 
de plastification, qui, dans la solution numérique, est affectée par la présence du niveau 
de la surface naturelle du sol. 


0 w 0 o7 
Pression de confinement, P [MPa] 


| Réaction du sol au processus d'expansion 
_— Solution analytique proposée par Carer et al. (1986) 
Figure 5.13 Comparaison entre la courbe de réaction du sol au processus d'expansion de La résine 


polyuréthane, calculée à l'aide du logiciel Plaxis 2D, et celle définie à pari de l solution anahtique 
proposée par Carter et al. (1986). 


Enfin, l'intersection entre la courbe de réaction du sol au processus d'expansion de la 
résine polyuréthane er la courbe théologique du Auide injecté a permis de définir à la fois 
la pression caractéristique (P4) et le rapport de gonflement du projet KÇ (tableau 5.2 et 
figure 5.14). 


Tableau 5.2 Caractéristiques du projet de traitement étudié. 


Volumeinal | Pression de gonfiement Tauxde | Volume finalde résine 
dela résine | caractritique dela rsine pour | gonfiement du | à la suite du processus 
iniectée, Vatmr] | “le traitement étudié P«IMPa) |" rofet Re [1 | d'expansion, (ml 
0075358 vaoins 
(75 lives) His F7 (401 res) 


À partir des résultats obtenus, une analyse a été mise en œuvre, à l’aide du logiciel Plaxis 
2D, afin de calculer les variations de l'état de contraintes et de déformations induites par 
un traitement en colonne. Toutefois, il convient de souligner qu'il n'a pas été possible de 
déterminer l'état de déformations final, dans la mesure où la pression limite de confine- 
ment du sol a été atteinte avant que ne soit atteint l'état d'équilibre avec la pression de 
gonflement de la résine, Par conséquent, la seule information que l'on peut déduire de 
l'analyse numérique est celle relative à la variation de l'état de contraintes dans des condi- 
tions limites. 
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Figure 5.14 Définition de L pression de genflement caractéristique de la résine et du rapport de gonflement 
pour le aitement objet de l'étude, 


Les figures 5.15 à 5.19 présentent les résultats obtenus avec l'analyse numérique. 


Figure 5.15 Schème typologique des déplacements totaux du sal (efarmed mesh, amalage déformé ). 
Échelle déformé. 


Modèles numériques 


Figure 5.16 Schéma typologique des déplacements totaux du sol (arrous, flèches»). 


Figure 5.17 Contraintes verticales (shadings, dégradés). 
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Figure 5.18 Contraintes horizontales (shadings, « dégradés»). 
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Figure 5.19 Points de plastfication du sol. 


Modèles numériques 


Sur les figures 5.17 et 5.18, il peut être observé que le traitement en colonne induit dans 
Le sol environnant à la fois une variation des contraintes verticales et une variation des 
contraintes horizontales, Dans l'ensemble, il est clair que ces variations affectent une zone 
d'extension maximale d'environ 1,5 à 2,0 m à partir du bulbe d'injection 

En ce qui concerne la pression d'expansion à l'intérieur de la résine, il est clair également 
dans ce cas qu'elle n'est pas répartie de façon homogène le long de la surface du bulbe 
d'injection, mais que les valeurs maximales s'observent dans la partie inférieure du bulbe 
et ont tendance à diminuer progressivement au fur et à mesure que l'on se déplace vers le 
haut. 

Enfin, en référence à l figure 5.19, il est possible d'observer qu'une grande parte de la 
formation sablonneuse se trouve dans des conditions de plastfication du simple fait du 
processus d'expansion de la résine. Cet élément de preuve confirme que le sol entourant 
le bulbe a atteint une condition limite au-delà de laquelle il n'est plus en mesure d'oppo- 
ser une résistance au processus d'expansion de la résine, cette dernière poursuivant son 
expansion jusqu'à atteindre l'équilibre avec La pression limite. 


2.3. Effets induits dans le sol par des injections rapprochées 


Le troisième exemple de traitement des sols par injection solide (compaction grouting) à 
pour objectif de simuler ls effets induits par des injections rapprochées réalisées dans un 
milieu granulaîre homogène isotrope en l'absence de nappe (c'est-à-dire sable meuble ou 
légèrement compact). 

Le comportement mécanique de la formation sablonneuse et du bulbe d'injection a été 
mis en œuvre en appliquant des lois de comportement identiques à celles adoptées dans 
les exemples illustrés dans les paragraphes précédents (tableau 5.1) 

Le modèle géomérrique adopté dans l'analyse a été défini dans des conditions de défor- 
marions planes (plane strain), de sorte que la section du projet doit être imaginée comme 
étendue à l'infini dans le sens longitudinal, Cetre schématisation ne permet pas de simu- 
ler avec précision les effets liés au processus d'expansion de la résine injectée localement, 
au niveau des deux points d'injection rapprochés, assimilables à l'expansion de deux cavi- 
és de forme sphérique, mais peut être représentative d'une section transversale d’une 
intervention plus étendue dans laquelle des injections alignées le long de deux lignes 
parallèles ont été effectuées (figure 5.20). 

Dans l'analyse, on a supposé que les points d'injection sont situés à une profondeur de 
3m à partir du niveau de la surface naturelle du sol er que les bulbes d'injection peuvent 
être assimilés à des sphères ayant un rayon égal à 0,2 m. 


La figure 5.21 est une illustration de la géométrie de la section de référence, 
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Figure 5.20 Vue en plan de l'intervention de traitement. 


Figure 5,21 Modèle géométrique de référence pour l'analyse. 


Modèles numériques 


De la même manière que dans les exemples précédents, on a procédé en premier lieu à la 
définition de la courbe de réaction du sol au processus d'expansion de la résine polyuré- 
thane, Cette analyse a été cffecruée en référence à l'une des deux injections effectuées, 
parce que la symétrie géométrique et la charge du modèle permettent de supposer que 
l'état de contraintes du matériau au niveau des parois de chaque bulbe d'injection est sen- 
siblement équivalent. Le résultat obtenu est représenté sur la figure 5.22. 


Rapport de gonflement, Va /V, 
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_— Réaction du sol au processus é'expansion 


Figure 5.22 Courbe de réaction du sol au processus d'expansion de la résine polyuréthane, 
caleulée par le logiciel Plaxis 20, 


Comme on peut l'observer, il suffit d'une expansion réduite de la résine pour que soit 
atteinte la pression critique de confinement dans le sol environnant les bulbes d'injection, 
égale à environ 0,133 MPa. Ce scénario peut être corrélé à un phénomène diffus de plas- 
tification du matériau, en raison de la proximité en surface des points d'injection, qui pro- 
voque une chute rapide de la résistance du sol, qui n'est plus en mesure de s'opposer au 
processus d'expansion de la résine polyuréthane, Par conséquent, une fois la pression eri- 
tique arreinte, le Auide injecté conrinue à se dilater à pression constante jusqu'à artcindre 


les conditions d'équilibre dictées par la courbe rhéologique (figure 5.23). 


Le tableau 5.3 présence de façon synthétique les principales caractéristiques du traite- 


ment effectué. 


‘Tableau 5.3 Caractéristiques du projet de traitement étudié. 


Volume ni de | Pression de gonflement | Taux de gonflement | Volume final de résine 
la résine injecté, | caracérctique delarésne | duprojet pour | à la suite du processus 
pourchaque | pourletalementén, | chmue injection | … d'expansion pour 
injection, V1 PelMPa] Ref chaque injections, 
0010053 : ones 
(@ 10e) LS se (= H7 live) 
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Figure 5.23 Définition de La pression de gonflement caractéristique de La résine et du rapport de gonflement 
pour le traitement objet de l'étude. 


À partir des observations effectuées précédemment, une analyse numérique a été réalisée 
afin d'évaluer l'état de contraintes-déformations du sol au niveau de la zone traitée. 
Cependant, il convient de préciser que l'obtention de la pression limite de confinement 
du matériau naturel entourant les bulbes d'injection n'a pas permis de définir avec préci- 
sion les déplacements liés à l'expansion complète de la résine, parce que le logiciel n'est 
pas en mesure de simuler des expansions volumétriques majeures comparativement à 
celles ayant conduit à l'obtention de La pression limite. 


Les figures 5.24 à 5.28 (pages 54 à 82) montrent les principaux résultats obtenus 
avec les analyses effectuées à l'aide du logiciel Plaxis 2D 

Sur les figures 5.24, 5.25 et 5.28, on peut observer que dans la zone où sont effectuées 
Les injections, les conditions de plastification sont atteintes dans un grand nombre de 
points, dictées par l'obtention de La pression limite de confinement. Ce phénomène, qui 
est interprété par le logiciel de la même façon qu'une cinématique de rupture en cas de 
dépassement de la capacité portante du sol doit, en réalité, être considéré comme un indi- 
cateur du fait que la pression limite de confinement de la formation sablonneuse est 
atrcinte, au-delà de laquelle la résine continue à gonfler librement jusqu'à atteindre l'équi- 
libre associé à la courbe rhéologique. 


Modèles numériques 


Figure 5.24 Schéma typologique des déplacements totaux du sol (deformed mesh, «maillage déformé). 
Échelle déformée. 


Figure 5.25 Schéma typologique des déplacements totaux du sol (arrous, «lèches»). 
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Figure 5.26 Contraintes verticales (shadings, dégradés) 


Figure 5.27 Contraintes horizontales (shadings, dégradés), 


Modèles numériques 


Figure 5.28 Points de plastfication du sol. 


Surles figures 5,26 et 5,27 sont représentées les variations de l'état de contraintes induites 
par le processus d'expansion de la résine polyuréthane dans le sol environnant, qui, 
lorsque la valeur limite de la pression de confinement est atteinte, restent inchangées bien 
que le fluide injecté continue à gonfler jusqu'à atteindre l'équilibre associé à la courbe 
rhéologique. En parteulier, on peut observer que la formation sablonneuse est vouchée 
par l'effet du traitement jusqu'à une distance d'environ 2 à 3 m par rapport aux points 
d'injection. On obtient la valeur maximale de la variation des contraintes au niveau du 
matériau situé entre Les deux bulbes, lequel subit une compression significative pendant Le 
processus d'expansion de la résine. 

En conclusion, on peur affirmer que le traitement réalisé a un fer significatif sur le maté- 
riau traité, dans une zone autour du point d'injection dont l'étendue maximale est de 
l'ordre de 2 à 3 m. Cet effet provoque une augmentation des contraintes de confinement 
initiales et, par conséquent, une amélioration des caractéristiques mécaniques du fait de 
la densification rapide du sol, 


3. Intervention de récupération 
des affaissements 


Les interventions par injections de compensation sont destinées à compenser partielle- 
ment les affissements excessifs susceptibles d'apparaître sous des fondations superf- 
cielles, voire même semi-profondes. 


Intervention de récupération des affaissements 


Cet objectif est atteinr grâce à l'expansion de la résine polyuréthane à l'intérieur de frac- 
tures préexistantes ou qui se sont formées à la suite du processus d'injection à proximité 
des structures faisant l’objet d'un traitement. 
Les approfondissements théoriques rapportés dans Les chapitres précédents à propos de la 
théorie de la fracturation permettent de supposer que les interventions faisant appel à des 
injections de compensation ont une plus grande efficacité dans les sols saturés et caracté- 
risés par de faibles valeurs de conductivité hydraulique (par ex.: argiles, limons). 
Des résultats optimaux sont obtenus lorsque le matériau au-dessous de la fondation est 
fracturé dans la direction horizontale (c'est-à-dire parallèlement au niveau de la surface 
naturelle du sol), étant donné que, dans ces conditions, le gonflement de la résine déter- 
mine Le soulèvement maximal de la structure 
Les sols de fondation dans lesquels, en règle générale, il est nécessaire de procéder à des 
interventions par injections de compensation sont des sols normalement consolidés ou 
peu consolidés (OCR & 1) et caractérisés par un coefficient de pression des sols au repos 
Ki = 9 0,0) inférieur à 1. Dans ces conditions, l'injection de la résine dérermine l’ou- 
verture d’une fracture dans la direction des contraintes principales minimales, ou dans la 
direction sub-horizontale, avec en conséquence la propagation de la fracture dans la direc- 
ion sub-verticale, Le gonflement de la résine provoque une augmentation des contraintes 
horizontales, jusqu'à provoquer une rotation des contraintes principales. Une fois atteint 
un état de contraintes dans lequel la contrainte verticale est la contrainte principale 
mineure, l'exécution ultérieure d'injecrions provoque la formation de fractures dans le 
plan horizontal. 

Par conséquent, le schéma d'exécution du traîrement prévoit généralement les phases sui- 

vantes: 

1. une première série d'injections de résine polyuréthane dans le sol au-dessous de la fon- 
darion. Ces injections aboutissent généralement à la fracturation du sol dans la direc- 
tion verticale; l'ouverture de ces fractures est à l'origine d'une inversion des contraintes 
principales aux alentours de la zone traitée; 

2. une seconde série d'injections de résine polyuréthane, avec formation de fractures horizon- 
tales dont l'ouverture aboutit à un soulèvement significatif des fondations superficielle. 

À la lumière de ce qui précède, une analyse a été mise en œuvre, à l'aide du logiciel Plaxis 
2D, afin de simuler les mécanismes d'ouverture des fractures dans Le sous-sol. 
Le modèle géométrique adopté dans l'analyse a été défini dans des conditions de défor- 
mations planes (plane sain), de sorte que la section du projet doit être imaginée comme 
éendue à l'infini dans le sens longitudinal. Cette simplification ne permet pas de simuler 
avec précision la géométrie effective des fractures prenant naissance dans Le sous-sol, les- 
quelles ont une érendue finie dans la direction longitudinale, mais peut représenter la sec- 
tion transversale d’une intervention dans laquelle des injections alignées le long d'une 
parallèle à ka direction longitudinale ont été effectuées. En outre, l'analyse en déformation 
plane est de toute façon représentative des variations de l'état de contraintes et de défor- 
mations qui se produisent aux alentours d’une seule fracture. 

Le matériau utilisé dans l'analyse est une argile saturée normalement consolidée sous une 

nappe, dont la loi de comportement est linéaire élastique avec critère de rupture de Mohr- 

Coulomb. Les paramètres mécaniques caractéristiques de ce matériau sont indiqués dans 

le tableau 5.4 (page suivante). 


Modèles numériques 


Tableau 5.4 Paramètres géotechniques atwibués dans l'analyse. 


Paramètres de déformabité Paramètres de résistance 
De module | coefcent | cohésion | AE | Ange de 
de Young £ | de Poisson, v | effective, c | %létistanee Au | dilatance, 1, 
IMPa] G] LOS Boum 
Aile saurée normale: 
ment consolidée L Eos 9 É 8 


La figure 5.29 présente la géométrie de la section de référence, dont la taille a été choisie 
de manière à évirer l'apparition d'effets de bord (largeur égale à 12 m, hauteur égale à 
9 m). Le point d'injection a été localisé à une profondeur de 2 m par rapport au niveau 
de la surface naturelle du sol. 


e 
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Figure 5.29 Modéle géométrique adopté dans l'anahse. 


Il a été attribué au matériau un comportement de type non drainé, mais l'analyse a été 
effectuée en termes de contraintes effectives, étant donné que le logiciel est en mesure de 
caleuler Les surpressions incerstitielles générées par les charges et les déplacements imposés. 
On a supposé que la fracture était caractérisée par une longueur fixe de 2 m, dans laquelle 
10 kg de résine Uretek sont injectés. Ainsi qu'il a été décrit dans les chapitres précédents, 
il est impossible de prédire en théorie le processus de propagation d’une fracture dans Le 
sol, car il existe des incertitudes insurmontables concernant le comportement effectif du 
matériau naturel. C'est pourquoi il fallu se référer à l'expérience des essais sur Le terrain, 
au cours desquels des excavations ont été effectuées pour évaluer la longueur des fractures 
ocasionnées par la résine. 


Intervention de récupération des affaissements 


Étant donné que le coefficient de pression au repos du matériau utilisé dans le modèle est 


inférieur à L (K, = 0,5), on a procédé en simulan les phases suivantes: 

— Phase 1: expansion de la résine à l'intérieur d'une fracture verticale 

= Phase 
expansion de la résine à l'intérieur d'une fracture horizontale 


— Phase 3: consolidation avec dissipation des surpressions interstitielles. 


après inversion des contraintes principales autour de la fracture verticale, 


Dans la phase 1, on a simulé des ouvertures croissantes de la fracture verticale, en calcu- 
lan les pressions de gonflement correspondantes. À partir de ces résultats, il a été possible 
de déterminer la pression et Le rapport de gonflement à l'équilibre de la résine, en réfé- 


rence à la courbe rhéologique du fluide injecté (tableau 5.5, figure 5.30). 


‘Tableau 55 Caractéristiques du projet de traitement étudié. 


Volume initial Pression de gonflement Taux de | Volume final de la résine 
dela résine | caractéristique de la résine pour | gonflement du | à la suite du processus 
injectée, Valm] | _le traitement étudié, P<[MPa] | _ projet, Rc F1 | _ d'expansion, V, [m°] 


uoussas 0159346 
vo: ve 
(9 litres) hé LS {2 159 litres) 
so 


250 


200 


È 


D O1 0? 03 6e 08 06 07 08 ns 1 11 12 13 té 
Fression de confinement et de gonflement, P [MPa] 


+ Réaction du sol au processus d'expansion 
—— Courbe rhéolagique de la résine polyuréthane 


Figure 5.30 Phase 1 - Détermination de la pression de gonflement caractéristique de La résine 
et du rapport de gonflement pour le traitement envisagé. 


Dans des conditions d'équilibre, le logiciel a permis de caleuler les variations de l'état de 


contraintes et de déformations induites dans le matériau environnant. 


Sur les figures 5.31 à 5.37 (pages 86 à 89), on peut voir les éléments suivants: mail- 


lage déformé, contraintes roales verticales ct horizontales, surpressions intestirclles et 


points de plastfication. 
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Figure 5.31 Phase 1 - Modèle déformé du sol (defermed mesh, «malage déformé). Échelle déformé. 


Figure 5.32 Phase 1 - Déplacements horizontaux du sol (shadings, « dégradés») 
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Figure 5.33 Phase 1 - Déplacements venicaux du sol (shadings, « dégradés») 


Figure 5.34 Phase 1 - Contraintes venicale totales dans le sol (shadings, «dégradés). 
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Figure 5.36 Phase 1 - Surpressions interstitielles (Shadings, «dégradés 1) 
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Figure 5.37 Phase 1 — Points de plastification (plastic points). 


Les résultats obtenus montrent une inversion des contraintes principales à proximité du 
point d'injection, ce qui nous à permis de passer à la simulation de la phase 2. En appli- 
quant un mode opératoire analogue à celui de la phase 1, dans ce cas également, on a 
recherché la condition d'équilibre entre la pression de gonflement de la résine contenue 
dans la fracture horizontale et la contrainte de confinement du sol environnant 
(figure 5.38 ct tableau 5.6, page suivante). 
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100 


Rapport de gonflement Vi /V, 1 


9 O1 02 09 04 05 96 07 08 0$ 1 1 12 12 14 1 
Pression de confinement et de gonflement, P [MPa] 


_—e= Réaction du sol au processus d'expansion 
— Courbe rhéologique de la résine polyuréthane 


Figure 5.38 Phase 2 - Détermination de la pression de gonflement caractéristique de La résine 
t du rapport de gonflement pour le traitement objet de l'étude. 
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Tableau 5.6 Caractéristiques du projet de traitement étudié. 


Volume ina Pression de gonflement Taux de | Volume final de la résine 
de a résine | caractéristique de la résine pour | gonflement du | à la suite du processus 
injectée, V{mv] | “le traitement étudié, PeIMPa] | “projet Re [1 | | d'expansion, [n°] 
0009346 0179345 
CE 9 tres) Le er LE 179 res) 


Sur les figures 5.39 à 5.45, on peut voir les éléments suivants: maillage déformé, 
contraintes totales verticales et horizontales, surpressions interstitelles er points de plasti- 


ficaion. 


Il convient de souligner que les déplacements indiqués dans la phase 2 ont été calculés 
après remise à zéro de ceux de la phase précédente, À l'inverse, l'état de contraintes tient 
également compte de la phase précédente. 


Figure 5.38 Phase 2 - Modèle déformé du sol (deformed mesh, «maillage déformé). Échelle déformée. 
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Figure 5,40 Phase 2 — Déplacements horizontaux du so! (shadings, « dégradés»). 


CORRE: 


Figure 5,41 Phase 2 - Déplacements verticaux du sol (shadings, « dégradés) 
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Figure 5.42 Phase 2 — Contraintes verticales totales dans le sol (shadings, « dégradés s). 


Figure 5,43 Contraintes horizontales totales dans le sol (shadings, «dégradés +). 
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Figure 5,44 Phase 2 - Surpressions interstitielles (shadings, «dégradés s). 


Figure 5,45 Phase 2 — Points de plastification (plastic points). 
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Le soulèvement maximal caleulé par rapport au niveau de la surface naturelle du sol est 
égal à 4,5 em (voir figure 5.41), soit environ 50 % de l'ouverture moyenne de la fracture 
horizontale (à savoir: 17/2 = 8,5 cm). 

Enfin, on est passé à la phase 3, au cours de laquelle on a procédé à une simulation de la 
consolidation du matériau, à la suite de laquelle il s'est produit une dissipation totale des 
surpressions interstitielles générées par les processus d'expansion des fractures. 


Les figures 5.46 à 5.50 présentent les résultats obtenus en termes de maillage déformé, 
déplacements totaux, surpressions interstitielles er points de plastification. 

Dans ce cas, les déplacements tout comme l'état de contraintes ont été calculés en tenant 
compte des résultats obtenus dans la simulation de la phase 2. 


Figure 5,46 Phase 3 - Modèle déformé du sol (deformed mesh, «maillage déformé). Échelle déformée. 
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Figure 5,48 Phase 3 — Déplacements verticaux du sol (hadings, « dégradés:). 
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Figure 5.50 Phase 3 - Points de plastification (plastic points). 


Condusion 


La figure 5.51 représente également les déplacements subis par le sol au cours du proces- 
sus de dissipation des surpressions interstitelles sans cenir compte de ceux caleulés pen- 
dant la phase précédente, de manière à mettre en évidence les variations de l'état de 
déformations associé à la consolidation. 


Figure 5,51 Phase 3 - Modèle déformé du sol avec remise à zéro ( 


falisaton) des déplacements calculés 
dans les phases précédentes (deformed mesh, «mailage déforné 


Échelle déformé. 


Il peut être observé que la phase de consolidation comporte une réduction du soulève 
ment obtenu dans des conditions non drainées (phase 2) d'ampleur contenue (l'abaisse- 
ment du niveau de la surface naturelle du sol varie de 0,45 em, au niveau de l'axe de la 
fracture, à 1,1 em, au niveau des bords du modèle). 

Les résultars obrenus avec le modèle numérique font ressortir l'importance des conditions 
de plastfication pour l'évaluation de l'état de contraintes et de déformations aux alen- 
tours de La fracture. Il en résulte que les modèles purement élastiques, généralement utili- 
sés pour l'étude de la propagation des fractures de façon analytique, conduisent à des 
résultats très éloignés du comportement réel des sols 


(| 4. _ Conclusion 


L'analyse numérique est l'outil le plus approprié pour l'analyse des effets induits par les 
injections de résine polyuréthane Urerck dans Les sols, car elle permet de tenir compte des 
effers de bord (c'est-à-dire proximité du niveau de la surface naturelle du sol au point d'in- 
jection), des conditions de plastification du sol ainsi que de l'hétérogénéité er l'anisotro- 
pie susceptibles de caractériser Le site étudié. 


Les analyses exposées dans ce chapitre sont des exemples de Ia façon dont peuvent être 
étudiés divers problèmes d'application où les effets produits par les inj 
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sont très différents les uns des autres (injections solides [ermpaction grouting] ou injections 
de compensation [compensation groutingl). Dans tous les cas, les résultats obtenus ne 
peuvent pas être estimés avec des modèles analytiques simplifiés, fondés sur des hypo- 
thèses simplificatrices, comme celles d'un milieu infiniment écendu ou d'un comporte- 
ment purement élastique. Il convient, toutefois, de noter que le degré de précision très 
pointu de l'analyse numérique s'accompagne obligatoirement d’une caractérisation géo- 
technique extrêmement déraillée. Par ailleurs, il ne fait aucun doute que l'analyse numé- 
rique est bien plus complexe que l’utilisation de solutions analytiques. Il incombe, par 
conséquent, au concepteur d'évaluer le niveau de précision de l'analyse qu'il souhaite 
effectuer et de toujours se ménager des marges de sécurité suffisantes pour tenir compte 
des simplifications adoptées 


Risque de liquéfaction 


1: _ Introduction 


Le succès des interventions de consolidation par le biais d'injections de résine Uretek doit 
être vérifié a poseriori de manière adéquate, grâce à une comparaison entre les paramètres 
obtenus lors des essais effectués avant et après l'intervention de remise en étar. Ces essais 
peuvent inclure des essais en laboratoire ou, plus fréquemment, des essais in situ, tels que 
par exemple des essais de pénétration statique et dynamique 

À partir de ces essais, il est possible de procéder à une vérification de la sécurité du sol de 
fondation traité avec des résines expansives, dans l'optique de confirmer son adéquation 
par rapport aux charges transmises par les fondations dans des conditions sismiques. En 
particulier, il est possible d'évaluer Le risque de liquéfaction du sol sous l'action d'un 
séisme, conformément aux normes sismiques européennes (norme EN 1998, 2005, 
Eurocode 8). 


Il convient de noter que le phénomène de liquéfaction concerne uniquement certains 
types de dépôts ou formations géologiques. Nous n'étudions ici que les sols à grains gros- 
siers (sables et graviers), meubles, saturés. En effet, ces sols, même lorsqu'ils sont secs, 
peuvent subir un compactage lors d'un événement sismique lié à l'action des cycles de 
charges-décharges. L'ampleur de ce phénomène est liée au degré de compacrage initial du 
sol considéré, aux contraintes de confinement (profondeur), à la magnitude de la secousse 
sismique er à sa durée (Lai er al, 2009). 

Dans le cadre des études d'ingénierie, il est nécessaire d'évaluer la sensibilité à la liquéfac- 
tion d'un dépôr de sol er son risque cfecrif de survenue, par rapport à un mouvement sis- 
mique attendu, érant donné que les conséquences de ce phénomène, en termes 
d'affaissements ou d'éventuels phénomènes d'instabilité, sont particulièrement graves. 


[2._ Évaluation du risque de liquéfaction 


Dans un sol à grains grossiers, peu compact, saturé, le phénomène de liquéfaction induit 
par une secousse sismique implique une perte partielle ou totale de la résistance au cisail 
lement et de la rigidité, à la suite de l'augmentation de la pression intersttille qui peut 
s'accompagner d'une forte réduction, voire même de l'annulation, des efforts effectifs du 
sol, ce qui entraîne un effondrement ou des déformations permanentes significatives. 

Étant donné que l'application d'un effort de cisaillement entraîne, dans des conditions 
non drainées, le développement de pressions interstitielles positives dans Les sols où il se 
manifeste, ct, dans des conditions drainées, une diminution de volume, la sensibilité à la 
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liquéfaction peut être mise en rapport avec l'état initial du sol, identifié par les valeurs de 
l'indice des vides () et par la tension moyenne efficace (p'), sur la base de la théorie de 
l'état critique. La courbe d'état critique, montrée par la figure 6.1, sépare Les états initiaux 
pour lesquels le sol, s'il est soumis à des charges déviatrices en conditions drainées, est 
caractérisé par un comportement contractant où dilatant, En particulier, les états initiaux 
représentés par les points qui se trouvent au-dessus de la ligne d'état critique sont ceux 
pour lesquels se manifestent des diminutions de volume en conditions drainées et se déve- 
loppent des pressions intersttielles positives en conditions non drainées. Par conséquent, 
ces états sont ceux dans lesquels le sol est sensible à la liquéfaction quand il cst soumis à 
des efforts de cisaillement cycliques. 


serabe 


Lime état ent 


CT 


Figure 6.1 Sensibilité à la liquéfaction en relation avec les paramètres d'état initiaux, e = indice des vides, 
= tension moyenne efficace (modifiée de Kramer, 1996) 


Les essais en laboratoire ont montré que le phénomène de liquéfaction était susceptible de 
se produire pour les matériaux sableux saturés, uniquement sous certaines conditions 
intrinsèques et de charge (Ishihara, 1993; Kramer, 1996). Outre les paramètres d'état ini- 
taux déjà cités, à des fins d'évaluation de la sensibilité à la liquéfaction, il est d'autres élé- 
ments particulièrement importants, tels que l'historique temporel des sollicitations 
sismiques auxquelles est soumis un élément de sol, l'historique de l'état de contraintes, sa 
teneur en fine et son indice de plastification. Les dépôts sableux propres, meubles, récents, 
à granulomérrie uniforme avec des particules arrondies, er peu profonds sont les plus sen- 
sibles à la liquéfaction (Laï er al, 2009) 

L'évaluation de la sensibilité à la liquéfaction des sols granulaires est effectuée par des 
méthodes empiriques fondées sur des essais géotechniques é situ, étant donné qu'il est 
extrémement difficile de procéder à un échantillonnage non perturbé de ces rypes de 
matériaux er qu'en outre, il est impossible de se fier à des échantillons reconstitués en 
laboratoire, qui ne reproduisent pas correctement la structure du sol én si ni même l'in- 
fluence de l'historique de l'état de contraintes. 

Les normes en vigueur précisent que l'évaluation du coefficient de sécurité à la liquéfac- 
ion peut être réalisée avec des procédures historico-empiriques, dont les plus courantes 
sont certainement celles relatives à la contrainte cyclique induire par un tremblement de 
terre, selon l'approche proposée initialement par Seed & Idriss (1971), fondée sur 'inter- 
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prétation des résulrats des essais de pénétration dynamique (SPT, Srendard Penerration 

Jia, réalisés sur des sites ayant déjà subi dans le passé des phénomènes de liquéfacrion. 

L'historique temporel des accélérations induites par l'événement sismique est traité 

comme une excitation cyclique d'amplitude constante, er le risque de liquéfaction est éva- 

lué en comparant les sollicitations induites avec un paramètre de résistance cyclique, 

brenu de manière empirique par des études de cas de liquéfaction survenus dans le passé 

ou par des essais en laboratoire (NTC 08, article 7.11.3.4.3). 

À une profondeur donnée par rapport au niveau de la surface naturelle du sol, on estime 

les variables suivantes : 

+ lasollicitation sismique d'une couche de sol, exprimée en termes de CSR (Gele Siress 
Ratio, «taux de contrainte cyclique») ; 

+ la capacité du sol à résister à la liquéfaction, exprimée en termes de CRR (Gel Resis- 
une Ratio, «taux de résistance cyclique»), qui correspond au taux de contrainte 
cyclique auquel se produit La liquéfrction. 


On suppose qu'il est possible d'avoir une liquéfaction, dans le cas où la sollicitation sis- 
mique (CSR) est supérieure à la résistance (CRR). Le coefficient de sécurité à la liquéfac- 
tion, FSL, est donc défini comme suit 

CRR 

GSR 
Du point de vue des performances, il est impératif que la pertinence du facteur de sécu- 
rité soit évaluée en fonction du niveau de connaissance disponible à propos du site et des 
risques associés pour les ouvrages qui y sont érigés (voir article 7.11.3.4.3 de la norme 


NTC 08). 


L'Eurocode 8 (norme EN 1998-5, 2005), en revanche, quantifie les marges de sécurité et 
arme qu'un sol doit être considéré comme sensible à l liguéfacrion lorsque la contrainte 
de cisaillement générée par un tremblement de terre à une profondeur donnée est supé- 
rieure à 80 % de la contrainte critique ayant provoqué la liquéfaction lors de séismes pas- 
sés à la même profondeur. Un niveau de contrainte de cisaillement égal à 80 % de la 
valeur critique correspond à un facteur de sécurité de 1,25 (Lai er al, 2009). 


FSI 


611 


2.1. Évaluation du taux de contrainte cyclique (CSR) 


L'état de contrainte induit par un séisme dans le sol suit une loi de variation dans Le temps 
qui est irrégulière, difficile à reproduire dans les essais en laboratoire: par convention, cet 
historique est représenté par une série équivalente de cycles de contraintes sinusoïdaux 
(Lai et al., 2009). En supposant des cycles d'une amplitude égale à 65 % de la contrainte 
de cisaillement maximale, Seed er al. (1975) ont déduit une relation entre le nombre de 
cycles équivalents qui produisent une augmentation de la pression interstitiele égale à 
celle de l'historique irrégulier associé à des tremblements de terre enregistrés et à la magni- 
tude correspondant à ces enregistrements (figure 6.2, page suivante). 


101 


Risque de liquéfaction 


Moyenne + 1 
dar 


Moyenne 1 
écarte 


Aombre de cycles équivalents pour 0.65 v,. 
3 __# 


5 7 ë D 
Megniuce 


Figure 6.2 Relation entre le nombre de cycles équivalents etla magnitude (d'après Lai et al, 2009, 
et selon Seed et a, 1975). 


En prenant comme référence un séisme de magnitude 7,5, la sollicitation induite par 
l'historique temporel correspondant peut être assimilée à une sollicitation cyclique avec 
un nombre de cycles équivalents égal à 15 et une amplitude +, égale à: 

2,206, 1621 
où ts représente la valeur maximale de la contrainte de cisaillement, qui peut être obte- 
nue à partir d'une étude de la réponse sismique locale du dépôt de sol, ou, en variante, qui 
peut être estimée à l'aide de certaines considérations concernant l'équilibre d'un élément 
de sol soumis à une accélération horizontale, avec la formule simplifiée suivante 


a 
en = 006 16.3] 
& 


+ amtst l'accélération horizontale maximale à la surface générée par le tremblement de 
terre, qui peut être estimée en suivant l'approche simplifiée proposée par les normes 


NTE 08, (comme $ 4) ; 

+ gest l'accélération de la pesanteur ; 

+ 0, est la pression géostatique verticale totale ; 

+ rgest un coefficient de réduction de la contrainte, qui tient compte des phénomènes 
d'amplification sismique induits par la déformabilité du sol. La valeur moyenne de r, 
peut être calculée comme (Youd et al, 2001) : 

r,=10-0,00765:z z<9,15m [6.4] 
1 =1,174-0,0267:z 915m<z<2%m 
Le rapport de contrainte cyclique CSR est donné par l'amplitude de la contrainte cyclique 
de cisaillement normalisé par rapport à la contrainte géostatique verticale effective O9: 
z. au © 
CSR = Le 20,65. 70. (6.5| 
4 dan 1 
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2.2. Évaluation du taux de résistance cyclique (CRR) 


La méthode idéale pour évaluer Le taux de résistance cyclique (CRR) consiste à prélever et 
A soumettre à des essais en laboratoire des échantillons de sol intacts (non perturbés). Tou- 
tefois, dans la pratique courante, compte tenu des difficultés rencontrées pour obtenir des 
échantillons inticts (non perturbés) de sols à grains grossiers, on utilise des corrélations 
empiriques, fondées sur les résulrats des études géotechniques in su. 


Les essais es plus couramment cfecrués dans ce contexte sont les suivancs : 


+ essai standard de pénétration (SPT) ; 
+ essai de pénétration statique (CPT) ; 
+ essais géophysiques d'estimation de La vitesse de propagation des ondes de cisaillement 


a. 


Les corrélations empiriques pour l' 


aaluation du CRR à partir des résultats des essais de 
pénétration SPT et CPT se basent sur un grand volume de données et sur un passé riche en 
études de cas et, pour certe raison, donnent des résultats plus fiables. Étant donné que les 
procédures d'évaluation du CRR fondées sur des essais in sit font référence à un tremble- 
ment de rerre d'une magniude de 7,5, ce qui correspond à quinze cycles de charge équiva- 
lents, l'estimation doit être corrigée pour tenir compre de la magnitude du séisme attendu. 
Cette correction est généralement introduite en appliquant à l'estimation du CRR corres- 
pondant à une magnitude de 7,5 (CRR;s) un Facteur d'échelle MSF (Magnitude Scaling 
Factor) évaluable de manière prudente avec la formule (Youd er 4, 2001) suivante 


107% 


[6.6] 


où M, est la magnitude du moment pour un tremblement de terre atrendu sur le site 


d'inrérèr. 


En variante, l'Eurocode 8 (EN 1998-5, 2005) suggère de se référer aux valeurs proposées 
par Ambraseys (1988) et présentées dans le tableau 6.1, qui, toutefois, sont moins pré 


cautionneuses que les valeurs fournies par l'équation [6.6] 


(Lai er al, 2009). 


Tableau 6.1 Facteurs de correction à appliquer au CRR (norme EN 19985, 2005). 


Mu ms Cycles équivalents 
55 286 35 
60 220 4 
55 169 65 
70 130 10 
#0 067 2 


En tenant compte du facteur de correction pour la magnitude, la formule pour l'évalua- 


tion du facteur de sécurité à la liquéfaction (équation [6.1]) est modifiée comme suit : 


FSI 


CRR, 


CSR 


-MSF 


[6.71 
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Le coefficient MSF dans l'équation [6.7] peut aussi être interprété comme un dividende 
à appliquer au CSR pour obtenir une valeur équivalente à la valeur de référence de la 
magnitude de 7,5. 

La base de données sur les cas historiques de liquéfaction concerne uniquement les dépôts 
superficiels (d'une profondeur inférieure à 15 m), caractérisés par de faibles valeurs de 
contrainte de confinement. Les résultats des essais eycliques en laboratoire montrent que 
la résistance à la liquéfaction augmente au fur et à mesure de l'augmentation de la 
contrainte de confinement, Toutefois, cette croissance n'est pas linéaire er il s'avère donc 
nécessaire d'introduire un facteur de correction , pour des valeurs de contraintes géos- 
tatiques élevées (Youd et al., 2001) 


K, (&) CET (681 
di Sie Spa 

où 6, est la contrainte géostatique verticale effective, p, la pression atmosphérique 
(exprimée dans la même unité de mesure que 0/4) et f'un paramètre en fonction des 
conditions du site, comme la densité relative, l'historique des contraintes et l’âge du dépôt 
(tableau 6.2). Avec l'introduction du facteur de correction K,, le coefficient de sécurité à 
la liquéfaction devient : 


CRR _ CRR 
CSR  CSR 


-MSF'K, 16.9) 
Pour des profondeurs supérieures à 15 m par rapport au niveau de la surface naturelle du 


sol (pour lesquelles il existe peu de données historiques, l'extrapolation souffre de la pré- 
sence de nombreuses incertitudes. 


Tableau 62 Valeurs du paramètre f'en fonction de la variation de la densité relative (24) (Youd et al, 2001) 


F Da Le] 
08 <40% 

0.8 à 0,005 (D,-40) 40% < D, < 80% 
mm Di> 80% 


2.2.1. Évaluation du CRR dans des essais SPT 


Les méthodes pour l'évaluation du CRR à partir des résultats de l'essai de pénétration 
(SPT) sont décrites de manière détaillée par Lai er al. (2009). On trouvera ci-après les pas- 
sages fondamentaux qui portent à cette détermination. Les relations font référence à 
(i)gos c'est-à-dire à la valeur du nombre de coups par pied Vsrrnormalisée par rapport 
à La pression de confinement et rapportée à un rendement énergétique des équipements 
de l'essai égal à 60 %, qui peut être estimée en utilisant la formule proposée par Youd #t 


al. (2001): 
IN) 


sr Ca Cr"Cn° CC 16.10] 
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ombre de coups mesurés dans l'essai SPT pour une pénétration de 30 em; 
+ Gy= facteur de correction pour la profondeur d'exécution de l'essai; 


+ G= facteur de correction du taux d'énergie de la masse; 

+ G= facteur de correction du diamètre du trou de forage; 

facteur de correction de la longueur du train de tiges; 

facteur de correction de la méthode d'échantillonnage (carottier). 


Les valeurs supposées des différents facteurs de correction pour les cas les plus courants 
sont présentées dans le tableau 6.3. Par ailleurs, signalons que, pour des profondeurs 
inférieures À 3 m, la valeur mesurée de la résistance pénétrométrique pr devrait être 
réduite de 25 % (EN 1998-5, 2005). 

L'utilisation de la résistance pénétrométrique Ar est susceptible d'entrainer des erreurs 
dans l'évaluation de la vulnérabilité à la liquéfaction, en raison de sa variabilité intrin- 
sèque, de sa sensibiliré aux procédures er des incertitudes liées à la réalisation de l'essai. 


La corrélation pour l'estimation du CRR Sobrient par une représentation dans un gra- 
phique des valeurs CSR et de (M)sg associées à des événements sismiques du passé 
(figure 6.3, page suivante). Pour normaliser les données issues de plusieurs études de cas, 
les valeurs du CSR ont été ramenées à une valeur de référence de la magnitude 7,5. La 
figure 6.3 permet une évaluation directe du risque de liquéfaction en identifiant un point 
dont les coordonnées sont données par la valeur de (Ni mesurée expérimentalement er 
le CSR attendu (après correction avec le facreur MSF pour la magnitude attendue). La 
ligne de séparation entre les cas où la liquéfaction a eu lieu ou non, représente la condi- 
ion limite pour laquelle CRR = CSR er peut donc être utilisée pour estimer le CRR sur 
la base de la valeur de (No mesurée in situ. 


Tableau 6.3 Facteurs de correction de la résistance SPT 
d'après Lai et al, 2009; modifié de Youd et al, 2001) 


Facteur Caractéristique de l'équipement Conection 
Psion lasatique Cy : V7 
d4<G<17 
Z Marteau de sécurité 
Rapport Cr Marteau à anneaux 
65-115 mm 
Diamètre du sondage Gy 150 mm 
200 mm 
<3m 075 
34m 08 
Longueur de tige Ca 46m 085 
6-10m 095 
10-30m 10 
Échanéllonneur standard 10 


Méthode d'échantillonnage CS 


Échancllonnage sans protection 1113 
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Puisque la liquéfaction d'un sol est également influencée par la présence de particules 
fines (c'est-à-dire ayant un diamètre inférieur à 0,075 mm ou passant au tamis ASTM 
200), les corrélations indiquées sur la figure 6.3 ont été différenciées pour diverses valeurs 
de teneur en fines. Youd er al. (2001) ont proposé une relation pour dérerminer une 
valeur équivalente de nombre de coups par pied, (V,)sges à partir des valeurs de (V)59 
mesurées dans le sable propre 


Aoees = A+B- (NV, a (6.1 
où & er ff sont des coefficients à déterminer avec les relations indiquées dans le tableau 6.4, 


dans lequel la teneur en fines FC correspond à la fraction dont la taille de particules est 
inférieure à 0,075 mm (passant au tamis ASTM 200). 


‘Tableau 6.4 Facteurs de correction pour l'évaluation de la valeur équivalente (Ni Jancs 
rapportée au sable propre (de Lai et al, 2008, d'après les modifications de Voud et al, 2001). 


Teneur en fines FC a 0 
<5% Ü 10 
ss 
5%-35% a+ [#© 
{1000 
235% s 12 


Smile UE 


Mc 


CsAoucæR 


Figure 6.3 Relations entre le CRR et les résultats des essais SPT (de Lai et al, 2009, d'après les modifications 
de Youd et al, 2001). 


Après avoir calculé la valeur de ()socss il est possible de déterminer le CRR en appli- 
quant la formule suivante (Youd er al, 2001): 


1e, 50 1 


CRR,; = = Du 
34- (Nues 135 [0 {Macs +457 200 


16.12] 


Évaluation du risque de liquéfaction 


Cette équation est uniquement valable pour (Ni)scs < 30, étant donné que les sols gra 
nulaires très denses, pour lesquels (Vçgcs > 30, ne sont pas soumis au phénomène de la 
liquéfaction, On trouvera ci-dessous un exemple de vérification du risque de liquéfaction 
à partir des résultats d’un essai SPT, selon la méthode simplifiée de Youd er al. (2001). On 
suppose un séisme théorique de magnitude M = 7,5 et une accélération horizontale 
maximale 4 = 0,16 g. Les valeurs Mr de l'essai pénétrométrique dynamique, dont 
l'efficacité énergétique est de 72 %, sont présentées dans le tableau 6.5. La nappe phréa- 
tique se trouve à 1,50 m du niveau de la surface naturelle du sol, comme représenté sur la 
figure 6.4; le poids volumique du matériau (sable propre, avec un pourcentage de fines 
inférieur à 5 %) est supposé égal à 19 KN/m°. 


Figure 6. Stratigraphie de référence. 


Tableau 6,5 Résultats de l'essai pénétrométrique dynamique. 


] rs (eoups/pieds) am) Mer (eoups/pleds) 
12 8 ne u 
22 5 m2 m 
32 5 n2 E 
42 Ü m2 E 
52 7 152 5 
62 8 162 2 
72 8 2 6 
32 55 182 7 
92 15 12 5 
102 10 202 » 


On commence, tout d'abord, par déterminer la sollicitation sismique, CSR (ég. (6.51), en 
fonction des contraintes géosratiques totales er effectives er de l'accélération maximale au 
sol (qui tient compte d'éventuels phénomènes d'amplification sismique). 

Ensuite, on calcule la capacité de résistance à la liquéfaction du sol, dans le cas d'un séisme 
de magnirude 7,5, CRR, ;, en fonction de la valeur de Mipr normalisée, évaluée en urilisant 
un facteur de correction Cj pour la profondeur et un facteur de correction G; pour Le taux 
d'énergie de la masse, conformément aux indications du tableau 6.6 (page suivante). En 
outre, puisqu'il s'agit d'un sable propre, les coefficients liés à la reneur en fines sont supposés 
avoir Les valeurs suivantes: & = 0 et B = 1, où les valeurs de (M)ug et de (i)époscoïncident. 
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Les résultats de la vérification du risque de liquéfaction pour le sice examiné sont présen- 
tés dans Le tableau 6.6, ainsi que dans les figures 6.5 et 6.6. 


Tableau 6.6 Vérification du risque de liquéfaction pour le site examiné. 


Profondeur ( 


2 


25 


Figure 65 Évaluation du FSL (- par rapport à la profondeur, z (m). 
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Contrainte de cisaillement (kPa) 
ge 10e m0 mo #0 50 emo mo mo wo 100 


Profondeur (m) 


—Sollicitation (kPa} Résistance (kPa) 


Figure 6.6 Comparaison entre la sliaton cyclique de ésallement, 1, induite par un séisme équivalent 
de 15 cycles etla résistance clique avec un nombre de cycles équivalents égal à 15 et une amplitude +, 
la résistance cyclique au csailement, ra. 


2.3. Indice du potentiel de liquéfaction 


La liquéfaction peut se produire localement, sans impliquer nécessairement l'effondre- 
ment ou une perte de fonctionnalité des structures en interaction avec le sol (Lai ef al. 
2009). Pour cette raîson, il st indispensable d'associer à l'évaluation ponctuelle de la sen- 
sibilité à la liquéfaction à diverses profondeurs, réalisée par exemple avec la méthode 
décrite au paragraphe « Évaluation du taux de résistance cyclique (CRR) », une évaluation 
globale de l'extension du phénomène et de ses conséquences. Une évaluation approxima- 
tive peut être effectuée, par exemple, en utilisant l'indice du potentiel de liquéfaction LPT 
(Liquefaction Potential Indes) (lwasaki et al, 1978), un paramètre intégral qui tient 
compte de l'épaisseur de la couche liquéfiable et de sa proximité avec la surface libre, mais 
aussi de la distance par rapport à la valeur unitaire du facteur de sécurité à la liquéfaction. 
L'indice LPI est défini par la formule suivante (Iwasaki er al, 1978): 


LPI = [Fte)o(2) +de (6.13 
où z est la profondeur par rapport au niveau de la surface naturelle du sol, en mètres. F 


est une fonction du facteur de sécurité à la liquéfaction, FSL, à la profondeur z donnée 
par: 


f —FSL(z) pour FSL(e)<1 
aÛ pour FSL(z)> 1 
et w(a) est une fonction de la distance de l'élément de volume par rapport au niveau de 
la surface naturelle du sol, donnée par: 

u(z)=10-0,5-2 [6.15] 
La fonction intégrale [6.13] est conventionnellement évaluée sur une épaisseur totale de 


20 m, parce qu'il a rarement été rapporté de cas dans lesquels le phénomène de liquéfac- 
tion est survenu à de plus grandes profondeurs. Le risque d'effets superficiels de La liqué- 


16.14] 
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faction, en rermes d'affaissements et d'éventuels phénomènes d'instabilité, augmente avec 
l'augmentation de l'indice de liquéfaction potentielle, L'évaluation du risque associé à la 
liquéfaction peut se fonder sur les indications fournies par Ivasaki ee al. (1982) et résu- 
mées dans le tableau 6.7; ces indications sont essentiellement en accord avec les résultats 
de l'étude de Toprak 8 Holzer (2003). Sur la base de ces indications, des manifestations 
superficielles de la liquéfaction pour des valeurs de LPI 3 > 5 et la manifestation de phé- 
nomènes d'expansion latérale pour LPI 3 > 12 (Lai er al, 2009) peuvent avoir lieu. 


Tableau 6,7 Corélation entre l'indice du potentiel de liquéfaction LPI te potentiel de rupture. 


LP Potentiel de rupture 
<5 Faible 

5-15 Élvé 

>15 Fès élevé 


En référence à l'exemple du paragraphe «Évaluation du CRR dans des essais SPT» 
(tableau 6.6), en utilisant l'équation [6.13], il est possible de dérerminer l'indice du 
potentiel de liquéfaction LPL: dans ce cas étudié, le LPI est égal à 14,36 (voir figure 6.7) 


et est associé à un fort potentiel de rupture. 


LI 
$s ss nées mes ses 


FRS 


Figure 67 Évaluation de l'indice du potentiel de liquélaction. 


3. Exclusion de la vérification du risque 
de liquéfaction 
Les normes NTC 08 (article 7.11.3.4.2) établissent que la vérification du risque de li 


faction des sols de fondation peut être omise, en présence d'au moins une des situations 
suivantes: 


1. Les événements sismiques attendus sont d’une magnitude M inférieure à 5. 
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2. Les accélérations maximales atrendues au niveau de La surface naturelle du sol en l'ab- 
sence d'ouvrages (conditions en champ libre) sont inférieures à 0,1 g. 

3. La profondeur moyenne saisonnière de la nappe est à plus de 15 m par rapport au 
niveau de la surface naturelle du sol, pour Le niveau du sol sub-hori 
tures avec des fondations peu profondes. 


tal er les strue- 


4. Les dépôts se composent de sables propres avec une résistance pénétrométrique nor- 
malisée (MX > 30 ou guy > 180. 

5. La distribution granulométrique à l'extérieur des zones indiquées dans l figure 6.8(a) 
dans le cas des sols avec un coefficient d'uniformité U, < 3,5 et la figure 6.8(b) dans 
Le cas des sols avec un coefficient d'uniformité U, > 3,5. 


Diner déni 


mon] She] Graver 
z* 
i Léon 


Dimère émmi 


Figure 6.8 Courbe granulométique des sols sensibles à la iquéfaction: (a) coefiient d'uniformité 
inférieur à 3,5; (b) coeficient d'uniformité supérieure à 3,5 (norme NTC 08). 


Siles conditions L et 2 ci-dessus ne sont pas remplies, la norme stipule que les expertises 
géorechniques doivent être conduites de manière à permerrre au moins la vérification des 
conditions 3, 4 et 5. 

Les deux premiers critères peuvent être corrélés à des critères historiques. Le critère 3, rela- 
üf à la profondeur de la nappe phréatique, se justifie sur la base de preuves historiques 
ayant démontré que, pour des profondeurs supérieures à 15 m par rapport au niveau 
naturel du sol, les conséquences sur la surface libre sont généralement d'ampleur limitée. 
Le critère 4 se base sur les corrélations empiriques qui permettent d'exclure le risque de 
liquéfaction en cas de valeurs suffisamment élevées des indices pénétrométriques. Enfin, 
le critère 5 définit la sensibilité à la liquéfaction associée à des critères granulométriques, 
en différenciant les courbes granulométriques de référence en fonction du coelficient 
d'uniformité du matériau défini comme: 
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dans quelle Dyg er Djsont les diamèrres correspondant à un ramisar (passant), respec- 
tivement, de 60 % et de 10 % dans la courbe granulométrique cumulée (Lancellotra, 
2004). 

Si aucune des conditions susmentionnées n'est satisfaite et si le sol de fondation com- 
prend des couches érendues où denses et liches de sables meubles sous la nappe, les 
normes NTC 08 prescrivent la nécessité d'évaluer le coefficient de sécurité à la liquéfac- 
tion aux profondeurs d'intérêt, 


4. Mesures d'atténuation du risque 
de liquéfaction 


Sie sol est sensible à la liquéfaction et que les effets qui en découlent semblent de nature 
à affecter les conditions de srabilité (pentes ou ouvrages), il convient de procéder à des 
interventions de consolidation du sol et/ou de transférer la charge sur des couches de sol 
qui ne soient pas sensibles à la liquéfaction [article 7.11.3.4.1 de la norme NTC 08]. 

Le choix optimal d'amélioration des sols pour empêcher la liquéfaction dépend de divers 
facteurs, notamment: 

+ la localisation, l'étendue, la profondeur et le volume des sols concernés; 

le niveau de gravité du tremblement de terre attendu; 

+ la sypologie structurale des ouvrages susceptibles d'être couchés; 

+ les conditions d'accessibilité du site. 


Il ressort clairement de la description des mécanismes à la base du phénomène de li 
faction que, pour prévenir son apparition ou au moins limiter les dégâts qu'elle o 
sionne, il est nécessaire de compacter le terrain et/ou de créer des mécanismes de drainage 
en mesure de dissiper la surpression interstitelle plus rapidement qu'elle ne se forme. Les 
autres méthodes pour réduire Le risque de liquéfaction consistent À améliorer Les proprié- 
tés mécaniques du dépôt de sol en procédant au remplissage des vides interstities 

Les interventions de consolidation au moyen d'injections de résines polyuréthanes expan- 
sives Urerek ont pour objectif d'augmenter la densité relative, D), des sols à grains grossiers 
par ex, sables er graviers), saturés er meubles, afin de limiter leur risque de liquéfacrion. 
La densité relative, ainsi qu'il est connu, est fonction de lindice des vides selon la relation 
suivante: 


(617) 


sol en question. L'indice des vides est, en outre, corrélé à la déformation volumérique €, 
au moyen de la relation suivante: 


et em représentent respectivement l'indice des vides minimum et maximum du 


[6.18] 


[6.19] 
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[6.20] 


où ep est l'indice des vides initial. 
La déformation volumérrique est donnée, dans des conditions de déformations planes, 
par la somme de la déformation radiale £, et de la déformation tangentielle cp: 


8, 2E,+6 16.21] 


dans laquelle: 
[6.22] 


16.23 


où # est le déplacement radial par rapport au centre de la cavité. 
Le déplacement radial #, est fonction de la distance radiale, », à partir du centre de la 
cavité er de la pression d'expansion, ?, selon l'expression suivante: 


BSP<R 


mn P)= ex CI (0) «el P>Pnt<n(?P) 16.24] 


P>Bnr>n(?) 


où: 

+ (2) = rayon plastique, évaluable de la façon suivante: 
\o 

n@)= P 

BL sin (pl 


+= rayon initial de la cavité ou du bulbe de sol traité; 
+ Py= pression initiale, supposée égale à la contrainte géostatique horizontale, dans le 
cas d'une eavité cylindrique, et à la contrainte géostarique moyenne, dans le cas d'une 
cavité sphérique: 
+ P,= pression d'avant la plastficarion, évaluée par l'expression suivante: 
bare], _n 
MEN 
+ _m = coefficient géométrique étant supposé égal à 1 dans le cas d'une cavité cylin- 
drique, égal à 2 dans le cas d'une cavité sphérique: 
_ 2-0 
U I-sind 


2. 


16.26] 
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tant (6.271 
ET) 
. cohésion; 
+ G= module de cisaillement, évalué comme G 5 
+ E'= module de Young; QUES] 
+ v = cocfficient de Poisson; 
+ eg = coefficient pouvant être exprimé sous la forme: 
Re 16.28] 
N+1 2G 
+ À B, C= coefficients évalués de la manière suivante: 
(6.29 
[6.30] 
16.31] 
dans lesquels: 
[6.32] 
(6.33) 
16.34] 
16.35] 
[6.36] 
[6.37] 
1=v) (+ MN) (M+ 
GENE EECAEN] 16381 


MN 
Par conséquent, la déformation radiale €, et la déformation tangentielle € peuvent être 
réécrites sous la forme: 


(nr 


B<P<R 


er) nt) 


vs 
con | sel P>Anr<n(?) 16391 


mR<P<R 
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BSP<R 


a t2)=l— [4 (2) (#2) pe P>Bor<n(?) [640] 


R:sin(é) (2) 
7 267 
Il convient de souligner qu'en présence du sol au-dessous de la nappe, toutes les corréla- 
tions présentées font référence à des analyses en condition drainée. Par conséquent, les 
contraintes indiquées doivent être comprises comme étant efficaces, 


P>Por>n(?) 


La procédure décrite permet donc de déterminer la densité relative à la suite d'une injec- 
tion de résine polyuréthane expansive sous la forme: 


1641] 


La densité relative finale sera fonction de la distance par rapport au centre de la cavité, », 
etàla pression d'expansion, 2. Elle dépendra, en outre, des paramètres physiques et méca- 
niques du sol, notamment de l'angle de frottement, g, et de l'angle de dilatance, 1}. La 
pression de gonflement P, qui, pendant le processus d'expansion, tend à être égale à la 
pression de confinement à la profondeur d'intérêt, est calculée conformément au mode 
opératoire décrit dans le chapitre 4. 

La densité relaive D) er la valeur (M)4g sont corrélées par l'intermédiaire de l'expression 
suivante proposée par Skempton (1986): 


= 16.42] 
Par conséquent, en inversant l'équation [6.42], il est possible d'obtenir une nouvelle 
valeur de (N) suite au traitement à l'aide de résines, puis une nouvelle valeur de Mipy, 
avec laquelle on peut répéter la procédure décrite au paragraphe «Évaluation du risque de 
liquéfaction» pour estimer l'augmentation du facteur de sécurité à la liquéfaction FSL et 
la réduction de l'indice du potentiel de liquéfaction, LPI. Il convient de souligner Le fair 
que, dans certe nouvelle phase de calcul, aux valeurs de contraintes géostatiques torales et 
effectives il faut ajouter les augmentations de tensions verticales otales et effectives, res- 
pectivement AO, er AG}, qui coïncident entre elles et qui peuvent être évaluées comme 
suit: 


A9, = AG, + AG,) (6.431 


où AO, et AGgreprésentent la variation des contraintes respectivement radiales et rangen- 
telles, pour la détermination desquelles on se reportera au chapitre 4 

On suppose que l'on effectue à diverses profondeurs (de 2,2 à 7,2 m), dans le sable de 
l'exemple cité au paragraphe « Évaluation du CRR dans des essais SPT », six injections de 
résine polyuréthane Uretek, On suppose également que l'expansion de la résine se produit 
chaque fois selon une cavité cylindrique d'une hauteur = 1 m, que le rayon initial du 
bulbe est égal à r, = 0,1 m er que la hauteur totale du traitement est égale à 6 m. Les para- 
mètres physiques et mécaniques du sable en question sont les suivants 


Risque de liquéfaction 


+ poids volumique du sol, y = 19 kN/m?: 
+ module élastique de Young, F'= 10 MPa; 
+ coefficient de Poisson, = 0,25; 

+ cohésion effective, e'= © kPa; 

+ angle de frottement, g = 35°; 


+ angle de dilatance, dr 
= degré de surconsolidarion, OCR 
+ indice des vides minimum, és 
+ _ indice des vides maximum, 6x = 


À chaque profondeur d'injection, on évalue, grâce aux résultats de l'essai de pénétration 
dynamique SPT, la densité relative initiale, D, l'indice des vides initial eye la porosité 
initiale, mg. En outre, on calcule le coefficient de pression des cerres au repos, &, la pres- 
sion horizontale initiale effective, Pj, erla pression d'avant la plastficarion, P. 


En utilisant la méthode décrite au chapitre 4, on estime la pression d'équilibre totale 
P,jque l'on obtient lorsque la pression de gonflement de a résine est égale à la contrainte 
de confinement. 

Les relations [6.39], [6.40], [6.20], en posant r= 7, et P = P., permettent de déterminer 
respectivement la déformation radiale, la déformation tangentielle er la déformation volu- 
métrique qui sont obtenues au niveau du centre de l'injection en présence de la pression 
d'équilibre. Enfin, avec l'éq. [6.41], il est possible de déterminer la densité relative finale, 
Dies Pas à la suite de l'intervention ec la comparer à la valeur initiale obtenue avec 
l'ég 16.42] à partir des valeurs de (Ni). Les résultats obtenus pour les différents points 
d'injection à des profondeurs différentes, sont présentés dans le tableau 6.8. 


‘Tableau 68 Estimation de la densité relative finale au niveau du point d'injection. 


CR Do “ mo] Ft) | Pate) | 4er) 
22 5 o4r 065 | 0 | 10 | 2461 | 04 
32 5 os | 065 | oëo | 1e | 271 | 04 
#2 Ü o4 | 06 | 0% | 273 | 22 | où 
52 7 ou | 06 | os | 2665 | 3560 | 04 
62 O 08 | oc os | 305 | 3854 | 04 
72 8 ou | 06 | os | sw | ao | o4 


Par la suite, avec les résultats obtenus, il est possible d'évaluer, au moyen de l'ég. [6.42], 
les nouvelles valeurs de (My et, avec elles, de réitérer la procédure décrite au para- 
graphe « Évaluation du risque de liquéfaction», pour estimer le nouveau facteur de sécu- 
rité à la liquéfaction, FSL, et l'indice du potentiel de liquéfaction, LPI. Les résultats de 
l'analyse sont présentés dans le tableau 6.9, ainsi que sur Les figures 6.9, 6.10 et 6.11. 


Mesures d'atténuation du risque de liquéfaction 


Tableau 6:9 Vérification relative à la liguéfaction postintevention. 


FL 


Wal re [ns Se | ce 


2al io) as] 34] 6557 lossl 17 | aa 110] 0109 | 129 | 129 | 118 


32/10] 608 | 458] 7652 [08] 107 | 08 139] oué | 0120 | 144 | 103 


42] 0] 78] 528] 242 [om] ni] os] 164] ot | 0123] 167 | 102 


52/15 088 lens] au26 loso] 125 | 096! 187 | 0124 | 0136 | 204 | 1,09 


62/15 {ul 708 | 9303 [077] 134 | 095[ 209 | 0127 | ou48 | 242 [16 


22[ 15 {1368/7098 | o82 [o7s| 132 | 4] 251] 0130 | 0143] 254 | so 


Contrainte de cisaillement (kPa) 
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z 


— Sollitation (kPa) — Résistance avant l'intervention (kPa 
— Résistance après l'intervention (kPa) 


Figure 6.9 Comparaison entre la sollicitation cydique de cisaillement 1, induite par un séisme équivalent 
Ga 15 ace atla réstance ojcique avec un nombre de cycles équalons égal à 15e uni amplitude c. 
et la résistance cydique au cisaillement, ru 


FsLt) 
ps os 19 130 20 2% ag 35 40 
5 
£ 
Su 
È 
FI 
2 
# 
—FSL avant l'intervention FSL après l'intervention 
Figure 6.10 Évaluation du FL { par rapport à a profondeur, 2 (m). Comparaison des valeurs du facteur 


de sécurité à la liquéfaction avant et après l'intervention, 


Risque de liquéfaction 
ss 
î | 
Hd | l 
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= Fa ar 


—LPI avant l'intervention — LPI après l'intervention 


Figure 6.11 Évaluation de l'indice du potentiel de liquéfaction. Comparaison des valeurs de LPI avant 
et après l'intervention. 
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Introduction 


Ce chapitre présente deux exemples concepruels d'injections de résines polyuréthanes 
expansives dans des sols à grains grossiers er à grains fins. En effet, comme Le montrent les 
chapitres précédents, l'approche conceptuelle est différente en fonction de La granulomé- 
trie du milieu. 


Dans les sols à grains grossiers, par suite d'un effet de pénétration initial imputable à des 
valeurs de perméabilité élevée, les injections entraînent La formation d'un bulbe de sol 
traité, dont le gonflement détermine la compaction et la densification du sol environnant. 
En outre, le bulbe proprement dit constitue un élément de renfort, qui est en mesure 
d'augmenter la résistance au cisaillement er La rigidité du sol traité 

La mise en œuvre en condition non drainée provoque, en revanche, la formation de frac- 
tures dans Les sols à grains fins. L'ouverture des formations horizontales peut être exploi- 
tée pour compenser les affaissements complets et/ou différentiels des fondations 
superñcielles. 


2. injections de résines expansives 
dans les sols à grains grossiers 


2.1. Injections de résines fortement expansives 


On souhaite concevoir une série d'injections de résine polyuréthane expansive, de type 
Uretek Geoplus”, dans un sable moyen/grossier, afin de densifier le matériau et d'amélio- 
rer ses propriétés mécaniques. 

Les paramètres physiques et mécaniques du sable en question sonc les suivants: 

9 KN/m ; 


+ poids volumique du sol, y 
+ module élastique de Young, 


+ coefficient de Poisson, v 
+ cohésion effective, e'= 0 kPa 
+ angle de frottement, ç = 35 
+ angle de dilatance, = 0°; 
+ degré de surconsolidation, OCR = 1; 

+ indice des vides minimum, emin = 0,1; 
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+ indice des vides maximum, emax = 1; 

+ pourcentage de fines, FC < 5 %. 

11 ny a pas de nappe, er il n'y a pas non plus de surcharge (g = 0 kPa). Étant donné que 
le sol est normalement consolidé, le coefficient de pression des terres au repos, Kj, peut 
être calculé comme suit: K, = 1=sin(g), ce qui donne K, = 0,426. 

On dispose, en outre, des valeurs Ar (tablean 7.1) issues d'un essai de pénétration 
dynamique, dont le rendement énergétique est de 72 %. 


‘Tableau 7.1 Résultats de l'essai de pénétration dynamique. 


2m) Aer (eoups/pieds) 2m) in (eoups/pieds ) 
12 8 na m1 
22 5 122 4 
32 4 132 2 
42 4 142 2 
52 7 152 25 
62 n 162 2 
72 E 172 6 
82 E 182 7 
92 15 192 5 
102 1 202 E] 
1 pied = 0,3048 mètre 


On souhaite effectuer, à diverses profondeurs (de 2,2 m à 7,2 m), six injections de résine 
polyuréthane Uretek. On suppose que l'expansion de la résine se produit chaque fois 
selon une cavité cylindrique d'une hauteur A = 1 m, que le rayon initial du bulbe est égal 
à 45 = 0,1 m et que la hauteur totale du traitement est égale à 6 m. 

À partir des valeurs Mspy on évalue la densité relative du sable aux différentes profondeurs 
d'intérêt, ainsi qu'il est illustré au chapitre 6: après avoir caleulé la contrainte verticale 
totale et efficace, on détermine (M) c'est-à-dire la valeur de Niprnormalisée par rap- 
port à la pression de confinement et rapportée à un rendement énergétique des équipe- 
ments de l'essai de 60 %, en supposant que tous les coefficients de correction sont égaux 
à l'unité, sauf Cy. Par la suite, a densité relative initiale D, peut être dérivée de la valeur 
(No en utilisant l'expression [7.1] proposée par Skempron (1986): 


N, 
De [Me 7.1 
co 
Les résultats de la procédure décrite ci-dessus sont présentés dans le tableau 7.2. 
À partir de la définition de la densité relative 
[72] 


on obtient la valeur de l'indice initial des vides 4 et celle de Ia porosité initiale my, ainsi 
qu'il est indiqué dans le tableau 7.3. 


Ce. b 


injections de résines expansives dans es sols à grains grossiers 


Tableau 7.2 Calcul de a densité relative initiale D), à pari des valeurs de l'essai de pénétration dynamique. 


2m) Mars to (hPa) | 0° (kPa) Ge CA De 
22 5 A8 A8 252 mn 039 
32 4 cos cos 126 60 02 
42 4 m4 m4 10 53 030 
52 7 o88 ve 05 83 037 
62 ïn 78 78 050 ns 045 
7 E 1368 1368 084 13,1 047 


‘Tableau 7.3 Valeurs de l'indice initial des vides et de La porosité initiale. 


2m) a me 
22 05 039 
32 072 042 
42 073 042 
52 067 040 
62 060 038 
72 0,58 037 


On passe ensuite au calcul de la contrainte initiale aux alentours de la cavité avant l'ex 
pansion P, qui, dans le cas de l'expansion d’une cavité cylindrique, est égale à la contrainte 
initiale horizontale, 04, obtenue en multipliant la contrainte verticale initiale, 0,9, par le 
coefficient de pression des terres au repos, Kç; puis, plus tard, au caleul de la contrainte 
d'avane la plastfication P3, qui est égale à: 


BR =R -(+sing) [73] 


‘ableau 74 Valeurs de la contrainte initiale et de la pression d'avant la plastification. 


2m) Pa Pa) Pi (kPa) 
22 17.83 2805 
32 25,93 4080 
42 3403 5355 
52 4213 66,30 
62 5023 79,05 
72 5834 9179 


Le volume initial, V, du bulbe injecté à différentes profondeurs est toujours égal à 
3,14107 m?, tandis que le volume de résine initial, V,, est obtenu en multipliant V, par 
la porosité initiale #. Sachant que la densité de la résine liquide, p,, est égale à 1 070 kg/mÿ, 
il est possible de calculer la masse de résine initiale, A4, à injecter aux différentes 
profondeurs. 


a 
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Tableau 7.5 Valeurs du volume de résine initial et de la masse initiale de la résine. 


2m) Va) M Ge) 
22 124102 1324 
32 231102 toi 
42 133102 1422 
52 126-102 1343 
62 118102 12,59 
72 115102 12,34 


Le volume final du bulbe, Ve peut être déterminé à partir du rayon initial de la cavité, en 
fonction de la pression de confinement ?. En outre, V,, étant connu, il est possible de 
déterminer la variation de volume subie par Le bulbe en fonction de la pression ?. En sup- 
posant que cette variation de volume, AV, soie attribuable à l'expansion de volume de la 
résine, il est possible d'obtenir Le volume final de la résine V,+ 

Le volume de résine infectée étant connu, on peut, par conséquent, déterminer le taux de 
gonflement de la résine V,/V en fonction de la contrainte de confinement du sol et Le 
comparer avec la loi rhéologique en fonction de la pression de gonflement de la résine : 
l'intersection des deux courbes permet d'ebrenir la valeur d'équilibre final *, entre la 
contrainte de confinement du sol er la pression de gonflement de La résine, telle que repré- 
sentée sur la figure 7.1 


> ul 
> 1 


Er ee 
TT pression de gonfoment 


PRont rab d kabi bkpab rt ob, 


book ab 4 Kai Bkpab fe ob D! A ob I KakO bp 
bre mind gent ed UUR 


Figure 7.1 Rapport de gonflement de la résine en fonction de La pression de gonflement de la résine 
rouge) et dela contrainte de confinement du sol (igne bleue) ; évaluation de a pression 
d'équilibre P à une profondeur de 2,2 m par rapport au niveau de la surface naturelle du sol. 


Le tableau 7.6 présente les valeurs de la pression d'équilibre P* que l'on obtient aux dif 
férentes profondeurs d'injection. Pour une pression donnée, il est également possible de 
calculer les valeurs du rayon final, 4, de la cavité, du rayon plastique, , et du volume final 
de la résine, V+ 


injections de résines expansives dans es sols à grains grossiers 


Tableau 7.6 Valeurs de la pression d'équilibre, Pr, du rayon final de La cavité, a, 
du rayon plastique, b, et du volume final de La résine, Vs, 
aux diférentes profondeurs d'injection. 


20) PrPs) a tn) BG) Val) 
22 253 10 218 009 
32 DA 133 17 ous 
42 387 176 146 007 
52 4178 16 125 007 
62 2510 162 109 0062 
72 485,0 15 08 0058 


Avec ces valeurs de pression d'équilibre, il est possible d'évaluer l'évolution de la densité 
relarive, Dr, aux alentours du trou d'injection, à partir du caleul des déformations radiales 
et tangentielles, respectivement £, et Eg, ainsi qu'il a été montré dans le chapitre 6. 
Comme on peut Le voir, par exemple, sur le graphique de la figure 7.2 dans le cas d’une 
injection à 2,2 m par rapport au niveau de la surface narurelle du sol, la densité relative a 
tendance à décroitre à partir du centre de l'injection, et prend la valeur initiale à une dis- 
tance supérieure au rayon plastique, dès lors que les déformations volumériques sont 
nulles. 


DE 
si 
Eu 
à 
Lu 1 2 E] 4 


Figure 7.2 Tendance de la densité relative aux alentours du point d'injection situé à une profondeur 
de 2,2 m par rapport au niveau de la surface naturelle du sol. 


Le tableau 7.7 (page suivante)présente les estimations des valeurs de la densité relative 
finale, D, que l'on obtient, aux différentes profondeurs d'intérêt, à la pression d'équilibre, 
à une distance » = 4 par rapport au point d'injection. En outre, il convient de souligner 
qu'à certe distance, en moyenne, on peut obtenir des améliorations de la densité relative 
d'environ 10 % de la valeur initiale. 


D a 
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Tableau 7.7 Estimation des valeurs de a densité relative finale, en fonction de la pression d'équilibre, 
à une distance de l'axe d'injection égale au rayon intl de la cavité, et valeurs correspondantes de N sv 


ztm) pra) DV= a) | Dr 29/0, | Na (oups/pieds) 
22 2653 van 109 10 
32 DA 0360 on 8 
42 war 0346 7 ë 
52 A4 043 om 3 
62 4510 va 12 19 
72 4350 052 22 2 


Ces valeurs de densité relative peuvent être transformées, grâce à la corrélation empirique 
[7.1] proposée par Skempton (1986), en valeurs de (IV) 59 et, par conséquent, Npr. 
(figure 7.3) à comparer avec les résultats d'éventuels essais de pénétration dynamique 
post-traitement de consolidation afin de comparer et de valider les résultats de la méthode 
de caleul décrite dans ce document. 


Essai de pénétration dynamique SPT 
Aer In° colpipiede] 


5 CT | | | 
So : 
". TL | | | 
marre As irauan 


Figure 7.3 Comparaison entre les valeurs de SPT aux profondeurs d'injection, avant et après l'intervention. 


2.2. Injections de résines moins expansives 


Afin de mettre en évidence Les effets positif du traitement de consolidation avec des injec- 
tions de résine polyuréthane expansive, de type Uretek Geoplus”, on rapporte ensuite les 
résultats de la conception du même type d'intervention avec une résine caractérisée par une 
pression de gonflement maximale d'environ 500 KPa et par La loi théologique suivante: 


10,5 


ee [7.4] 
21 +0” +P (MPa)] 


<S 


Les différences entre la courbe rhéologique des résines Urerek Geoplus* er celle des résines 
moins expansives peuvent être appréciées sur la figure 7.4. 


injections de résines expansives dans es sols à grains grossiers 


sine rio = sine moins exparae 


Rapport de gonflement, V,/ V, 


Pression de gonfloment, P [MPa] 


Figure 7.4 Tendances des courbes théologiques respectives des résines Uretek Geoplus® rouge) et 
des résines moins expansives (bleu), dent la pression de gonflement maximale est de l'ordre de 500 kPa. 


Dans ce cas aussi le volume de résine injectée étant connu, on peut déterminer le rapport 
de gonflement de la résine V,{V,en fonction de la contrainte de confinement du sol et le 
comparer à La loi rhéologique [7.4] en fonction de la pression de gonflement de La résine: 
l'intersection entre les deux courbes permet d'obtenir la valeur d'équilibre final, 2, entre 
la contrainte de confinement du sol et la pression de gonflement de la résine, telle que 
représentée sur la figure 7.5, Comme on peut le constater, dans ce cas, la pression d'équi- 
libre P* est nettement plus faible que dans le cas d'injections de résine polyuréthane 
expansive. Le tableau 7.8 (page suivante) présente les valeurs de la pression d'équilibre P* 
que l'on obtient aux différentes profondeurs d'injection. Pour une pression donnée, il est 
également possible de caleuler les valeurs du rayon final, 4, de la cavité, du rayon plas- 
tique, b, et du volume final de la résine, Ve 


So NO 10 20 AD 30 9 20 40 50 50 
Pébood k ab d kel bkpab tr ob bp! kr ab ] Kaki bep 
sb tb mag dl gor ndlr db, / (F/RI 


Figure 7.5 Rapport de gonflement de la résine en fondation de sa pression de gonflement ligne rouge) 
et dela contrainte de confinement du sol (ligne bleue); évaluation de la pression d'équilibre P* 
à une profondeur de 2,2 m par rapport au niveau du terrain naturel. 
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Tableau 7.8 Valeurs de la pression d'équiibre, P, du rayon final de la cavité, a, du rayon plastique, b, 
et du volume final de La résine, V, aux diférentes profondeurs d'injection. 


24m) P° (kPa) a (em) 5m) Va) 
22 1547 120 104 0.026 
32 1870 n7 os 0025 
42 230 15 067 0.023 
52 2310 m3 055 0021 
62 2460 ni 048 0019 
72 2610 no 042 001 


Les valeurs de pression d'équilibre étant connues, on évalue la tendance de La densité rela- 
tive, D,, aux alentours du trou d'injection: comme on peut le voir, par exemple, sur Le gra- 
phique de la figure 7.6 dans le cas d'une injection à 2,2 m par rapport au niveau du 
terrain naturel, la densité relative a tendance à décroitre à partir du centre de l'injection, 
et prend la valeur initiale à une distance supérieure au rayon plastique, dès lors que les 
déformations volumérriques sont nulles, 


FT 


a 


se üs ï E 2 
212433 a 6some7 120 


Figure 7.6 Tendance de la densité relative aux alentours du point d'injection situé à une profondeur 
de 2,2 m par rapport au niveau du terrain naturel. 


Le tableau 7.9 présente les estimations des valeurs de la densité relative finale, D), que 
l'on obtient, à différentes profondeurs d'intérêt, à la pression d'équilibre, à une distance 
1 = ay par rapport au point d'injection. En outre, il convient de souligner qu'à certe dis- 
tance, en moyenne, l'ampleur des améliorations s'avère moindre que celles obrenues avec 
des injections de résine Uretek. 


injections de résines expansives dans es sols à grains grossiers 


Tableau 7.9 Estimation des valeurs de a densité relative finale, en fonction de la pression d'équilibre, 


à une distance de l'axe d'injection égale au rayon initial de la cavité, et valeurs correspondantes de NSPT. 
2 (m) P* kPa) D,tr= a) | D4r=a9/0,s | Nam (eoups/pieds) 

22 1547 0,408 105 8 

32 1870 033 107 6 

42 2130 0321 108 6 

52 2310 0,3% 107 10 

62 2460 0,469 105 15 

72 261,0 0,490 105 E 


Ces valeurs de densité peuvent être transformées, grâce à la corrélation empirique [7.1] 
proposée par Skempton (1986), en valeurs de(M)s et, par conséquent, MGpr. 

Par ailleurs, sur le graphique de la figure 7.7, on peut voir également la tendance des 
valeurs de Meyri comme on peut le constater, les injections de résine Uretek s'accom- 
pagnent d'une densification majeure du sol, liée à la présence de pressions d'équilibre, P*, 
de valeurs sensiblement plus élevées et, par conséquent, à des valeurs de Apr. plus élevées 
comparativement à un traitement avec l'autre type de résine, caractérisé, quant à lui, par 
une faible capacité de gonflement. 


Essai de pénétration dynamique SPT 
Nar{n° colpipiede] 
9 5 10 16 2m 3% sm 3% 4 4 


Profondeur] 


J_Æ | 


5 L 
+ Avant intervention 
— Après intervention avec résine Uratek 
— Après intervention avec résine moins expansive 


Figure 7.7 Comparaison entre les valeurs de Ms aux profondeurs d'injection, avant et après l'intervention, 
avec les deux types de résines, 
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3. Injections de résines expansives 
dans les sols à grains fins 


On souhaite concevoir une intervention consistant à injecter des résines expansives 
Uretck, dans un sable fin, à une profondeur z = 3,5 m par rapport au niveau du terrain 
naturel, dans le but de compenser laffaissement d'une fondation carrée, dont les côtés 
L= 1,5 m, soumise à une charge 4 = 150 kPa. La nappe phréatique se trouve à une pro- 
fondeur de 0,5 m par rapport au niveau de la surface naturelle du sol. 

Les paramètres physiques et mécaniques du sable en question sone les suivants: 

+ poids volumique du sol, y = 18,5 kN/m?; 

+ module élastique de Young, E°= 1,17 MPa; Æ = 1,4 MPa; 

+ coefficient de Poisson, v’ = 0,25; v, = 0,5; 

+ angle de frottement, @ = 25°; 

+ degré de surconsolidation, OCR = 1. 


Étant donné que le sol est normalement consolidé, le coefficient de pression des terres au 
repos, &,, peut être caleulé comme suit: K, = 1—sin(q), ce qui donne #, = 0,577. 
Pour évaluer l'affaissement, #, de la fondation, il est possible de se servir de l'expression 
suivante 


[74] 


dans laquelle 8, représente la base équivalente de la fondation, qui est égale au diamètre 
d’une fondation circulaire ayant une surface équivalente et, par conséquent, peut être cal- 
culée de la manière suivante: 


Vr 
Dans le cas présent, l'faissement west de 20,3 em. 


La variation de l'état de contraintes agissant à la profondeur z peut être déterminée à l'aide 
des formules suivantes: 


[7.6] 


(771 
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z 


ps] 


dans lesquelles AO, A0 , et A, sont respectivement les différences de contrainte verticale, 
de contrainte horizontal et de contrainte moyenne à la profondeur z. 

Pour le cas présent, la variation de contrainte verticale à la profondeur de 3,5 m est de 
12,3 kPa, tandis que la variation de contrainte horizontale est de —0,9 kPa. Par consé- 
quent, la contrainte verticale totale 0, est de 77,0 kPa, tandis que la contrainte horizon- 
tale totale 0, est de 43,6 kPa. 

On veut tout d'abord injecter une masse de résine M, = 10 kg. Sachant que la densité 
initiale de la résine est égale à 1 070 kg/m”, le volume initial de résine est V,,= 9,35: 10 m3. 
Comme prévu, puisque le coefficient de pression des terres au repos est inférieur à 1, on 
peut s'attendre, tout d'abord, à là formation d'une fracture sub-verticale en forme de 
pièce de monnaie (penmy-shaped crack). En référence au mode opératoire décrit dans Le 
chapitre 4, en supposant une demi-longueur de la fracture £ = 0,77 m, on atteint une 
pression d'équilibre de l'ordre de 123,5 kPa, comme on peut le voir clairement sur le gra- 
phique de la figure 7.8. En outre, dans l'hypothèse d'une ténacité à la rupture 
,08 MPa: re la pression critique, associée à la demi-longueur de la fracture, est 


17.81 
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Le on n 
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Figure 7.8 Rapport de gonflement de la résine en fonction de sa pression de gonflement ligne rouge) 
et dela contrainte de confinement du sol (ligne bleue); évaluation de la pression d'équilibre P* 

à une profondeur de 3,5 m par rapport au niveau du terrain néturel, dans l'hypothèse 

d'une demiongueur de fracture de 0,77 m. 
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Bremples de calcul 


On peut donc en conclure que pour le sable considéré, en supposant une ténacité à la 
facture K, = 0,08 MPa Tr demilongueur L = 0,77 m correspond à la valeur 
exacte, qui donne la congruence entre la pression critique dérivée du critère d’Invin ex la 
pression d'équilibre. 

La formation de La fracture sub-verticale suivie de son ouverture produit localement une 
augmentation de la contrainte horizontale du sol, égale à la pression d'équilibre nette, et 
done une inversion des contraintes principales: la pression d'équilibre nette dans le cas 
présent (cas étudié) est égale à la différence entre La pression critique, évaluée par le critère 
d'Irwin, et la contrainte horizontale agissant à la profondeur examinée: elle est de 80,8 kPa. 


Une nouvelle injection de 10 kg de résine peut conduire à la formation d’une fracture 
sub-horizontale, dont l'ouverture permet de soulever la fondation et de compenser une 
parie de l'affaissement. 

Les prévisions concernant la demi-longueur de la fracture horizontale peuvent être 
conduites avec le même mode opératoire que celui décrit précédemment, ou en recher- 
chant l'état d'équilibre entre la pression de gonflement de la résine et la contrainte de 
confinement (dans ce cas, verticale) générée dans le sol: dans ce cas, en supposant une 
demi-longueur de la fracture de 0,72 m, on atteint une pression d'équilibre égale à 160 kPa, 
comme on peut le voir clairement sur le graphique de la figure 7.9. 


v 


—Prasion de gontement 
— Conte da confinement 


PRemi mab dl ka bkpab rt ob, 


Vo 25 50 75 100 125 160 VS 200 225 260 275 500 325 59 
bond k ab dl kabi bkpab rt ob bp KR ab 1! Kaki bp 
sb tb md dl gor nd db, / (F/RI 


Figure 7.9 Rapport de ganflement de la résine en fonction de sa pression de gonflement ligne rouge) 
et de la contrainte de confinement du so (igne bleue); évaluation de La pression d'équilibre P* 

à une profondeur de 3,5 m par rapport au niveau du errain naturel dans l'hypothèse 

d'une demiongueur de fracture de 0,72 m. 


L'ampleur du soulèvement, 1, peur être évaluée, dans un premier temps, en négligeant la 
compression de la couche de sol au-dessus du point d'injection, au moyen de la formule 
fournie par Sneddon & Lowengrub (1969) pour le calcul de l'ouverture de la fracture : 

4-2 VE 2 
ER —— [7.10] 


Le soulèvement maximal obtenu au niveau de l'axe du trou d'injection (à savoir: pour une 
distance x = 0 ex pour le cas en question) est de 7,3 em. 


ge A 


Ÿ Théorie de l'expansion 
d'une cavité dans le sol 


La théorie de l'expansion des cavités est une référence fondamentale pour la conception 
des interventions nécessitant l'injection de résines polyuréthanes expansives dans Le sol. 
Des solutions pour matériaux élastiques linéaires, parfaitement plastiques, ont été présen- 
tées par Bishop e al (1945) pour les sols purement cohésifs, et par Hill (1950), Ménard 
(1957), Chadwick (1959), Vesic (1972), Carter er al. (1986), Manassero (1989), Ya 8e 
Houkby (1991) pour Les sols froutants et cohésifs. 


Introduction et définition du problème 


Cerre théorie analyse en détail le problème lié à l'expansion d'une caviré cylindrique ou 
sphérique dans un milieu infini, caractérisé par une loi de comportement de type élasto- 
plastique parfait et d’un critère de rupture de Mohr-Coulomb avec une loi d'écoulement 
plastique non associé. 

Dans le cas partieulier d’une cavité cylindrique ('est-à-dire des traitements en colonne), 
il est possible de se référer au schéma de La figure A.L, qui illustre un milieu étendu à l'in- 
fini, traversé par une cavité cylindrique de rayon 4, dans un système de coordonnées de 
référence de type polaire, dont l'origine coïncide avec le centre du trou. 


: 


Figure A.1 Problème d'expansion d'une cavité cindrique dans un mieu élasto-plastique parait 
infiniment étendu. 


Les hypothèses iniiales prévoient des conditions de déformations planes (déformations 
nulles dans la direction parallèle à l'axe du trou (axe des z)). 

Compte tenu du milieu infiniment étendu er homogène, l'état de contraintes aux alen- 
tours de la cavité avant l'expansion peut être considéré comme étant isotrope et calculé 
comme suit 


+2 
Reg ma ne, =0et2 00 au 
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+0, = contrainte agissant dans la direction radiale par rapport à l'axe de la cavité; 
+ op = contrainte tangentielle; 

+0, = contrainte efficace agissant dans la direction verticale; 

+ 6,9 = contrainte géostatique verticale; 

horizontale. 


+ Gp = contrainte géost 
Dans l'analyse, on supposera que la pression initiale agissant sur les parois de la cavité est 
égale à P et qu'elle augmente progressivement par la suite jusqu'à atteindre une valeur 
générique P. Les conditions de symétrie du schéma de caleul adopté impliquent que certe 
augmentation de pression conduit à des déplacements uniquement dans le sens radial 
(u) 

Dans la phase initiale, l'augmentation de pression dans la cavité provoque des déforma- 
tions du sol environnant dans des conditions purement élastiques. À une valeur donnée 
de P (P,), le matériau adossé aux parois de la cavité atteint les conditions de plastification, 
qui s'étendent progressivement dans le sol environnant au fur et à mesure de l'augmenta- 
tion de P (figure A.2). 

Les paragraphes qui suivent présentent de manière détaillée les procédures de calcul per- 
mettant d'établir, dans des conditions élastiques et dans des conditions de plastification, 
l'état de contraintes et de déformations du matériau entourant la cavité en expansion. 


Pesenrsiei 
aime ns © ñ 


Figure A2 Conditions de plastfication aux alentours de la cavité dues à l'augmentation de P. 


Il convient de souligner que toutes les corrélations rapportées ci-dessous sont obtenues en 
se référant à des analyses en conditions drainées, dans l'hypothèse où le temps nécessaire 
à la consolidation du sol enourant la cavité est négligeable comparativement à celui 
nécessaire à l'application des charges. Par conséquent, les contraintes indiquées doivent 
être comprises comme étant efficaces. En particulier, la pression P'à l'intérieur de la cavité 
doit également être comprise comme étant efficace. Dans le cas où l'on doit cffectuer une 
analyse dans des conditions non drainées, le sol peut être assimilé à un milieu à phase 
unique. Dans ce cas, l'analyse doit se rapporter uniquement aux contraintes totales, en 
adoptant, dans le domaine élastique, le module de rigidité non draîné et une valeur de 
coefficient de Poisson égale à 0,5, er, dans Le domaine plastique, un critère de rupture 
purement cohésif. Dans les formules ci-dessous, pour le milieu cohésif-frottant, il est pos- 
sible de poser ç = 1 pour éviter des problèmes d'ordre numérique. 
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État de contraintes et de déformations dans des conditions purement élastiques 


Enfin, il convient de souligner que les corrélations présentées ci-dessous font référence 
aux conventions de signes propres à la géorechnique classique, de sorte que Les sollicita- 
tions qui induisent une compression de l'élément de sol examiné sont considérées comme 
constituant un signe positif. 


2. État de contraintes et de déformations dans 
des conditions purement élastiques 


2.1. Cavité cylindrique 


Dans un premier temps, l'augmentation de pression à l'intérieur de La cavité cylindrique 
détermine l'apparition de déformations dans le sol environnant de type purement élas- 
tique. 

Les conditions de symétrie radiale du schéma de caleul adopté permettent d'expliquer 
l'érar de contraintes de l'élément générique du sol au moyen d'une seule équation d'équi- 


libre: 


do 
eu [A2] 
En outre, sont également supposées les conditions aux alentours suivantes: 
+ limite d'étar valide le long des parois latérales de la cavité cylindrique 
9,(8) =P [A3] 
+ condition limite valide à une distance infinie de la cavité cylindrique 
5,(-)=0(+)=8 IAA] 


En ce qui concerne l'érar des déformations, il est possible d'expliquer les égalités sui- 
vantes: 


[A5] 


+, = déformation dans la direction radiale par rapport au centre de la cavité cylin- 
drique; 

+ 6o = déformation tangentielle. 

Enfin, compte tenu des conditions purement élastiques, la corrélation entre l'état de 

contraintes et l'état de déformations est donnée par la loi de Hooke: 


0,=A+2-GE+P, AG] 
= Ae+2 GE +2, [A7I 
dans laquelle: 
+ A = opérateur de Lamé 
pa Ne [AS] 
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module de Young du sol; 


coefficient de Poisson du sol; 


déformation volumétrique 

+ G= module d'élasticité tangentielle du sol. 

En combinant l'équation d'équilibre [A.2] avec l'équation [A.G], et en tenant compte de 
la corrélation entre les déformations et les déplacements [A.5], on obrient l'équation dif 
férentiele du second ordre ci-après, exprimée en fonction du déplacement dans la dircc- 
tion radiale 1: 


[A9] 
à partir de laquelle on obtient: 
G, 
a, = Cort IA.10] 
r 
avec: 
+ G, C; = constantes d'intégration. 
En imposant à [A.10] la condition aux limites [A.4], on obtient : 
G=0 IA] 
à partir de laquelle : 
G [A2] 


En substituant [A.12] dans [A.5], il est possible d'obtenir la déformation radiale €, à 
inclure dans l'équation [A.6]. De cette façon, on obtient une expression de la contrainte 
radiale 5, qui, en imposant la condition aux limites [A.3], permet de déterminer la valeur 
de la constante d'intégration G. 

En conclusion, à la lumière de ce qui précède, l'état de contraintes et de déformations du 
matériau entourant la cavité cylindrique dans des conditions purement élastiques peut 
être défini au moyen des équations suivantes : 


(=) à 

pe 1A13] 

(ne IA 14] 
2-G-r 

Aa) 

tAïGI 


Pour caleuler l'augmentation de contrainte dans la direction verticale (A0) induite par 
l'expansion de La cavité, on peut se référer à la relation dérivée de la théorie de l'élasticiré 


Ao, =v-(Ao, +A0,) AA7] 


CP 
ai 
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La relation [A.17] est valide dans l'hypothèse selon laquelle les déformations dans la direc- 
don verticale sont nulles, et la contrainte principale verticale G, a une valeur comprise 
entre celles de la contrainte radiale principale 0, et celles de la contrainte rangentielle 
principale o, (0, > 9, > y). Les deux conditions sont généralement remplies si l'on 
adopte des valeurs sypes des paramètres géotechniques caractéristiques des sols naturels. 


2.2. Cavité cylindrique ou sphérique 

Yu & Houlsby (1991) ont étendu les propos du paragraphe précédent, également au cas 
de la cavité sphérique. 

En particulier, en adoptant un schéma de référence analogue à celui présenté dans la 
figure AL, les équations suivantes ont été élaborées : 


@ D'OR IA19) 


[A.20] 


+ _m= valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le cas d'une cavité sphérique ; 

+ 49 = rayon initial de la cavité ; 

G= module d'élastciré rangentielle du sol. 


Même dans ce eas, le caleul de l'augmentation de la contrainte dans la direction verticale 
(AG induite par l'expansion de La cavité peut être effectué en référence à [A.17]. 


3. État de contraintes et de déformations 
dans des conditions élasto-plastiques 


Dans le cas où la pression à l'intérieur de la cavité (7) est uelle qu'elle est égale à la valeur 
de la « pression d'avant la plastification » (P)), le matériau formant les parois de cette 
cavité atteint les conditions de rupture. 

2) Critère de rupture de Mohr-Coulomb pour un matériau purement frottant 


f=5,-N:0,=0 [A21] 
avec 
1+sinÿ aa 
1-siné 


dans laquelle : 
— 4 = angle de résistance au cisaillement du sol 
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) Critère de rupture de Mobr-Coulomb pour un matériau frottant et cohéif 


f=0,-N:oœ-Y=0 [A23] 
— N= voir [A.22] ; 
Eu à 1A24] 
1= sing 
dans laquell 


= 


angle de résistance au cisaillement du sol ; 
-ohésion effective du sol. 


Dans le cas où il se produit une augmentation ultérieure de P (P> P,), on observe la for- 
marion d'une « région de plastfication », qui s'étend depuis le rayon de la cavité (a) 
jusqu'à une distance générique (a < r < 6), en fonerion de l'ampleur de ? (figure A2). 
En dehors de cetre région (r > 8), le matériau continue à avoir un comportement pure. 
ment élastique. 

À l'intérieur de la région de plastification, en revanche, le matériau se trouve dans le 
domaine plastique. Par conséquent, son état de contraintes-déformations est subordonné 
au critère de rupture et à la loi d'écoulement plastique qui le caractérisent. 

Il est important de souligner que la pression P ne peut augmenter que jusqu'à une valeur 
limite (2), au niveau de laquelle la cavité tend idéalement à s'élargir indéfiniment dans 
le sol. 

Les paragraphes suivants décrivent de manière détaillée l'analyse du comportement du 
matériau après que la pression d'avant la plastification est arteinte (P;) en se référant au 
domaine élastique, et en se référant au domaine plastique, jusqu'à atteindre la Pins 

Il convient de souligner que ce problème peut être résolu selon deux approches distinetes 
1. l'approche par «analyse du sol aux petites déformations; 

2. l'approche par «analyse du sol aux grandes déformations». 

L'approche 1 repose sur l'hypothèse selon laquelle, une fois P) atteinte, les déformations 
provoquées par l'augmentation progressive cominue de P sont de faible importance, non 
seulement dans le domaine purement élastique, mais également dans la région de plasti- 
cation. Cette hypothèse conduit à une simplification considérable du problème; mais, 
dans le cas Le plus fréquent où les déformations dans le domaine plastique ne sont pas 
pertes, elle se traduit par une moindre précision des résuleats obrenus, notamment pour 
le caleul de la pression limite (P 

L'approche 2, en revanche, repose sur l'hypothèse que les déformations affectant le maté- 
riau dans la région de plastification sont tellement élevées qu'il faut les corréler aux dépla- 
cements correspondants au moyen de relations de type logarithmique: 


[A5] 


[A.26] 
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dans laquelle: 
distance générique depuis le centre de la cavité après le processus d'expansion ; 
distance générique depuis le centre de la cavité initiale avant le processus d'expan- 


# 


3.1. Estimation de la pression d'avant la plastification 


À la lumière de ce qui précède, il est nécessaire, en premier lieu, d'estimer la pression à 
l'intérieur de la cavité qui est en mesure de provoquer la plastification du matériau consti- 
tuant ses parois (pression d'avant la plastification (P,)). 


3.1.1. Cavité cylindrique, sol purement frottant 


Dans le cas particulier d'une cavité cylindrique et d'un sol purement frottant, la pression 
d'avant la plastfication (2) peut être caleulée en appliquant les conditions de rupture 
[A21] aux contraintes calculées avec [A.15], [A.16] au niveau des parois de la cavité 
(= a). 


De cette façon, on obtient l'égalité suivante: 


2N 
P=P, sk [A.27 
: (2) {1 +sing) A 271 


+ N= voir [A.22]. 


3.1.2. Cavité cylindrique ou sphérique, sol frottant et cohésif 


Afin d'étendre les analyses présentées dans le paragraphe précédent à la cavité sphérique 
er aux sols caractérisés par un comportement frotrant et coésif, Yu & Houlsby (1991) 
ont proposé La relation ci-après pour calculer la pression d'avant La plastification (P): 
mLY+(N -1).2 
TER 
MEN 


(A28] 


+ m = valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le cas d'une cavité sphérique; 

+ N= voir [A2]; 

+ Ÿ= voir [A24] 


3.2. Approche par «analyse du sol aux petites déformations » 


3.2.1. État de contraintes-déformations du sol dans le domaine élastique 


Le présent paragraphe présente une analyse de la définition de l’état de contraintes-défor- 
mations du matériau entourant la cavité en expansion dans le domaine élastique (r > 4) 
après que la pression d'avant la plastification a été atteinte (P). 
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321.1. Cavité cylindrique, sol purement frottant 


Dans le cas particulier d'une cavité cylindrique et d'un sol purement frotrant, Davis & 
Selvadurai (2002) ont mis au point une procédure d'évaluation de l'état de contraintes 
déformations du matériau entourant la cavité en expansion dans le domaine élastique 
(r > b), après que la pression d'avant la plastification a été atteinte (2,). 

En premier lieu, il est possible d'affirmer qu'à l'interface entre le domaine purement élas- 
tique et la région de plastification (r = 4), la contrainte radiale peut être supposée égale à 
la pression d'avant la plastification [A.27]: 


5,(6)= R= R-(1+sint) [A.29] 


En outre, comme il s'agit du domaine élastique, il est possible de considérer comme 
valides Les relations [A.13] à [A.16], qui, en remplaçant les termes à et P respectivement 
par b et 0,(b), permettent de définir les déplacements, les déformations et les contraintes 
du matériau qui se trouve dans le domaine élastique (r > b): 


IA.30] 
[A1] 


IA 32] 


IA33] 


Le caleul de l'angmentation de la contrainte dans Ia direction verticale (AG) induire par 
l'expansion de la cavité peut être effectué en référence à [A.171. 


3.2.2. Cavité cylindrique ou sphérique, sol frottant et cohésif 
Yu & Houlsby (1991) ont étendu les analyses effectuées dans le paragraphe précédent, 


également au cas d'une cavité sphérique, en examinant un matériau de cype frottant et 
cohésif entourant cette cavité. 

À la lumière de ce qui précède, les déformations et les contraintes du matériau qui est 
situé dans le domaine élastique (r > D), peuvent être calculées en utilisant les équations 
suivantes: 


(Y+(N-1) 2 {ef 


#7 (rune &) Fe 
LEHN=DRE pul csv 

o, ns d- [A.351 
Cp db T4 (Ni) IE 

nr Ne # L IA.36] 


où: 


+ m = valeur numérique étant supposée é 


le à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le cas d'une cavité sphérique; 
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+ N= voir (A.22]; 

+ V= voir (A24]. 

Le caleul de l'augmentation de la contrainte dans la direction verticale (AG) induire par 
l'expansion de la cavité peut être effectué en référence à [A.17]. 


3.2.2. État de contraintes-déformations du sol dans la région de plastification 


Le présent paragraphe expose une analyse de la définition de l'état de contraintes-défor- 
mations du matériau entourant la cavité en expansion dans la région de plastfication 
(a < r < D) qui s'est formée après que la pression d'avant la plastifcation (P;) a êté 
atteinte et dépassée, 


3.22.1. Cavité cylindrique, sol purement frottant 


Dans le cas particulier d'une cavité cylindrique et d'un sol purement frottant, Davis & 
Selvadurai (2002) ont mis au point une procédure d'évaluation de l'état de contraintes 
déformations du matériau entourant la cavité en expansion dans la région de plastifica- 
don (a < » < b), en référence à l'approche simplifiée par «analyse du sol aux petites 
déformations». 

Dans des conditions de plastificarion, l'état de déformations du sol dépend de la loi 
d'écoulement plastique qui caractérise le sol proprement dit. En termes généraux, étant 
donné une loi d'écoulement plastique non associé (matériau se dilatant), on peut appli- 
quer l'égalité suivante: 


g=5,-M:0, IA37] 
où: 
+sinp 
L A38 
1=sin me 


dans laquelle: 
+= angle de dilarance du sol. 

En imposant les conditions de plastification [A.21] à [A.2], on peut définir l'équation 
d'équilibre valable dans La région de plastification:: 


(N= 
do,, 0, (N 20 


[A.39] 


où: 
+ N= voir [A2]. 

La solution de l'équation [A.39], au vu de ce qui est reporté dans [A.32], [A.33] dans l'hy- 
porhèse où »= 4, permet de définir l'écat de contraintes du sol entourant la cavité cylin- 
drique dans la région de plastification 


Des 

a,=P,-(+ sin0) ( LA 40] 
+ 
Des 


o=2-(1- no?) IA41] 
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dans laquelle: 

+ N= voir (A22]. 

Le caleul de l'augmentation de la contrainte dans la direction verticale (AG) induite par 
l'expansion de la cavité peut être effectué en référence à [A.17]. 

En ce qui concerne l'état de déformations dans la région de plastification, il est nécessaire 
de souligner que, dans l'hypothèse où Le sol entourant la cavité 6 par une loi 
d'écoulement plastique non associé [A.37], la relation reliant les contraintes aux déforma- 
tions est la suivante: 


£ caractéri 


ë : [mé 
EE LÉ IA 42] 
SF C-VOAN- VEN ln MN |e | 
+ G= module d'élaricité tangentelle du sol; 
+ = coefficient de Poisson du sol; 
+ N= voir [A2]; 
+ M= voir [A.38]. 
Le système [À.42] comprend une équation unique linéairement indépendante, qui, expli- 
citée en fonction de l'inconnu 6, et à la lumière de [A.S], peut être exprimée comme 
suit: 


di, à _(=V (+M-N)-V{M+N) 6, 


M. [A.43] 
der N 2:G L 1 
En dérivant [A.40] par rapport à la variable 6, on obtiei 
CE Gran) (® Je LA 44] 
N F3 


M: 


Cr Cé-v}G+men)-v(+ N)] 3e (£ 
dr G N 


r 


dont la solution, calculée en imposant la condition aux limites dérivée de [A.30]: 
Prsing-b 


u,(b)= [A.46] 
(= 146 
permet de caleuler le déplacement dans la direction radiale dans la région de plastifica- 
tion 
ï fr  peÿ au fre y] 
OP OOIE  DIEI 
2:G 6) N-1/&% b M+1/\6 r 


dans laquelle: 

3-2 [u-v) (+) -v (m4 n)] 
MEN 

A=2+8 IA.49] 


IA 48] 
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Yu & Houlsby (1991), en référence aux résultats des travaux de Carter er al. (1986), ont 
étendu les analyses effectuées dans la section précédente, également au cas de la cavité 
sphérique, en examinant un matériau de type frottant et cohésif entourant cette cavité. 


Cavité cylindrique ou sphérique, sol frottant et cohésif 


Dans le cas présent, l'éat de contraintes du matériau dans la région de plastification peut 
être calculé au moyen des équations suivantes : 


tm} N[Y+(N-}2] M 


F [A.50 

CN 1) + AN) ë pe 
Gt) {Y4N-0-2], 

CN =1) m4 N) AS 


N 

+ m = valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le cas d’une cavité sphérique ; 

+ N= voir (A2); 

# = voir [A24] 

Le caleul de l'augmentation de la contrainte dans la direction verticale (AG) induite par 

l'expansion de la cavité peut être effectué en référence à [A.17]. 

En ce qui concerne l'estimation du déplacement radial dans la région de plastification, il 

est possible de se référer à La relation 


Y+(N-1)-2 NME Min] (of. 
LE EN) G UNE MEN Mn) D) UE 
Min NME DRE 

TER res : 
e ; F HN -1)2 
1 ve (2m) | (4m) 
DrNNE., IA.S3] 
Cv) (W-1) M 
4 M+tm CRE) 
: EH As 4 à 


+ m = valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le eas d’une cavité sphérique; 

2 V= voir [A 24]: 

+ N= voir [A22]; 

+ M= voir [A8]; 


& 


141 


À 


Théorie de l'expansion d'une cavité dans le sol 


+ E= module de Young du sol; 
+ G= module d'élasticité rangentielle du sol; 
= cocfficient de Poisson du sol. 


3.2.3. Calcul de l'extension de la région de plastification 


Afin de compléter les analyses effectuées dans les paragraphes précédents pour l'évaluation 
de l'érat de contraintes ec de déformations du matériau qui se trouve dans la région de 
plastificatien, on trouvera ci-après un certain nombre de relations destinées aux estima- 
tions du paramètre b. 


323.1. Cavité cylindrique, sol purement frottant 


Si l'on analyse le cas de l'expansion d'une cavité cylindrique dans un milieu purement 
frottant, la valeur de 4 peut être calculée au moyen d'une relation obtenue à partir de 
AO, en faisant l'hypothèse qu'au niveau des parois de La cavité (r = 4), la contrainte 
radiale (a) est égale à la pression P: 


P=R {1+sind) 1.55] 


à partir de laquelle on obtient: 
[A:S6] 
avec: 


bé P I 
ÊE rl | 
+ N= voir [A.22); 


+ = angle de résistance au cisaillement du sol. 


3.23.2. Cavité cylindrique ou sphérique, sol frottant et cohésif 


Ya 8 Houlsby (1991) ont étendu les analyses effecruées dans le paragraphe précédent, 
également au cas d’une cavité sphérique, en examinant un matériau de type frottant et 
cohésif entourant cette cavité. 


Plus précisément, le paramètre 6 peut être estimé au moyen de l'équation suivante: 


IA57] 


Ont N)[Y+(N-0:2] 
Nm) [HN -)27] 
dans laquelle: 
+ m = valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylin 
égale à 2 dans le cas d'une cavité sphérique; 
+ = voir [A 24]; 
2 N= voir [A.22]. 
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3.3. Approche par «analyse du sol aux grandes déformations» 


Dans le cas où il s'avère impossible de simplifier le problème d'expansion de la cavité à 
l'aide de l'approche par l'analyse aux petites déformations, Yu & Houlsby (1991) ont mis 
au point une série de corrélations qui permettent de définir l'état de contraintes-déforma- 
ions du matériau entourant la cavité en présence de «grandes déformations» dans la 
région de plastification. 


L'un des aspects qui différencie surtout certe approche par rappor à celle trairée dans le 
chapitre précédent réside dans l'hypothèse selon laquelle, lorsque la pression d'avant la 
plastification (P;) est dépassée, l'augmentation de ? détermine une augmentation minime 
non négligeable de la dimension de la cavité, dont le rayon passe d'une valeur initiale 4 
à une valeur générique 4. 

Les paragraphes ci-après présentent les relations fondamentales permettant de calculer les 
contraintes et les déformations dans le sol entourant la cavité après atteinte et dépasse- 
ment de la pression d'avant la plastification (7), en référence au cas générique d'une 
cavité cylindrique ou sphérique er d'un matériau caractérisé par un comportement frot- 
tant et cohésif. 


3.3.1. État de contraintes-déformations du sol dans le domaine élastique 


L'état de contraintes-déformations du matériau entourant la cavité en expansion dans le 
domaine élastique (r > 4) après que la pression d'avant la plastification a été atteinte (7), 
peut être évalué en adoptant Les relations [A.34], [A.35], [A.36] et [A.17]. 

Dans ce cas spécifique, le seul changement par rapport aux informations figurant dans le 
paragraphe «Approche par «analyse du sol aux petites déformations »» réside dans l'éva- 
luation de l'extension de la région de plastification (paramètre 6) (paragraphe « Calcul de 
l'extension de la région de plastifications) 


3.3.2. État de contraintes-déformations du sol dans la région de plastification 


Après atteinte et dépassement de la pression d'avant la plastifcation (P)), l'état de 
contraintes du matériau entourant la cavité dans la région de plastfication (a < » < 4) 
peut être calculé au moyen des relations [A.50), [A.51] et [A.17]. 

Dans ce cas aussi, la différence par rapport aux informations figurant dans le paragraphe 
«Approche par «analyse du sol aux petites déformations »» réside dans l'évaluation de 
l'extension de la région de plastification (paramètre #) (paragraphe « Calcul de l'extension 
de la région de plastification») 


En ce qui concerne l'état de déformations du sol dans la région de plastification, il est pos- 
sible de se référer à l'égalité suivante: 


[A.59] 
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dans laquelle: 


ST) 

al 

N 

+ m = valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans Le cas d’une cavité sphérique; 


+ Y= voir [A24]; 
+ N= voir (A2); 


+ M= voir [A 38]; 
+ n= voir [A.54]; 
+ É= voir [A53]: 


+ G= module d'élasticité tangentielle du sol; 
+ R= voir [A.S8] 

Il est possible de démontrer que la série numérique [A.60] converge rapidement pour 
route valeur de Ver Matrribuée au sol. En règle générale, on peut affirmer que es dix pre- 
mières itérations sont suffisantes pour garantir un résultat de précision acceptable aux fins 
de l'analyse 

Après avoir calculé au moyen de (À.59] la valeur de 4, on peut estimer le déplacement 
radial et la déformation tangentielle au niveau des parois de la cavité: 


u(&)=2-4 [AG2] 
f 
&(&)= (2) [AG3 


où: 
+ 4 = rayon initial de la cavité. 
3.3.3. Calcul de l'extension de la région de plastification 


Après avoir caleulé le rayon final supposé de la cavité (a) à la suite de l'expansion détermi- 
née par la pression P (À.59], il est possible de calculer Le rayon de la région de plastifica- 
ion (4) au moyen de l'équation: 


IAA] 


+ m = valeur numérique érant supposée égale à 1 dans le cas d'une cavité cylindrique, 
le à 2 dans le cas d'une cavité sphérique; 

+ N2= voir (A2); 

+ R2= voir [A.S8]. 
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3.34. Calcul de la pression limite (Pa) 


Au moyen de l'approche aux grandes déformations, il est possible d'obtenir une relation 
qui permette de caleuler la valeur de la pression limite (A). 


Ainsi qu'il est indiqué dans l'introduction du paragraphe «État de contraintes er de défor- 
mations dans des conditions élasto-plastiques », il s'agit de I limite au-delà de laquelle la 
pression à l'intérieur de la cavité ne peut plus subir d'augmentation, dans la mesure où, 
lors de l'apparition de cerre condition, la région de plastfication tend à se dilater indéf- 
niment dans le matériau entourant la cavité. 


En particulier, en appliquant à [A.59] la condition: 


li {:) IAGS] 
on obtient La relation: 

Mn 
nr (N 1) (re a 


N'(M+m) 2(m+ N).G IAE] 


at) 


+ m = valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le cas d'une cavité sphérique; 

+ N= voir [A2]; 

2 = voir [A2A]: 

+ M= voir [A.38]; 

en = voir [453]; 

= voir [A.52]; 

+ G= module d'élasticité tangentielle du sol. 


obtenue en posant le dénominateur du second membre de l'équation [A.58] égal à 0. 

Au moyen de [A.66], on obtient la valeur de A (R.,5). 

Ensuite, au moyen de [A.60] et [A.61], on obtient le paramètre R,, défini par la relation: 
Ont NY TY4(N 21) Pau] 
Nm) HN] 


Enfin, en inversant l'équation [A.67], on obtient la valeur de pression limite (7 19) 


IA.67] 
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Les mécanismes 


de fracturation 
dans le sol 


L'injection de résine expansive dans le sol peut provoquer la formation de fractures qui 
créent un élargissement considérable du rayon d'influence du traitement effectué, compa- 
rativement au rayon prévu selon la théorie d'expansion des cavités. 


Introduction et définition du problème 


Pour cette raison, les paragraphes suivants illustrent quelques-uns des aspects les plus 
importants liés à la mécanique des fractures appliquée aux sols. En particulier, il convient 
de souligner que cette chéorie permet de définir les conditions nécessaires à la formation 
de fractures dans le sol, en condition drainée et non drainée, après l'injection d’un fluide 
sous pression. 

Le problème lié à la propagation des fractures dans le sol est sans aucun doute bien plus 
complexe, et les seules analyses formulées à son sujet sont des analyses qualitatives déri- 
vécs d'études expérimentales sur site et en laboratoire. 


L2._ Fracturation du sol en condition non drainée 


La formation de fractures dans un sol dépend essentiellement de sa résistance à a raction 
Co). 

Dans les matériaux naturels, il est raisonnable de supposer que celle-ci a des valeurs très 
petites ou proches de 0 kPa, sauf en présence de processus de cimentation du squelette 
solide du sol. 


Plus précisément, il peut être observé que le sol se fracture lorsque la contrainte principale 
mineure efficace (0; est égale à la valeur de la résistance à la traction propre au matériau 
érudié: 

Matériau non cimenté 5,=0 (BA) 


Matériau cimenté 5=-0! 1B2] 


Il convient de souligner que les corrélations présentées font référence aux conventions de 
signes propres à la géorechnique classique, de sorte que les sollicitations qui induisent une 
compression de l'élément de sol examiné sont considérées comme un signe positif. En 
outre, dans [B.2], la résistance à la traction (0°) est supposée avoir un signe positif. 

Dans un milieu saturé er dans le cas général d'un sol cimenté, l'application d'une force de 
traction dans la direction dans laquelle agit la contrainte principale mineure détermine 
l'apparition d'un mécanisme de rupture par cisaillement dans des conditions non drai- 
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nées er la fracturation se produit au moment où [B.2] est satisfaite, laquelle peut être réé- 
rite sous la form 


0 [8.3] 
où: 
+05 = contrainte principale totale mineure; 
+ 4 = pression interstitell initiale dans des conditions hydrostatiques; 


+ À, = surpression intersttille générée lors du mécanisme de ruprure par cisaillement 
qui conduit à la rupture du sol. 


Les surpressions interstitielles peuvent être estimées au moyen de la relation proposée par 
Skempton (1954): 
Au 


[A0 +4 (A0, -A0;)] 18.4] 


dans laquelle: 
+ AG = variation de la contrainte principale totale majeure; 

+ As = variation de la contrainte principale totale mineure; 

+ A, B= paramètres de Skempton. 

Dans Le cas d'un sol saturé (8 = 1), en substituant la [B.4] à la [B.3] on obtient: 


5, = 4, +[A0,+ 4-(40, -A0,)]-0' 18,5] 


La [B.5] peut être exprimée comme sui 


(o,+a)+Aa (B.6] 


avec: 
+ 03 = contrainte principale efficace mineure avant l'application de la force de trac 
tion: 


ss 


— A0, # 18.7] 
Dans un sol cohésif (par ex., argile), l'injection de résines expansives Uretek peut provo- 
quer la Fracturation du matériau: 

Conformément à la théorie de l'expansion des cavités (annexe A), dans l'hypothèse où 
l'on schématise le comportement du sol avec une loi de comportement de type linéaire 
élastique-plastique parfait, i peut être démontré que, dans des conditions purement élas- 
tiques, l'augmentation de la pression agissant à l'intérieur de la cavité dérermine à La fois 
une augmentation de la contrainte radiale totale (AG) et une diminution de la contrainte 
tangentielle totale (AO4) (figure B.1(a)). 
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Figure B.1 Fracturation du sol déterminée par l'injection d'un fluide sous pression dans une cavité 


a) rupture en traction dans des conditions non drainées et () rupture due une résistance au cisaillement 
dans des conditions non drainées (d'après Mitchell & Sogs, 2005). 


Entre ces deux valeurs, la relation suivante est vérifiée: 


A6, [B.8] 
en posant: 
A5, [8.9] 
« 
49, =40, [8.10] 


dans des conditions de déformations planes (déformations nulles dans la direction verti- 
cale) (4 & 0,5), la [B.6] peut être réécrite comme suit: 

46,=6/,+6/ [Bu] 
Puisque AO, peur également être exprimée sous la forme: 


Ac,=P,-6 18.12] 


où: 

+ 05, = contrainte principale totale mineure avant l'augmentation de pression à l'inté- 
rieur de la cavité (égale à Py qui figure dans la théorie de l'expansion des cavités, 
annexe A); 

+ Py= pression à l'intérieur de la cavité qui induit La fracruration du sol. 
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En substituant la [B.12] dans la [B.11], on obient: 
Pr=2-0u M +0 18.13] 


Au moyen de la [B.13), il est possible d'estimer la pression à l'intérieur de la cavité qui 
provoque la fracturation du sol dans des conditions de déformations purement élastiques. 
Il est important de souligner que [B.13] a été dérivée en supposant une expansion uni- 
forme de la cavité au fur et à mesure de l'augmentation de la pression agissant sur les 
parois de cette cavité (théorie d'expansion des cavités), ainsi qu'une loi de comportement 
du sol de type linéaire élastique. En réalité, les déformations de la cavité ne sont pas uni- 
formes et a fracturation du sol peut apparaître dans des zones localisées, à une pression 
inférieure à celle calculée au moyen de [B.13]. 

En variante de ce qui précède, il peur arriver que l'augmentation de la pression à l'inté- 
rieur de la cavité provoque la réalisation des conditions de plastification au niveau des 
parois de la cavité proprement dite avant que la contrainte tangentielle efficace ne dimi- 
nue au point d'être égale à la résistance à la traction du matériau (figure B.1(b)). Cette 
condition est vérifiée sk 


CR [B.14] 
avec: 
+0, = contrainte radiale totale aux alentours de la cavité; 
+ Où = contrainte tangentielle totale aux alentours de La cavité; 
+ 6, = résistance au cisaillement en condition non drainée du matériau; 
er: 
5,>-06! 18.15] 


En référence à la théorie de l'expansion des cavités (annexe A), il est possible d'affir- 
mer que [B.14] est vérifiée au moment où la pression à l'intérieur de la cavité atteint la 
valeur de pression d'avant la plastification (?), qui, dans des conditions non drainées, 
peut être estimée comme suit: 


+6, BJ6] 


+ m = valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le cas d'une cavité sphérique; 

+ oj = contrainte principale totale mineure avant l'augmentation de pression à l'inté- 
rieur de la cavité (égale à 2, qui figure dans la théorie de l'expansion des cavités, 
annexe A); 


+ e, = résistance au cisaillement en condition non drainée du matériau. 
Dans le cas où il se produit une augmentation ultérieure de la pression à l'intérieur de la 
cavité (P > 2), on observe la formation d'une « région de plastification », qui s'étend 
depuis le rayon de la cavité jusqu'à une distance , en fonction de l'ampleur de P 
(figure B.2). 

En dehors de cette région (r> 9), le matériau continue à avoir un comportement pure- 
ment élastique. 
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À l'intérieur de la région de plastification, en revanche, le matériau se trouve dans le 
domaine plastique. Par conséquent, son état de contraintes-déformations est subordonné 
au critère de rupture et à la loi d'écoulement plastique qui le caractérisent. 


Prsentsei 
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Figure B.2 Conditions de plastication aux alentours de la cavité dues à l'augmentation de P. 


Dans la région de plastification, il est possible d'avoir, localement, une concentration des 
déformations de cisaillement et que, par conséquent, des microfissures apparaissent dans 
le matériau, Dans ces conditions, il est raisonnable de supposer que les injections de résine 
Urerek provoquent l'ouverture de microfractures et favorisent leur propagation dans le sol 
traité. 

À la lumière de ce qui précède, il est possible de conclure que [B.13] permet d'estimer la 
pression à l'intérieur de la cavité qui induit la fracruration du sol dans des conditions de 
déformations purement élastiques, tandis que, dans des conditions de plastification du 
matériau, il faudra se référer à [B.16]. 

Il est important de souligner que l'expansion de la région de plastification avant la fracru- 
ration du sol ou l'augmentation de la résistance au cisaillement, en condition non draînée 
du matériau, due à l'augmentation de la pression de confinement, impliquent une aug- 
mentation de la pression de fracturation (P) par rapport à celle calculée au moyen de 
18.16]. Ainsi, plusieurs auteurs (Jaworski ef 4l., 1981; Yanagisawa & Panah, 1989) ont 
proposé une série de corrélations empiriques en alternative à [B.16] pour calculer }, qui 
peuvent se résumer au moyen de l'équation suivante: 


Demo tn 18.17) 


dans laquelle: 

+ 03, = contrainte principale totale mineure avant l'augmentation de pression à l'inté- 
rieur de la cavité (égale à P, qui figure dans la théorie de l'expansion des cavités, 
annexe A); 

+ mn = constantes multiplicatives propres au sol traité, 


En ce qui concerne la plage de variations des valeurs de m, Jaworski ef al. (1981) ont rap- 
porté des valeurs comprises entre 1,5 et 1,8, tandis que pour Panah 8 Yanagisawa (1989), 
m = 1,05 + 1,085 m. 


= 
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La figure B.3 présence les valeurs de Pobtenues pour divers cypes de sols en fonction de 
la pression de confinement initiale (04). 

D'après le graphique, il ressort clairement que la corrélation qui lie les deux grandeurs 
pour les différents types de matériaux est toujours comprise entre [B.13] et [B.16], qui 
peuvent, par conséquent, être considérées comme une limite supérieure et une limite infé- 
rieure à utiliser pour le caleul de P; 
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Figure 8.3 Pression de fracturation (P) en fonction de La pression initiale de confinement 
pour différent types de sols (d'après Mitchell & Soga, 2005) 


3. Fracturation du sol en condition drainée 


En condition drainée, de façon similaire à ce qui a été indiqué dans le chapitre précédent, 
il peut être observé que le sol se fracture au moment où la contrainte principale efficace 
mineure (9',) est égale à la valeur de la résistance à la traction propre au matériau étudié 
{[B.1] et [B.2)). 

Conformément à la théorie de l'expansion des cavités (annexe A), dans l'hypothèse où 
l'on schématise le comportement du sol avec une loi de comportement de rype linéaire 
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élastique parfaitement plastique, il peut être démontré que, dans des conditions purement 
élastiques, l'augmentation de la pression agissant à l'intérieur de La cavité dérermine l'ap- 
parition d'une rension radiale et d'une tension tangentielle dans le matériau environnant, 
qui caïncident respectivement avec la tension principale majeure et mineure et peuvent 
être calculées comme suit : 


[B.18] 


[B.19] 


où: 

+ m = valeur numérique érant supposée égale à 1 dans le cas d'une cavité cylindrique, 
égale à 2 dans le cas d'une cavité sphérique; 

+ Ph = étar de contraintes efficace isorrope du sol autour de la caviré avant l'expansion: 


H=5= 18.20) 


dans laquelle 
A 
— d'y = contrainte géostatique horizontale efficace 
— My = état de contraintes total isotrope du sol autour de La cavité avant l'expansion: 
PE Pitu 821] 


ontrainte géostatique verticale efficace; 


ave 
= ny = pressions interstitielles initiales dans des conditions hydrostatiques; 

— P = pression à l'intérieur de la cavité qui détermine l'expansion de celle-ci. En 
condition drainée, dans l'hypothèse où le temps nécessaire à la consolidation du 
sol entourant la cavité est négligeable par rapport à celui nécessaire à l'applica- 
tion des charges, toute La différence de pression P— agit sur le squelette solide 
du sol, sans provoquer de surpressions interstitilles. Si ces hypothèses ne sont 
pas satisfaites, il est nécessaire de se référer à l'analyse dans des conditions non 
drainées ; 

— à 

— r = distance générique dans la direction radiale par rapport à la cavité. 


ayon initial de la cavités 


Au niveau des parois de la cavité (7 = &), en substituant [B.19] à la condition [B.2], on 
obtient: 

2,222 -u +0) 1822] 
+ Pj= pression qui détermine la pression de fracturation du sol entourant La cavité dans 

des conditions purement élastiques ; 

+ 07,= résistance à la traction du matériau. 
Mëme dans des conditions drainées, le sol entourant la cavité peut atteindre les condi- 
tions de plastfication avant de subir la fracturation dans le domaine purement élastique. 


LS 
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En référence à la théorie de l'expansion des cavités appliquée à un matériau cohésif'e frot- 
tant (annexe À), il est possible de caleuler la pression en mesure de provoquer la plas- 
rification du matériau à proximité des parois de la cavité au moyen de la relation: 


a OL EP 1823) 
m+N 
où: 

Lsind 1824 

1-sin 

2 2:c"cosû 
Fe (8.25) 

dans laquelle: 


= q = angle de résistance au cisaillement du sol, 


‘ohésion effective du sol; 

= valeur numérique étant supposée égale à 1 dans le cas d'une cavité cylindrique, 
le à 2 dans le cas d’une cavité sphérique; 

+ Py= voir [B21]. 


En conclusion, à l'instar de ce qui est rapporté dans l'annexe B, il est possible d'affir- 
mer que [B.22] et [B.23] permettent de définir la plage de variations de la pression agis- 
sant sur les parois de la cavité en expansion, pouvant provoquer la facturation du sol 
entourant le point d'injection, dans des conditions drainées. 
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